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SECTION  1 


INTRODUCTION 


This  report  describes  the  research  conducted  on  the  characterization 
of  fatigue  crack  growth  in  bonded  structures.  In  this  volume,  analytical 
techniques  required  to  predict  fatigue  crack  growth  in  bonded  structures 
are  discussed.  The  specific  problem  is  that  of  a two-ply,  bonded,  metallic 
laminate  with  a through-crack  in  one  ply.  The  analytical  techniques  developed 
are  applicable  to  adhesively  bonded,  weldbonded  and  brazed  structures. 

The  first  step  in  crack  growth  analysis  is  to  obtain  the  stress  inten- 
sity factors  for  the  structure  under  consideration.  Two  methods  for  obtain- 
ing stress  intensity  factors  for  cracked  bodies,  namely,  the  finite  element 
method  and  the  mathematical  method,  are  discussed  in  Section  2.  A two- 
dimensional  finite  element  model,  in  which  plates  are  represented  by  mem- 
brane elements  and  the  adhesive  by  shear  elements,  is  developed,  The  finite 
element  analysis  is  carried  out  with  NASTRAN.  In  the  mathematical  method 
discussed  in  this  volume,  the  problem  of  obtaining  stress  intensity  factors 
for  a cracked  body  is  reduced  to  the  solution  of  integral  equations,  using 
complex  variables,  and  these  integral  equations  are  then  solved  numerically. 

The  application  of  each  method  to  a variety  of  structural  problems  is  shown  * 

The  results  obtained  from  each  method  for  the  same  problem  are  compared. 

In  Section  3,  the  application  of  the  analyses  to  multilayered,  adhesively 
bonded,  we ldbonded , and  brazed  structures  is  shown. 

In  Section  4,  the  influence  of  bending  caused  by  the  presence  of  a crack 
in  only  one  layer  of  a multilayered  structure  is  discussed.  A new  method  to 
correct  the  stress  intensity  factors  for  the  influence  of  bending  is  described. 

In  Section  5,  a criterion  for  the  cracking  of  a sound  layer  bonded  to  a 
cracked  layer  is  developed.  The  criterion  is  based  on  the  load  transferred  to 
the  sound  layer  from  the  cracked  layer. 
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The  propagation  of  a crack  in  an  adhesively  bonded  structure  is 
accompanied  by  the  propagation  of  debond  in  the  adhesive.  In  Section  6,  a 
criterion  for  the  propagation  of  debond  is  developed. 

The  guidelines  for  the  application  of  analysis  to  bonded  structures 
are  outlined  in  Section  7* 

The  computer  programs  developed  in  this  research  program  are  outlined 
in  the  appendices.  The  computer  programs  developed  for  the  mathematical 
method  of  analysis  for  three  classes  of  problems  are  as  follows: 

1,  A cracked  sheet  stiffened  by  a partially  debonded  stringer  at 
an  arbitrary  location, 

2,  A cracked  sheet  with  two  partially  debonded  stringers,  symmet- 
rically located  about  the  centerline, 

3,  A two-layer,  adhesively  bonded  structure  with  a debonding  crack 
in  one  layer. 

These  computer  programs  require  the  elastic  properties  of  the  adhesive, 
cracked  layer,  and  sound  layer  as  inputs  and  have  provisions  for  account- 
ing for  the  debond  in  the  adhesive*  The  computer  run  times  for  these 
programs  are  very  short  compared  to  those  for  finite  element  analyses, 
hence  the  programs  are  well  suited  for  parametric  studies. 
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SECTION  2 


ANALYSIS  OF  CRACKED  ADHESIVELY  BONDED  STRUCTURES 


The  analysis  of  crack  problems  in  adhesively  bonded  structures  has 
attracted  attention  In  recent  years  due  to  the  increased  application  of 
bonding  in  primary  aerospace  structures.  However,  very  little  analytical 
work  has  been  done  in  this  field  due  to  limited  availability  of  analytical 
tools  to  handle  the  complicated  problems.  While  one  of  the  most  common  tools 
used  for  analysis  of  complex  structures  is  the  method  of  finite  elements,  an 
adhesively  bonded  structure  is  inherently  three-dimensional . A three- 
dimensional  finite  element  analysis  of  cracks  in  adhesively  bonded  structures 
does  not  seem  to  be  feasible  at  this  time,  mainly  because  of  a lack  of  suit- 
able three-dimensional  cracked  elements  to  analyze  such  problems.  However, 
a modified  two-dimensional  finite  element  analysis  was  used  in  Reference  1 
to  analyze  a cracked  sheet  with  a bonded  stringer*  A good  correlation  was 
obtained  between  analysis  and  experiments,  which  suggested  that  a modified 
two-dimensional  finite  element  model  can  be  used  for  analysis. 

Finite  element  methods  are  very  versatile  tools  for  analyses.  Various 
stiffener  geometries,  finite  width  effects,  influences  of  holes,  etc.,  can 
be  studied  by  this  method.  The  major  disadvantages  are:  considerable  time 
is  required  to  model  the  structure,  computer  run  times  are  often  long,  and 
the  solutions  are  not  closed  form,  thus  cannot  be  easily  generalized.  If 
finite  element  techniques  are  used  to  do  parametric  studies  of  any  nature, 
the  cost  would  be  prohibitive,  hence  this  technique  Is  not  well  suited  for 
any  parametric  stress  intensity  factor  studies. 

The  use  of  mathematical  methods  to  study  crack  problems  in  adhesively 
bonded  structures  has  been  reported  by  several  researchers.  References  2 
through  4 have  used  a complex  variable  formulation  to  reduce  the  problem  to 
the  solution  of  a set  of  integral  equations  that  are  solved  numerically.  The 
mathematical  techniques  are  applicable  to  bonded  plates  that  have  infinite 
width.  For  bonded  stringer  cases,  certain  simplified  assumptions  regarding 
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stiffener  areas  are  made.  The  advantages  of  the  mathematical  method  are  that 
the  solutions  are  obtained  by  using  numerical  analysis*  which  generally  requires 
short  computer  run  times*  and  parametric  variations  are  permitted,  as  well  as 
generalized  solutions#  The  computer  programs  developed  in  this  study  for 
mathematical  analysis  have  the  elastic  properties  of  adhesives,  adherends,  and 
debond  sizes,  etc*,  as  inputs,  hence  parametric  studies  can  be  easily  made  at 
a greatly  reduced  cost#  The  main  disadvantages  are  that  certain  simplifying 
assumptions  have  to  be  made  to  account  for  stiffener  area  or  finite  boundaries, 
hence  the  results  obtained  must  be  carefully  reviewed  for  such  cases* 

The  principles  involved  in  the  two  methods  of  anaLysis  and  their  applica- 
tion, as  well  as  the  assumptions  made,  are  discussed  in  the  following  paragraphs* 

2.1  FINITE  ELEMENT  METHOD  OF  ANALYSIS 

The  presence  of  singular  stresses  ahead  of  a crack  tip  require  special 
care  in  finite  element  analysis.  Special  elements,  known  as  cracked  elements, 
have  been  developed  to  account  for  the  singularity  ahead  of  a crack  tip  in  a 
two-dimensional  finite  element  analysis*  For  three-dimensional  crack  problems, 
the  cracked  elements  are  still  in  the  development  stage.  The  analysis  of  adhe- 
sively bonded  structures  with  a three-dimensional  cracked  element  presents  a 
very  complicated  problem  that  appears  to  be  beyond  the  scope  of  the  present 
state  of  knowledge.  Although  increasing  amounts  of  attention  are  being  given 
to  this  problem  (Reference  5).  Hence,  a modified  finite  element  analysis  using 
the  available  two-dimensional  cracked  element  has  been  adopted  in  the  present 
program*  This  method  of  modeling,  along  with  corrections  for  bending  (discussed 
in  Section  4)  has  shown  a very  good  correlation  between  analytical  and  experi- 
mental results. 

The  method  of  finite  element  analysis  uses  the  available  two-dimensional 
cracked  element  ahead  of  the  crack  tips  in  the  cracked  layer.  The  cracked  and 
sound  layers  are  modeled  as  two-dimensional  structures.  The  adhesive  layer  is 
treated  as  a shear  spring  rather  than  as  an  elastic  continuum.  This  method  of 
treating  adhesives  has  been  used  before  in  the  analysis  of  bonded  joints  and 
cracked,  adhesively  bonded  structures  (References  1-4*6),  In  the  finite  clement 
analysis,  the  continuous  shear  spring  is  represented  by  shear  elements  (Refer- 
ence l).  With  these  assumptions,  it  is  possible  to  model  the  structure  as 
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two-dimensional.  Consider  the  adhesively  bonded  structure  shown  in  Figure  1 
where  a through-crack  exists  in  the  plate,  and  a debond  exists  in  the  adhe- 
sive around  the  crack  in  Material  1,  however  there  is  no  crack  in  Material  2, 
Although  the  debond  may  have  any  shape,  it  is  assumed  to  be  elliptical  here, 
for  convenience. 

The  bonded  plates  of  both  materials  are  modeled  as  rectangular  or  triangu- 
lar elements  as  shown  in  Figure  2*  The  enlarged  portion  of  Figure  3 identifies 
the  location  of  a cracked  element  provided  ahead  of  the  crack  tip  in  the  cracked 
plate*  The  x and  y coordinates  of  the  grid  points  in  the  bonded  regions  of 
the  two  plates  are  kept  the  same,  i.e,,  the  finite  element  model  in  the  bonded 

regions  of  the  plates  is  identical*  In  the  unbonded  region,  the  finite  ele- 

ment model  of  the  two  plates  need  not  be  the  same  since  no  connection  is  pro- 
vided between  the  plates  in  the  unbonded  region.  Around  the  periphery  of  the 
debond  in  the  bonded  region,  the  adhesive  is  subjected  to  high  shear  stresses, 
hence  a closely  spaced  mesh  is  provided  to  reduce  the  length  of  the  shear  ele- 
ments shown  in  Figure  3.  The  closely  spaced  shear  elements  will  give  an  accu- 
rate determination  of  shear  stresses  and  load  transfer  to  the  uncracked  layer. 
Thus  providing  for  more  accurately  computed  stress  intensity  factors. 

In  the  bonded  region,  the  grid  points  of  the  two  materials  are  connected 
by  the  shear  elements*  These  shear  elements  have  the  same  material  proper- 
ties as  the  adhesive.  The  thickness  of  the  shear  elements  perpendicular  to 

the  cracked-uncracked  plies  is  taken  as  the  thickness  of  the  adhesive*  The 

shear  elements  have  rectangular  prismatic  shape.  If  there  are  any  partially 
bonded  areas  around  the  crack,  they  can  be  incorporated  in  the  model  by  chang- 
ing the  properties  of  the  shear  elements  in  the  partially  bonded  areas. 

In  this  finite  element  analysis,  the  cracked  element  developed  by  Pro- 
fessor Plan,  of  MIT  (Reference  7),  has  been  used.  This  cracked  element  has 
been  incorporated  in  NASTRAN  by  Northrop.  Figure  4 shows  the  NASTRAM  plot  of 
a two-ply,  adhesively  bonded  panel  with  a crack  in  one  ply. 

The  adhesively  bonded  structures  having  different  widths  (e*g,,  a plate 
with  an  adhesively  bonded  stringer)  can  be  modeled  in  the  manner  described 
above.  In  this  case,  only  the  bonded  region  of  the  two  layers  will  have  a 
similar  grid  poitit  and  finite  element  model.  The  NASTRAN  plot  of  a cracked 
plate  with  an  adhesively  bonded  stringer  is  shown  in  Figure  5. 
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Figure  1,  Adhes  ively  Bonded  Strut  ture  wi  th  a Debond  and  Crack 
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Finite  Element  Model  of  an  Adhesively  Bonded  Structure 
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Figure  3 + Cracked  Element  Around  Crack  Tip  and  Fine  Mesh  Around  Debond  in  Finite  Element  Model 
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Figure  4.  NASTRAN  Finite  Element  Model  of  Finite  Width,  Center  Crack 
Tension  Bonded  Panel  (one^quarter  of  the  structure  shown) 


Figure  5.  NASTRAN  Finite  Element  Model  of  a Cracked  Plate  with  a Bonded 
Stringer  Geometry  (one-quarter  of  the  structure  shown) 

This  finite  element  analysis  can  also  be  used  for  bonded  structures  with 
more  than  one  cracked  Layer*  In  such  a structure*  cracked  elements  are  pro- 
vided ahead  of  each  crack  tip,  in  each  cracked  layer*  The  stress  intensity 
factors  for  each  crack  are  computed  from  displacements  and  stiffness  matrix 
of  each  cracked  element* 

2.2  MATHEMATICAL  METHODS  OF  ANALYSIS 

Mathematical  methods  of  analysis  are  very  useful  In  analyzing  crack 
problems  in  bonded  structures.  This  is  particularly  so  where  boundaries 
are  regular  and  the  crack  lengths  are  small  compared  to  the  dimensions 
of  the  structure.  However*  certain  simplifying  assumptions  have  to  be 
made  in  the  mathematical  modeling  of  bonded  structures  (Section  2,  Vol- 
ume I ) . 

A complex  variable  approach  to  reduce  the  problem  of  cracked  struc- 
tures to  the  solution  of  integral  equations  is  considered  convenient,  and 
is  used  in  these  investigations.  The  following  paragraphs  describe  the 
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mathematical  formulation  of  three  classes  of  structural  problems:  (1)  a 
cracked  plate  with  a partially  debonded  stringer,  (2)  a cracked  plate 
stiffened  by  two  partially  debonded  stringers  and  symmetrically  located 
about  the  centerline,  and  (3)  a two- layer,  adhesively  bonded  structure 
with  a debonding  crack  in  one  Layer.  The  plate  is  assumed  to  have  a through- 
the- thickness  flaw  and  a debonding  crack  in  the  adhesive  around  the  main 
crack,  which  propagates  under  fatigue  loads.  The  debond  in  the  adhesive 
also  propagates  as  the  crack  in  the  plate  propagates, 

2,2*1  A Cracked  Plate  Stiffened  by  a Partially  Debonded  Stringer  at  an 

Arbitrary  Location 

The  cracked  plate  with  an  adhesively  bonded  stringer  shown  in  fig- 
ure 6,  can  be  represented,  due  to  superposition  techniques,  by  the  two 
prob  Lems  of  Figures  7a  and  7b*  The  structure  in  Figure  7a  has  no  crack, 
hence  there  are  no  stress  singular! ties  or  shear  stress  in  the  adhesive, 
therefore  the  problem  of  determining  stress  intensity  factors  and  shear 
stresses  in  the  structure  of  Figure  6,  reduces  to  analyzing  the  perturba- 
tion problem  of  Figure  7b,  The  solution  of  the  problem  consists  of  the 
following  three  steps: 

1 , Formulation  of  the  problem 

The  perturbation  problem  is  formulated  under  the  assumption  of 
generalized  plane  stress.  The  thickness  of  the  adhesive  is 
assumed  to  be  small  compared  to  the  thickness  of  the  plate, 
hence  the  adhesive  is  treated  as  a shear  spring  (Reference  3), 

The  shear  stresses  transmitted  through  the  adhesive  are  treated 
as  body  forces  in  the  plate  analysis.  The  stringer  transmits 
body  forces  to  the  crack  plate,  as  shown  in  the  free  body  dia- 
gram of  Figure  8,  It  is  further  assumed  that  the  stringer  area 
is  concentrated  on  one  point  at  the  centerline  of  the  stringer. 
Therefore,  the  body  forces  transmitted  will  be  concentrated 
along  one  line.  Along  the  centerline  of  the  stringer,  the  dis- 
placement continuity  between  sheet  and  stringer  may  be  written  as 
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h 

v (y)  - vs(y)  = p(y)  . y on  l (i) 

^ S^a 


where  L denotes  the  line  x = d,  the  distance  between 

the  center  of  the  crack  and  the  center  of  the  stringer  from  b 
(the  debond  size)  to  “* 


The  terms  v^(y)  and  v^Cy),  are  displacements  in  the  plate  and 
stringer,  respectively,  ^ is  the  shear  modulus  of  the  adhe- 
sive, is  the  width  of  the  stringer,  and  P(y)  is  the  contact 
shear  force  per  unit  Length  of  the  stringer  at  the  y location. 

The  displacements  of  the  plate  v and  the  stringer  v can  be 

P s 

written  as 


o 


vs(y)  - 


ks(y^o>  P(yo}  dyo 


y on  L 


(2) 


LI 
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Figure  7.  Superposition  Technique  for  a Cracked  Plate  with  an  Adhesively  Bonded  Stringer 
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Figure  8.  Free- Body  Diagrams  of  the  Plate,  Adhesive,  and  Stringer 

where  k^(y)  is  the  solution  of  the  plate  with  a crack  subjected 
to  uniform  unit  pressure,  kp(y*Y0>  is  Green's  function  for  the 
pLate  with  a stress-free  crack,  subjected  to  concentrated  forces 
as  shown  in  Figure  8f  and  k (y,y  ) is  the  displacement  at  the  y 
location  (at  x = d)  due  to  the  forces  (P  = 1)  acting  at  (Ref- 
erence 4)  locations. 

For  the  stringer  it  is  found  that 


ks(y,y0)  = 


ks(0,yo)  = 0 (3) 


-X.- 

AsEs  f 

V7  * 


y<y, 


y>yt 


where  E and  A are  the  elastic  modulus  and  the  cross-sectional 
s s 

area  of  the  stringer,  respectively. 

In  Equation  (1),  kQ(y)  is  given  by  (Reference  2) 


ko<y)  - — 
p 


(k  4 1)  1m  / z'  - a2  4 (1  - <)y 


2y  Re 


/ z - a 


k^(o)  - 0 , | x | > a 


(4) 


where  k = ( 3 - v )/(i  4 v ),  u and  v are  the  elastic  constants 
P p p p 

for  the  plate,  a is  the  half-crack  length,  d is  the  distance  from 
the  centerline  stringer  to  the  midpoint  of  the  crack  and 


z = d 4 iy 


y on  L 


( V) 


To  obtain  the  kernel  k (y,  y ),  the  solution  given  in  Reference  8 

P ° 

for  a plate  with  a stress-free  crack  and  subjected  to  concentrated 
forces  X 4 iY  acting  a is  used.  For  a cracked  plate  subjected 
to  body  forces,  the  displacements  at  z are  given  by 


2p-  (u  4iv  ) = - kST  log(z  - z ) 4 log  (z  - z 
p p p o ° 

+ U)  - C 2 ) +(  Z “ ~ 

t \ z - z„ 


(z  - z)  $^(z)  4 rigid  body  displacement  (b) 


IU 


where 


2<t>0  ( z ) = S [log(z-z0)  - I(zQ)  I -j  { 2 , zQ ) - 1 2 ( z , zQ ) 

- < 1og(z-zQ)  + k I2(z,i0)  + k I(zQ)  I1(z,z0)] 

+ <Vzo>  1 [ -Jr-  + n-z0)  1 3(2. i0)  + 

z-zo 

+ J(z0)  ^(z.Zq)]  (7) 


Vz)  = 


ci (z)  - Kz0n  - -*f-  ckz)  - i(z0>] 


2/z-a^  ' 


+ <zo‘5o)  s [ 

(z-z„) 


I(z)  - I(iJ  J(zJ 


]} 


2 -Z, 


(8) 


and 


ii<2'2°)  ■ -77: 


log 


z -a 
o 


2 , /2  2 /T  2 

z z - a 4-  /z  - a /z  - a 
o y o v 

a(z-z  ) 
o 


1 9(z,  z ) = log  f—  + \/(z/a)2  - l"1  + z I (2,  z ) 

t O L 3 — 1 O L O 


. r / 2 s 

l 2, 2 ) = — ; — -rr~  — — + z I ! { 2f  2 } 

3 o 2 2 z-z  o 1 o 


a ~ z L o 
o 


! (z,z  ) = I1  (z^z  ) + z I (z,z  ) 

4 O L D O J O 


(9) 
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1(2)  = /z2-a2  - z 


J(z)  = —=  - 1 

2 = x + iy 

20  = xo  + iyo 

<.  - X + 1Y 
2irhp(  1 +k) 

Combining  Equations  (6)  through  (9),  the  displacements 

and  v can  be  written  as 
P 

2yp ( Up+i Vp ) = S {-  ic[log(z-z0)  + log(z-zo)] 

+ [®](z»Zq)  + 9-j  ( z , zQ ) ] - 2*  [ ® i ( z » Zq ) + < 

+ (SjijCie  e2(z)  - e2(i)]  + ( |2l!»  ) e5u,z0 

+ S [(  ) - 1 + k 83(2. z0)  - 03(z. 2„)  - 

z'zo 

+ k e4(z,z0)]  + rigid  body  displ. 


(10) 

8,(2, i„n 

n 

84(2, 20) 

(11) 
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where 


01(z,zo)  = log[zz0  - a2  + /zQ2-a2  /z2-a2  ] 


02(z)  = log[z  + /z2-a2  ] 


i z -z  /_2  ,2  A 2 .2 

9,(1. Zj  - i { -V2  )[1  ♦ /Z  '%  /Z?  ~a  ] 


z-z. 


1 / z-z 


2-2 
a “zo 


94(z»z0)  = 2 < 77" 

■o 


I(z)  - I(zJ 


z-z-  /P-I7 


^ ] 


e5(z.z0)  ■ e4(z-5o>  ■ 7 ItU  J(20> 


/FTa2 


C 12) 


Putting  Y = 1 at  z = d + iy  and  Y = 4 at  z = d - iy  (X^0)f 

and  taking  the  imaginary  part  of  the  final  expression  to  obtain 

k (y,  y ) which  relates  a unit  load  (Y)  located  at  a to  the  y 
p ' o o J 

direction  displacement  (v) 


kp(y,y0) 


1 

4TTliphp  ( 1 +K) 


<[log ( z - z 0 ) 


+ 1 og ( z- 2q  ) ] 
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+ 03^»zo^  * 93^,zo^  + ^ ^ ° ) 05^z,zo^ 

z"zo 

- C ^ ) e5(2,20)f 

z"zo  J 


(2  .ZQ)  ] 
9 3 ( z ’ 2 o ) 1 


(13) 
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Thus,  from  Equations  (1)  and  (2),  we  obtain  the  integral,  equation 


P(sO  + / k(y,y0)  P(y0)  dy0  - - 

Lp 


P k0(y)  » y on  l 


(U) 


a 


where 


(15) 


The  kernel  k(y,yo)  has  only  logarithmic  singularities  and  is 
integrable  on  L,  Equation  (14)  can  be  treated  as  a Fredholm 
equation  of  the  second  kind, 

2 , Stress  Intensity  Factors 

The  stress  intensity  factor  for  the  normal  stress  cr  is  defined 

y 

as 

k = lim[/ 2(x-a)l  a (x,0) 
a " y 

and  obtained  as  follows  (Reference  2) 


(17) 


where 


(18) 
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and 


a 

o 


a for  right  tip 

-a  for  left  tip 


z - d + iy 
o J o 


(19) 


3 . Numerical  Solution  of  Integral  Equations 

The  Fredholm  equation  given  in  Equation  (14)  can  be  solved  numeri- 
cally by  reducing  it  to  a system  of  algebraic  equations.  A simpLe 
way  of  doing  this  is  to  replace  the  integral  by  summation,  and  use 
collocation  to  obtain  an  n x n system  of  equations.  Hence, 


P(y^ ) + l k(y1  ,y j ) P(y^ ) Ay  ^ 

0 1 


dspa 

9 ko^yi  ^ » 1=1 N 

d 


(20) 


which  approximately  gives  P(y ) , j = 1,2,  ,,,  ,N#  Similarly,  for 
the  stress  intensity  factor,  Equation  (17)  can  be  written  as 


+ “ l o(yJ  P(yJ  Ay, 
ao  j=l  J J 


(21) 


The  kernels  in  the  integral  equation,  (14),  have  logarithmic 
singularities,  hence  the  singular  part  of  the  kernel  k(y^,  y^) 
in  Equation  (20)  is  evaluated  separately  as  an  integraL  along  the 
segment  Ay,,  The  number  of  collocation  points  N,  are  chosen  in 
such  a way  than  an  increase  in  N does  not  affect  the  results. 

The  computer  program  for  the  solution  of  this  problem  is  described 
in  Appendix  A. 
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2.2,2  A Cracked  Plate  Stiffened  by  Two  Partially  Debonded  Stringers, 

Symmetrically  Located  About  the  Centerline 


A cracked  plate  with  two  symmetrically  located,  partially  debonded 
stringers,  is  shown  in  Figure  9a,  The  actual  problem  of  determining 
stress  intensity  factors  and  shear  stresses  in  the  adhesive  for  a plate 
uniformly  loaded  far  from  the  crack,  reduces  to  the  perturbation  problem 
shown  in  Figure  9b,  similar  to  the  single-stiffener  problem  discussed  in 
Paragraph  2,2.1,  The  free  body  diagram  with  the  body  forces  acting  on  the 
pLate  is  shown  in  Figure  10,  The  body  forces  transferred  to  the  plate 
are  symmetrically  located  as  shown.  Due  to  symmetry,  the  problem  may  be 
formulated  for  the  unknown  shear  stresses  in  one  stringer  only. 

The  displacement  continuity  between  sheet  and  stringer  is  given  by 
Equation  (l).  The  displacement  in  the  stringer  v$(y)  is  given  by  the 
second  of  Equations  (2)  and  the  displacement  in  the  plate  v^(y),  by 

y?>  - pyy) + Jl,  y^y  p(y  % <22 

where  denotes  the  lineX  = d,  b ^ y < ® and  x = -d , b £ y < «,  Using  the 
displacements  given  by  Equation  (11),  the  integral  equation  for  this  problem 
can  be  written  as 


P(y>  + 


d m> 

s a 


/ 


y^y 


- / Kp(y,yo)P(yo)dyo 


d |o. 
s a 


k0(y) 


(23) 


where  L denotes  the  line  x = d,  b 5 y < ® 

The  kernels  k^(y,yo)  have  logarithmic  singularities  and  are  inLegrable  on 
. Equation  (23)  can  be  treated  as  the  Fredholm  equation  of  the  second 
kind.  The  stress  intensity  factors  in  this  case  at  both  crack  tips  are 
equal  and  are  given  by 

— = p + i/L  a(y0)p(yo)dyo  (2i° 

/ a 
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Figure  9P  Cracked  Plate  with  Two  Adhesively  Bonded  Stringers 
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Figure  10.  Free-Body  Diagrams  of  the  Plate,  Adhesive,  and  Stringers 


where  Of(yQ)  is  given  by 


2ttH  ( 1+  k ) a 

P 

-(140  L a ■ + 

a - Zq 

+ 0 + «>  r*-1-  wi  - a + *>  1 

a H-  Zq 

+ (zo  - Zq  ) 

a + 1(-2q  ) . J(^j  ) 

/ y3  a - Zb 

+ (zq  ~ zq  ) 

a + I ( - zb  ) , J(  - zq  ) 

/ * a + Zq 

(a  + Zb  ) 

+ ( Zq  - Zq  ) 

a + I(- Zo  ) ( J('Zo  ) 

a - zb 


a + 2b 
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+ (zo  - 26  ) 
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(25) 


where  zq  = d + iy0 

The  integral  equations  are  solved  using  collocation  as  described  in 
Paragraph  2.2.1,  and  stress  intensity  factors  are  computed. 

The  computer  program  for  this  problem  is  detailed  in  Appendix  B. 

2.2.3  Two-Layer  Adhesively  Bonded  Structure  with  a Debonding  Crack  in 

One  Layer 

The  analyses  discussed  in  the  previous  two  paragraphs  were  concerned 
with  a cracked  sheet  bonded  to  a sound  layer  of  small  width  compared  to 
the  sheet.  In  this  problem,  the  cracked  sheet  and  sound  Layers  are  of  the 
same  width,  and  infinitely  wide. 

Consider  the  adhesively  bonded  structure  of  Figure  11a,  consisting  of 
two  plates  with  thicknesses  hj  and  h2 , respectively,  and  bonded  through  an 
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adhesive  layer  of  constant  thickness  ha*  The  plate  is  subjected  to  forces 
Tx  and  Ty  Per  unit  length  of  the  plate,  Plate  1 is  assumed  to  have  a through- 
the- thickness  flaw  and  a debonding  crack  in  the  adhesive  around  the  main  crack. 

In  Figure  11a,  the  size  of  the  debond  is  shown  to  be  the  same  as  the 
length  of  the  crack.  Experiments  confirm  that  this  will  generally  be  the 
case  when  the  initial  flaw  has  grown  due  to  fatigue,  and  the  edge  of  the 
debond  will  propagate  with  the  crack  tip.  If  the  initial  debond  size  in  the 
adhesive  is  large,  it  is  possible  that  the  crack  length  will  be  shorter  than 
the  debond,  however,  this  problem  can  still  be  formulated  in  the  same  manner 
as  the  problem  of  Figure  11a. 

The  analysis  of  the  adhesively  bonded  structure  in  Figure  11a  will  be 
based  on  the  following  assumptions: 

1,  The  thickness  of  the  plates  (hi  and  h? ) is  small  compared  to  the 
inplane  dimensions  so  the  structure  may  be  considered  to  be  under 
a generalized  plane-stress  loading, 

2,  Plates  1 and  2,  and  the  adhesive  layer,  are  homogeneous  and  linearly 
elastic, 

3,  The  thickness  variation  of  stresses  in  Layers  1 and  2 is  neglected* 

4,  The  thickness  (ha)  of  the  adhesive  is  small  compared  to  the  thick- 
ness of  the  plates,  therefore  the  adhesive  may  be  treated  as  a shear 
spring  rather  than  as  an  elastic  continuum, 

5,  The  surface  shear  transmitted  through  the  adhesive  acts  as  a body 
force. 

Let  ux  , v1  and  u2  , v2  be  the  x,  y components  of  the  inplane  displace- 
ment vectors  in  Materials  1 and  2,  respectively,  and  T ^ and  be  the  compo- 
nents of  the  shear  stress  acting  on  the  adhesive.  Assumption  4 above,  gives 
the  following  continuity  conditions : 


Ui 


vi 


(26) 
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(27) 


(c)  Perturbation  problem 


Figure  LU  Superposition  Technique  for  an  Adhesively  Bonded  Plate 
with  a Debond in^  Crack 
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where  |^a  is  the  shear  modulus  of  the  adhesive.  The  body  forces  acting  on 
Plates  L and  2 are  given  by 


T t 


where  Xx  $Yl  are  the  body  forces  in  Plate  lf  and  X2,Y?,  the  body  forces  in 
Plate  2 at  the  point  having  coordinates  (x,y). 

The  solution  to  the  problem  in  Figure  11a  can  be  obtained  by  superimposing 
the  two  cases  shown  in  Figure  lib  and  11c,  In  Figure  lib,  the  disp lacements 
for  the  two  plates  are  such  that  ui  = u?  , and  vl  = v2  throughout  the  plate, 
hence  there  are  no  interface  shear  stresses  or  body  forces  in  the  plate. 

Thus,  the  stress  intensity  factors,  interface  stresses,  and  boundary  of  the 
debonding  crack  of  Figure  11a  will  be  given  by  the  perturbation  problem  of 
Figure  lie.  The  solution  to  the  problem  of  Figure  11a  is  carried  out  in  the 
following  steps, 

1 . Determination  of  Applied  Crack  Surface  Stress 

The  crack  surface  stress  p^  to  be  applied  to  the  crack  surface  in 
the  perturbation  problem  of  Figure  lie,  is  equal  to  the  stress 
taken  by  the  plate  of  Material  1 in  Figure  11b,  Let  T ^ and  T^2 
be  the  forces  per  unit  length  In  x direction,  taken  by  plates  of 
Materials  1 and  2,  respectively.  Similarly,  and  T^0  are  the 

forces  per  unit  length  in  y direction,  taken  by  plates  of  Mate- 
rials 1 and  2,  respectively.  The  equilibrium  conditions  in  x and 
y directions  give 


T , + T = T 
xi  x2  x 


T , + T - = T 
yl  y2  y 


(29) 

(30) 
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The  continuity  of  displacement  in  the  problem  of  Figure  lib,  gives 


U\  - Us 


or 


(Txl  - vl  Tyl,J 


hi  Ex 


Tx2  ' V2  Ty2 

h2E2 


(31) 


Similarly 


(Tvl  ~ h Txl)  (T»i  • "i  T«z) 

hlEl  h2E2  (32) 

Equations  (29)  through  (32)  can  be  solved  for  TXx , Tyg , Tyx  , 
and  Ty^ . The  stress  p0  is  equal  to  / h x and  is  given  by 


P0  = 


fa.T  + a„T  ) 
1x2  y/ 


(33) 


where 


V e2  “ih  \ 

al  - h^(h1E1  * h2E2  - „1h2E2  7 jpqrJ 

Vs  E1 

a2  hiEi  + h2E2 

V2E1  * VlE2  hlEl  + H2E2 

a3  - hjEj  * h2E2  - - VlL2 


(34) 
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2. 


Formulation  of  Per tur bat ion  Problem 


The  perturbation  problem  will  be  looked  upon  as  consisting  of 
Plate  2 with  body  forces  in  Region  D of  Figure  12a,  and  Plate  1 
with  body  forces  and  loading  on  the  crack  surface,  as  shown  in 
Figure  12b, 

The  stresses  and  displacements  in  terms  of  complex  functions  cc, 
fi,  are  given  by  (Reference  9) 


<7+0 

x y 


2\0(Z)  +0(2)] 


(35a) 


a - cr  + 2ir  = 2[  (z  - z)0'(z)  - Q(z)  + R{z)] 
j x xy 


(35b) 


2pt<u 


f f 

+ iv)  = k / ®{z)dz  - / 

■'n  Jn 


£2<z)dz  - (z  - z)O(z) 


(35c) 


a.  Displacements  in  Uncracked  Plate 

For  an  infinite  plate  under  the  body  forces  X2 ,Yg , acting  at 
at  z3  (>fe,yo)f  the  complex  functions  are  given  by  (Reference  8) 


0{z)  = 


S 


z 


z 


0 


£2(9)  = 


S = 


+ iY2 
2 71(1  + K) 


(36) 
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(a)  Plate  of  Material  2 with  body  forces 


x 


(b)  Plate  of  Material  1 with  body  forces  and  crack 

surfaces  loaded 


Figure  12*  Loadings  on  Plates  of  Materials  l and  2 in 
Perturbation  Problem 
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Substituting  Equation  36  in  Equation  35c,  the  displacements 
in  the  uncracked  plate  with  body  forces  (plate  of  Material  2) 
can  be  obtained.  These  displacements  are  given  by 


2m^(u?+  tvs  ) = - <S[  log (z  - zq  ) + log(z  - To )] 
+ < Sf !og(z  + Z0 ) + log { 2 + Zq  )1 


1 

N 

<5" 

_ 1 

c 

Z + z0 

Z - z 0 

* I 

- 

— 5 

1 

-t 

1 

+ 

IN 

- <S[log(z  - Zo)  + log(z  - Zo)] 
+ <S[log(z  + z0 ) + log ( z + zc )] 


Z - Zo 


Z + Z0 


+ s 


-1 


Z - z0 


- s 


- 1 


z + Zq 


(37) 


Assuming  that  the  body  forces  XP t are  continuous  functions 
of  Zq  or  (kq,  y0 ) as  defined  in  region  Df  the  displacements  at 
point  2 = x + iy  in  the  plane  are  given  by 

«e(*»y)  = JT  [Kn  (xty;xo  ,y0  )X2(xc  (y0) 

D 


+ ^12  (x,y  ;xo  ,y0  )Yg  (xq  .yo  )]dx<j  dy0 

vs(x,y)  = JJ  [ Ka  i ( x , y ; xq  yo  ) X2  ( Xq  , y0  ) 

D 

+ K2a  (x,y;xcy0  )YP  (xb  ,y0  )]dxbdy0 
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where  the  kernels  K (i,j  = 1,2)  are  given  by  Green's  functions 
(Equation  (37)), 

b,  Displacements  in  a Cracked  Plate 

The  displacements  in  the  cracked  plate  will  consist  of  two 
components,  one  due  to  uniform  applied  pressure  on  the  crack 
surface  (Reference  9)  and  the  other  due  to  body  forces  acting 
on  the  plate  surface.  The  displacements  due  to  uniform  applied 
pressure  are  given  by 


uu(*.y>  - - 1,Re 


/ 2 2^ 
{ z - a ) 


1/2 


- 2ylm 


2 2. 

z - a ) 


172 


+ (1  - K )X 


= P0fl.<X*  y> 


V 


11 


(X’y)  = 4^ 


f 

2 2 1/<2 

- 2yRe 

z 

<«  + l)Im 

z - a 

2 2 1/2 

l 

- a ) J 

+ (i  - K)y 


= P0f2(x,y) 


(39) 


where  < - (3  - va) / ( 1 + v2) , and  vx  are  elastic  constants* 

To  find  the  solution  of  the  cracked  plate  under  distributed  body 
forces,  the  solution  for  the  concentrated  forces  Xx ,YL  acting 
at  z0  is  used  as  Green’s  function*  This  solution  is  given  by 


*(z) 


V2) 


Q(a)  = 


{z0  " Vs 
f~rj * V 

(z  - z0) 


(40) 
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Viili 

2tt(1  + «) 


Vz>  ■ r/*2-*2) 


-1/2 


z - z. 


[l(Z)  - I(z0) 


(41) 


Substituting  Equation  (40)  into  (35c),  the  displacements 
are  obtained * These  displacements  in  the  cracked  plate  due 
to  body  forces  acting  at  an  arbitrary  location,  are  given  by 
Equations  (LL)  and  (12)*  Assuming  again,  that  the  body  forces 
X1  , Y1  are  continuous  functions  of  (x^  , y0 ) ) defined  in  region  D, 
using  Equations  (11)  and  (12)  as  the  Green’s  function,  the  dis- 
placements in  the  cracked  plate  due  to  body  forces  Xl  and  Yi  , 
may  be  expressed  as 

Ul2(x,y)  = H [Hu(x,y;x0,y0)Xl(x0Jy0) 


+ 


H12(x,y; 


0’ 


y0)Yl(Vy0)]dViy0 


(42a) 
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v12(x,y>  = / /[H2l(x,y;x0)y0)X1{x0,y0) 


H22(x,y;  VWVVK^O 


The  total  displacements  in  the  base  plate  of  Material  1 may  then 
then  be  obtained  by  adding  Equations  (39)  and  (42). 

ui(x,y)  = uxl  + u13,  vj(x,y)  = vn  + v13  (43) 


c.  Integral  Equations  for  Tx,  Ty 

The  x and  are  the  x and  y components  of  the  shear  stress 
on  the  adhesive  layer.  The  relationship  between  body  forces 
and  shear  stresses  is  given  by 


X 


1 


(44) 


X 


2 


Use  of  the  displacements  ui  , u2 , V!  , and  ve  in  Equations  (26) 
and  ( 27 ) yields  the  following  system  of  integral  equations 

for  the  unknown  functions  t and  t 

x y. 
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h 

it  y**» + w [ kn<x-y>  vww 


+ kl2<x’y;xOyO)Ty,xO'yO)]dxOdyO  = pOfl(X’y) 


h 

57Vx,y)  + y 

+ k22(x>y,xOyO)Ty<X0'yO)]dxOdyO  " P0£2(X'y) 


(x,y)  € D (45) 

The  functions  f; , fa  are  given  by  Equation  (39)  and 


H11 

k = — — + 
Kll  ht 


!u  , Hl2  ^12 

h2  * K12  ' hj  h2 


H 


*21 


21 


K 


21 


H 


22 


K 


22 


22 


(46) 


The  kernels  K,^,  (i,j  = 1,2),  which  have  logarithmic  singu- 
larities, are  known  and  are  square  integrable  in  region  D. 
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3* 


The  Stress  Intensity  Factors 


In  general,  the  stress  intensity  factors  in  plane  stress  problems 
may  be  expressed  in  terms  of  the  Kolosov-Muskhel lishvili  function 
$ (z)  as  follows  {Reference  8) 


ki  ’ 


ik„ 


= xl™a  1 2(x  - a)] 


1/2  | 


ffy(x»  0)  - irxy(x,0) 


1/2 

= zlima  2[2(z  - a)]  *{z)  (47) 

In  the  problem  considered  here,  the  shear  component  of  stress 
intensity  factor  k2  , is  zero  due  to  the  symmetry  in  Loading  and 
geometry.  The  cleavage  component  ki  , Is  determined  by  adding  the 
effects  of  the  crack  surface  pressure  po  and  the  body  forces  X^  and 
Yi  . The  stress  intensity  factor  k*  can  be  written  as 

kl  = P</a  * y [ hl(*0,y0,xi<*0,y0) 


+ h2(VVYl<X0’y0)|‘Vy0 


(48) 


where  (xq  ,y0  ) and  h2  { xq  , y0  ) are  the  cleavage  components  of  the 
stress  Intensity  factor  due  to  the  concentrated  body  forces  Xl  and 
Y| , respectively  (Figure  12)  and  are  obtained  by  substituting  Equa- 
tion (40)  into  Equation  (47),  Substitution  of  Equation  (40)  into 
Equation  (47)  gives 


ki  - ike  = 


77 


s — - I^z°  n + <s  ra  + liin  )n 

a - z 0 - 

a - zq 


- S(zq  -To)  a+l(J*)  i 

(a  - z0  Y a - Zq  . 


(49) 
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hi{xo,yo)  and  hg  <*o  *Yq)  are  coefficients  of  X1  and  Y1  respectively, 
in  real  part  of  Equation  (49)* 

4*  Solution  of  Integral  Equations 

The  system  of  integral  equations  given  by  Equations  (45),  are  the 

Fredholm  type  and  may  be  solved  by  using  the  standard  numerical 

techniques.  In  this  case,  it  is  done  by  dividing  the  region  D 

into  smaller  cells,  and  unknown  functions  I and  T are  assumed  to 

x y 

be  constant  in  each  cell.  Thus,  using  a numerical  integration 
scheme  to  evaluate  the  integrals,  the  integral  equations  are  reduced 
to  a system  of  algebraic  equations*  The  kernels  in  the  integral 
equations  have  logarithmic  singularities,  hence  the  singular  part  of 
the  kernels  is  evaluated  separately,  in  closed  form.  In  the  actual 
integration,  a telescopic  grid  is  used.  A typical  grid  layout  used 
in  integration  is  shown  in  Figure  13,  where  only  one-quarter  of  the 
plate  is  shown  because  of  symmetry.  The  cell  size  is  kept  small 
above  the  crack  plane  for  a distance  of  about  a half-crack  length, 
as  the  shear  stresses  are  high  in  this  region  and  are  maximum  at  the 
boundary  of  debond.  The  cell  size  is  progressively  increased,  as 
shown.  The  debonded  region  is  approximately  represented  by  straight 
lines  for  integration  purposes.  The  boundary  of  the  domain  of  inte- 
gration goes  to  infinity,  hence  the  size  of  Region  D is  restricted 
in  numerical  analysis  such  that  the  stress  intensity  factors  are  not 
appreciably  affected , 

A generalized  program  to  obtain  stress  intensity  factors  could  not 
be  developed,  as  convergence  of  the  solution  is  dependent  on  the 
thickness  and  material  properties  of  adherends  and  adhesive.  The 
convergence  of  the  solution  also  depends  on  crack  length.  In  the 
integral  equations  (45),  if  h^a  Is  very  small  (less  than  0,1  x 1QL  ) , 
the  solution  will  not  converge,  and  shear  stresses  will  oscillate. 

This  can  be  avoided  by  decreasing  cell  size  further  and  increasing 
integration  points.  Two  computer  programs,  one  for  small  crack 
length  (a  £ 0.4)  and  the  other  for  large  crack  lengths  (0*5  ^ a £ 1.0) 
are  given  in  Appendices  C and  D,  respectively. 


36 
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Figure  13,  Subdivision  of  Bonding  Area  D for  Numerical  Analysis 
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2 ■ 3 COMPARISON  OF  FINITE  ELEMENT  AND  HATEHMATICAL  ANALYSES 


In  Paragraphs  2.1  and  2.2,  the  finite  element  and  complex  variable 
approaches  to  ana Lyzing  crack  prob Lems  were  discussed , In  each  of  the  analysis 
methods,  certain  simplified  assumptions  were  made  and  the  results* obtained 
by  the  two  methods  for  the  same  problem, were  compared.  In  the  following 
paragraphs, the  problems  of  a cracked  plate  with  a bonded  stringer,  and  two 
adhesively  bonded  plates  with  a crack  in  one  plate,  are  considered  for  com- 
parison purposes. 

2,3*1  A Cracked  Plate  with  a Single,  Centrally  Located,  Adhesively  Bonded 

S L ringer 

Figure  14  describes  the  structure  (with  dimensions)  used  for  analysis 
comparison.  Complex  variable  analysis  assumes  that  plate  is  infinitely  large. 
The  finite  element  model  used  in  the  analysis  is  shown  in  Figure  5.  The  plot 
of  normalized  stress  intensity  factors  K^/c/rfa  versus  the  half-crack  length  a, 
obtained  by  two  methods  of  analysis,  is  shown  in  Figure  15.  It  is  seen  that 
the  results  of  the  finite  element  analysis  and  complex  variable  approach  are 
close  only  for  small  crack  lengths  (a  ^ 0.1  inch).  For  half-crack  lengths 
larger  than  0.1  inch,  the  stress  intensity  results  obtained  by  the  complex 
variable  technique  are  much  larger  than  those  obtained  using  the  finite 
element  analysis.  In  the  complex  variable  formulation,  the  load  transfer 
for  the  stringer  is  concentrated  along  the  centerline  of  the  stringer,  hence 
the  width  of  the  stringer  in  Load  transfer,  for  short  cracks  (under  the 
stringer),  from  the  cracked  plate  to  the  stringer  is  not  fully  effective. 

In  the  finite  element  analysis,  the  entire  width  of  the  strap  is  modeled, 
allowing  for  load  transfer  over  the  entire  width  to  occur,  and  thereby  yields 
lower  stress  intensity  factors  compared  to  the  complex  variable  approach  for 
short  cracks  (under  the  stringer). 

To  obtain  better  solutions  using  the  complex  variable  approach,  it  was 
assumed  that  the  area  of  the  strap  is  concentrated  at  two  points  that  are 
symme trica 1 ly  located  about  the  centerline  of  the  plate,  as  shown  in  Figure  9, 
These  areas  were  assumed  to  be  concentrated  at  the  crack  tip  if  2a  ^ B 
{B  = width  of  the  strap),  and  at  the  edge  of  the  strap  if  2a  > B,  The  prob- 
lem was  formulated  using  a complex  variable  approach.  Due  to  symmetry,  the 
problem  was  reduced  to  the  solution  of  an  integral  equation  for  unknown 
shear  stresses  in  one  stringer  only. 
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Figure  14.  Geometry  of  a Cracked  Plate  with  an  Adhesively 
Bonded  Stringer  (no  debonding) 

SmalL  Crack  (2a  < dg) 


39 


Normalized  Stress  Intensity  Factors  K._/  a 


1.0 


Figure  J5.  Comparison  of  Finite  Element  Determined  Stress  Intensities  with 
Integral  Equation  Solutions  (Cracked  Plate  with  Bonded  Stringer) 
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The  anaLysis  of  the  probLem  was  carried  out  using  the  complex  variable 
approach  discussed  in  Paragraph  2*2,2.  The  plot  of  K ^ / a / rr  a versus  the  half- 
crack length  a for  the  two  symmetrically  placed  stringers  is  also  plotted  in 
Figure  15.  It  is  seen  that  up  to  a half-crack  length  of  about  0.3  inch,  the 
results  given  by  the  finite  element  analysis  and  the  integral  equation  approach 
are  very  close.  Beyond  this  point,  the  integral  equation  approach  gives  higher 
values  of  the  stress  intensity  factor*  The  difference  between  Lhe  results  of 
the  finite  element  analysis  and  the  integral  equation  approach  increases  as 
the  crack  length  increases,  up  to  a half-crack  length  of  about  one  inch. 

Beyond  a one- inch  crack  length,  the  difference  decreases.  At  a half-crack 
length  of  one  inch  (total  length  equal  to  the  strap  width),  the  results  of 
the  integral  equation  approach  are  about  33  percent  higher  than  those  of  the 
finite  element  analysis.  At  a half-crack  length  of  two  inches,  the  results 
of  these  two  methods  of  analysis  differ  by  only  two  percent.  The  resuLts  of 
Figure  15  indicate  that  the  complex  variable  approach,  assuming  a concentrated 
stringer  area  at  one  point,  may  be  used  for  very  small  or  very  large  crack 
lengths  only*  For  larger  crack  lengths,  either  the  finite  element  method  or 
the  two-stringer  formulations  may  be  used.  The  computer  run  time  for  the 
finite  element  analysis  is  about  120  seconds  (CPU  time),  whereas  the  integral 
equation  approach  computer  run  time  is  approximately  8 seconds  (CPU  time) * 

Thus,  it  can  be  seen  that  a considerable  saving  in  computer  run  time  is 
obtained  by  using  a complex  variable  approach.  It  should  be  noted  however, 
that  more  accurate  solutions  can  be  obtained  using  the  complex  variable  method 
if  the  stringer  area  is  subdivided  in  two  or  more  concentrated  strips, 

2,3*2  Two  Adhesively  Bonded  Plates  with  a Crack  in  One  Plate 

The  dimensions  of  the  structure  are  shown  in  Figure  16  for  the 
problem  of  a cracked  plate  adhesively  bonded  to  an  uncracked  plate* 

The  adhesive  is  assumed  to  have  an  elliptically  shaped  debond  with  a 
ratio  of  minor  axis  to  major  axis  equal  to  0.1*  The  leading  edge  of 
the  debond  is  assumed  to  coincide  with  the  leading  edge  of  the  crack. 

This  shape  of  debond  was  predicted  from  analysis  as  discussed  in  Section  6 of 
this  volume  and  has  also  been  observed  in  experiments  (Section  4,  Volume  D* 
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Figure  16.  Geometry  of  Cracked  Plate  Adhesively  Bonded  to  an 
Uncracked  Plate 
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The  finite  element  model  used  in  the  analysis  of  the  cracked  structure  is 
shown  in  Figure  4.  In  the  finite  element  analysis,  the  plate  width  is 
assumed  to  be  12  inches,  and  in  the  complex  variable  approach,  the  plate 
width  is  assumed  to  be  infinite.  Up  to  a half-crack  length  of  one  inch 

cl 

— = 0.0833,  the  finite  width  effect  is  negligable.  It  can  be  noted  from 
Figure  17  that  the  normalized  stress  intensity  factors  obtained  by  the  two 
methods  of  analysis  agree  very  well. 

The  computer  run  times  for  finite  element  analysis  are  about  240  sec- 
onds (CPU  time),  and  for  mathematical  methods,  the  computer  run  times  are 
48  seconds  (CPU  time),  for  small  crack  length  programs,  and  between  60  and 
90  seconds  (CPU  time)  for  long  crack  length  programs,  depending  on  the  de- 
bond size*  Thus,  with  the  mathematical  solutions,  a considerable  saving  in 
computer  run  times  is  achieved  when  a study  of  various  material /geometrical 
sensitivities  are  required* 


Normalized  Stress  Intensity  Factor  K^/a/n  a 


Finite  Element  h = 0,008fT 


Figure  l?. 


Comparison  of  Finite  Element  Determined  Stress  Intensities 
with  Integral  Equation  Solutions  (two  plates  bonded) 
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SECTION  3 


APPLICATION  OF  ANALYSIS  TO  MULTILAYERED  AND  OTHER  BONDED  STRUCTURES 


Finite  element  and  mathematical  methods  of  analysis  were  discussed  in 
Section  2,  primarily  in  relation  to  two-ply,  adhesively  bonded  structures. 
However,  these  methods  are  also  applicable  to  cracked,  multilayered,  bonded 
structures,  as  well  as  to  rivet  bonded,  we  Id bonded,  and  brazed  structures. 
The  application  of  the  methods  to  these  structures  is  discussed  in  the  fol- 
lowing paragraphs* 


3* 1 APPLICATION  OF  ANALYSES  TO  CRACKED  STRUCTURES  CONSISTING  OF  MORE  THAN 

TWO  LAYERS 

The  finite  element  and  mathematical  analysis  methods  discussed  in  Sec- 
tion 2 can  be  extended  to  multilayered,  adhesively  bonded  structures  with 
cracks  in  more  than  one  layer.  The  application  of  two  analysis  methods  to 
cracked  multilayered  structures  is  discussed  in  the  following  paragraphs. 


3*1,1  Application  of  Finite  Element  Analysis  to  Cracked  Multilayered 

Structures 

The  modified  two-dimensional  finite  element  method  described  in  Para- 
graph 2.1  can  easily  be  applied  to  cases  of  more  than  two  bonded  plates. 
Consider  the  structure  of  Figure  18,  which  has  through- the- thickness  cracks 
in  Plates  1 and  2,  and  de bonding  in  adhesive  Layers  1 and  2.  The  plate  is 
subjected  to  Force  T , in  y direction,  per  unit  length  of  the  plate*  In  the 
finite  element  modeling  of  this  structure,  a cracked  element  is  provided  ahead 
of  each  crack  tip  in  Plates  1 and  2,  If  only  one  of  the  plates  has  a crack, 
then  cracked  elements  are  only  provided  for  the  crack  tips  in  the  cracked  plate 
Plates  1 and  2,  as  well  as  Plates  2 and  3,  are  connected  by  shear  elements. 

The  shear  elements  are  assumed  to  run  continuously  throughout  the  length  of 
the  plates  in  the  bonded  regions*  In  these  regions,  the  finite  elements  are 
closely  spaced  around  the  periphery  of  the  debond  as  discussed  in  Paragraph  2,1 
and  shown  in  Figure  19,  The  case  of  three  adhesively  bonded  plates  is  shown 
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Figure  18.  Adhesively  Banded  Plates  Consisting  of  Three  Layers 


PLATE  2 
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in  Figure  18,  however,  if  there  are  more  than  three  plates  in  the  structure, 
this  method  of  modeling  can  be  extended  in  a manner  similar  to  that  described  above. 


3 , 1,2  Application  of  Mathematical  Method  of  Analysis  to  Cracked  Multilayered 

Structures 

The  mathematical  technique  described  in  Paragraph  2.2  can  be  extended  to 
the  case  of  a bonded  structure  consisting  of  more  than  two  Layers,  like  the 
structure  shown  in  Figure  18.  Each  layer  is  assumed  to  be  infinite.  The 
solution  to  the  loading  condition  shown  in  this  figure  can  be  obtained  by 
solving  the  perturbation  problem  with  crack  surface  loading  as  discussed  in 
Paragraph  2,2  The  body  forces  and  crack  surface  loading  on  each  plate  in  the 
perturbation  problem  is  shown  in  Figure  20,  Let  Txi  and  Tyx  be  the  shear  stresses 
in  adhesive  layer  1,  and  Ty^  and  Tys  be  the  shear  stresses  in  adhesive  layer  2, 

If  ui  , v1  , and  it  t Vq  are  the  x-y  components  of  the  in-plane  displacement 
vectors  in  Materials  1 and  2 respectively,  then  the  continuity  conditions  for 
Adhesive  Layer  1 are  given  by: 


Ui  -u? 


h4 


T 

Xl 


Vl 


- v2 


T 

yi 


Similar  conditions  for  adhesive  layer  2 are  given  by: 


^ - u3 


Sa 


T 

X2 


- v3 


T 

Y2 


(50) 


(51) 


where  u3  , v3  , are  x-y  are  components  of  displacements  In  Plate  3.  As  shown 
in  Figure  20,  the  formulation  for  displacements  in  Plates  1 and  3 will  be 
simiLar  to  the  case  of  two  bonded  plates  discussed  in  Paragraph  2,2,  The 
displacements  u?  and  v2  in  Plate  2 will  consist  of  3 parts:  1}  due  to  crack 
surface  loading;  2)  due  to  body  forces  Xs  and  ; and  3)  due  to  body  forces 
X3  and  Y3  , In  Green's  function  formulation  of  displacements,  the  body  forces 
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US 


Figure  20.  Plates  with  Body  Forces  and  Crack  Surfaces  Loaded 


X2  and  Y?  are  integrable  in  domain  D-Li  and  body  forces  X3  and  Y3  are  inte- 
grable  in  domain  D-L^ . The  relation  between  body  forces  and  shear  stress  is 
given  by 


Using  the  displacement  continuity  conditions  (Equations  (26)  and  (27)),  a 
system  of  integral  equations  for  t iyX  and  , iy2  can  be  obtained*  The 
integral  equations  are  of  the  Fredholm  type,  with  kernels  having  Logarithmic 
singularities*  These  equations  can  be  solved  numerically,  and  shear  stresses 
and  stress  intensity  factors  computed  as  discussed  in  Paragraph  2.2. 

3*2  APPLICATION  OF  ANALYSIS  TO  RIVET  BONDED,  WELDBONDED,  AND  BRAZED  STRUCTURES 

The  modified  two-dimensional  finite  element  analysis  discussed  in  Para- 
graph 2.1  is  applicable  to  rivet  bonded,  weldbonded,  and  brazed  structures. 

In  rivet  bonded,  weldbonded,  or  brazed  structures,  the  method  of  stringer 
attachment  is  modeled  as  a shear  element,  continuous  throughout  the  Length  of 
the  attached  portion.  In  this  type  of  modeling,  the  attached  stringer  may  be 
of  any  shape* 

Consider  the  structure  of  Figure  21,  where  a strap  is  spotweLded  and 
bonded  to  the  sheet,  using  the  weldbonding  process.  There  is  a debond  in  the 
adhesive  around  the  spotweld  caused  by  lack  of  adhesive  flow,  as  noted  in  this 
figure.  The  flaw  in  the  weld  may  be  caused  by  the  entrapped  air.  This  is 
typical  of  situations  that  might  possibly  occur  in  weldbonded  structures.  The 
through-crack  in  the  sheet  may  have  originated  as  a part-through  flaw  that  had 
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DE  BON  DING  CRACK 
IN  THE  ADHESIVE 


Figure  21,  WeJdbonded  Structure 

grown  to  a through- the- th ickness  flaw  during  fatigue  cycling.  The  precise 
finite  element  analysis  of  this  type  of  structure  would  require  several  orders 
of  magnitude  of  complication  over  that  previously  discussed,  therefore  certain 
simplifying  assumptions  have  to  be  made.  This  structure  can  be  easily  modeled 
using  a two-dimens ional  finite  element  model  based  on  the  principles  discussed 
in  Paragraph  2 , L , Both  the  spotweld  and  the  adhesive  are  represented  as  shear 
elements.  No  shear  elements  are  provided  in  the  portion  of  the  adhesive  where 
there  is  debonding,  nor  in  the  portion  of  the  weld  where  there  is  a void,  as 
shown  in  Figure  22.  In  this  model,  the  properties  of  the  shear  elements  rep- 
resenting the  spotweld  and  adhesive  will  be  different,  since  the  stiffnesses 
of  the  adhesive  and  spotweld  are  different,  The  shear  elements  are  closely 
spaced  in  the  vicinity  of  the  crack  and  debond.  The  grid  point  interval  for 
shear  elements,  and  consequently,  the  size  of  the  plate  elements  in  both  skin 
and  strap,  are  gradually  increased  away  from  the  crack  and  debond.  The  strap 
is  assumed  to  be  bonded  throughout  its  length  by  having  shear  elements  along 
the  length.  In  a cracked  structure,  the  first  six  rivets  or  welds  adjacent  to 
the  crack  surface  will  be  effective  in  load  transfer  for  crack  lengths  up  to 
about  2,0  inches,  hence  in  the  structure  of  Figure  21,  the  first  six  spotwelds 
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— ADHESIVE  THICKNESS 

t — H 


Figure  22.  Modeling  of  Weldbonded  Structure  Showing  Shear  Elements 
Representing  Spotwelds  and  Adhesive 


would  require  special  modeling.  The  shear  elements  representing  the  ■'.dhesive 
have  the  same  elastic  properties  as  the  adhesive.  The  shear  elements  rep- 
resenting the  rivet  or  weld  are  proportioned  to  have  the  same  deflection  as 
that  of  the  riveted  or  welded  joint  (Reference  1),  Let  A be  the  deflection 
of  a welded  lap  joint  under  unit  load.  In  the  idealized  shear  element  In  the 
finite  element  model,  the  def lection  under  unit  load  is  given  by 


6 = 


h 


A ^ 
s a 


(53) 


where 

the  area  of  the  shear  element 
shear  modulus  of  the  weld  material 

distance  between  the  points  connected  by  the  shear  element 
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This  deflection  of  the  shear  element  should  be  equal  to  the  deflection  of  a 
spotwelded  lap  joint*  or 


A - 6 = 


A ju 
s a 


A 


s 


Aju 

a 

h 


<54) 


This  method  of  modeling  was  used  for  riveted  and  bolted  panels  in  Reference  1 
with  excellent  correlation  obtained  between  analytical  and  experimental  data 
for  crack  openings  and  load  transfer  from  cracked  member  to  uncracked  member. 
It  should  be  noted  that  this  method  of  modeling  will  require  deflection  data 
from  lap  joint  tests. 

Next*  consider  the  rivet  bonded  structure  of  Figure  23,  The  shear  ele- 
ments representing  rivets  will  be  specially  modeled  so  that  the  deflection  of 
the  shear  element  is  the  same  as  the  deflection  of  the  riveted  joint.  As 
discussed  in  Reference  1,  the  deflection  in  riveted  aluminum  alloy  sheets  is 
given  by 


6 


E d 
a 


CS5) 


where 

6 - deflection 
P = applied  load 

E = modulus  of  aluminum 

a 

d = rivet  diameter 
f = constant 

For  aluminum  alloy  rivets* 


f = 5.0  + 0.8 

*a  hi 

where  hi  and  h$  are  the  thicknesses  of  the  joined  sheets. 


(56) 
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PLATE 


Figure  23,  Rivet  Bonded  Structure 


In  the  idealized  shear  element,  the  deflection  is  given  by 


6 


Ph 


V 


a 


where 

= the  area  of  the  shear  element 

V = shear  modulus  of  aluminum 
a 

h = distance  between  the  grid  points  connected  by  the  shear  element 
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Equating  the  two  deflections 


Pf  Ph 


Ed  AM 


a 


s a 


and  solving  for  area  of  the  shear  element * 


If  the  materials  connected  by  the  attachments  have  different  moduli*  the 
empirical  constant  in  the  fastener  deflection  Equation  (55)  is  given  by 


(58) 


where  Ex  and  Eg  are  Young’s  moduli  of  the  two  connected  sheets. 

This  type  of  finite  element  modeling  can  also  be  used  for  brazed  struc- 
tures.  In  order  to  proportion  shear  elements  for  brazed  structures*  it  will 
be  necessary  to  establish  a deflection  equation  of  the  type  represented  by 
Equation  (53),  This  will  have  the  form 


6 = KP 


(59) 


where 


6 = deflection  per  unit  length 
P = applied  load 
K = constant. 


The  constant  K can  be  obtained  from  load  displacement  data  obtained  from 
simple  tests  on  double  lap  joints. 
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SECTION  4 


EFFECT  OF  BENDING 

4. 1 INTRODUCTION  AND  GENERAL  COMMENTS 

The  finite  element  and  mathematical  analyses  of  adhesively  bonded 
structures  discussed  in  Section  2 were  described  for  extensional  type  loading 
where  it  was  assumed  that  neither  the  cracked  layer  nor  the  sound  layer  have 
bending  stiffness.  In  reality,  however,  the  presence  of  a crack  in  one  Layer 
of  a two-ply,  adhesively  bonded  structure  or  a plate  with  a bonded  stiffener 
gives  rise  to  out-of-plane  deformation  due  to  unsymmetry.  This  causes  bending 
in  the  structure.  The  influence  of  this  out-of-plane  bending  will  be  to  in- 
crease the  tensile  stresses  in  the  cracked  layer  and  thus  increase  the  stress 
intensity  factor.  Thus,  the  stress  intensity  factors  obtained  by  two- 
dimensional  analysis  need  to  be  corrected  for  the  influence  of  bending.  In 
this  section,  a method  to  do  this,  based  on  load  transferred  to  the  sound 
layer  of  a two-ply  structure  is  described. 

4.2  DEVELOPMENT  OF  CORRECTION  FOR  CENTER  CRACKED  STRUCTURE, 

Consider  a single  cracked  layer  shown  in  Figure  24a,  The  solution  Lo  this 
problem  is  obtained  by  superposing  the  cases  shown  In  Figures  24b  and  24c. 

From  Figure  24c,  the  crack  has  a uniform  stress  cr®*  The  total  force  at  the 
crack  plane  is  area  x stress  = (2ahx  Ja™,  This  force  is  taken  entirely  by  the 
stress  singularities  ahead  of  the  crack  tips.  Next,  consider  a two-ply  ad- 
hesively bonded  panel  with  a crack  in  one  layer  as  shown  in  Figure  25a.  In 

00 

this  structure,  the  same  force  release  (2aa  hi)  at  the  crack  plane  is  taken 
partly  by  the  cracked  layer  in  the  form  of  stress  singularities  similar  to  a 
single  layer,  and  the  remainder  of  the  force  is  transmitted  to  the  adjacent 
sound  layer  through  the  adhesive.  Let  the  force  transmitted  to  the  sound 
layer  be  2ah!a  (area  x stress)  where  is  the  stress  transferred  to  the 

sound  layer.  The  total  force  acting  on  the  end  of  each  layer  is  shown  in 
Figure  25b,  Due  to  the  force  transferred  to  the  sound  layer,  the  force  in 
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a 


(a)  Cracked  plate 

with  remote  stress 


00 

a 


(b)  No  crack  (c)  Crack  with  uniform 

stress  field 


Figure  24*  Superposition  Technique  for  a Single-Layer  Crack 


each  plate  at  the  plane  of  the  crack  will  be  different  from  that  at  the  ends 
as  shown  in  Figure  25c*  The  net  internal  unbalanced  force  between  the  two 
layers  at  the  crack  plane  is  equal  to  that  which  is  transferred  from  the 
cracked  to  the  sound  layer,  i.elf  2ah^Ot*  This  unbalanced  force  causes  the 
bending  (Figure  25d). 

Assuming  the  thickness  of  adhesive  to  be  small  compared  to  plate  thick- 
ness, the  couple  generated  by  the  load  transfer  between  the  two  plies  is  given 
by 


C = 2ahi  0 1 (h^hg) 

where  is  the  distance  between  the  centerlines  of  the  pLates. 


or 


C — ah^  (hi  ) O' 


(60) 


Assuming  that  the  bending  is  resisted  by  the  entire  plate  width  Wt  the  maximum 
bending  stress  in  the  Layer  is  given  by 

C y 


= 


max 


(61) 
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(a)  Two-ply,  adhesively  bonded  panel  with  a 
Win  a"  WI^ct*  crack  in  one  PLy 


(Wha”-2ahiCJt)  <wl*°  +2ahlCTt) 


(c)  Section 
at  the 
crack 
plane 


Free-body 

diagram 


2ahi  a\ 


'2ahi  cr 


(d)  Net 

forces 
caus ing 
bending 


( b)  Side  view  showing  the 
total  force  at  the  ends 


Figure  25*  Bending  Correction  Eor  Adhesively  Bonded  Panel  with  a Crack  in  One  Ply 
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where  y is  the  distance  of  extreme  fibers  of  the  cracked  plate  from  the 
neutral  axis  and  I is  the  moment  of  inertia  of  the  section.  For  computing 
the  values  of  y and  the  moment  of  inertia  I of  the  section,  the  thick- 
ness  of  the  adhesive  may  be  neglected  and  the  calculations  based  on  the 
initially  uncracked  cross-sec tion  of  both  plies.  Thus,  the  cross-section 
of  Figure  26a  may  be  represented  by  the  one  shown  in  Figure  26b.  For  com- 
puting the  neutral  axis  and  moment  of  inertia  of  the  section  based  on  the 
uncracked  section,  the  equivalent  section  (Reference  10)  is  shown  in  Fig- 
ure 27,  where  it  is  assumed  that  E?  is  greater  than  Elt 
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(a)  Section  of  an  adhesively  bonded  panel 


— W m- 

Es 

SOUND 

Ei 

CRACKED 

(b)  Idealized  section  for  bending 
analysis 

Figure  26.  Section  for  Applying  Bending 
Correction 
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Figure  27,  Equivalent  Section  for  Bending  Stress 

Analysis  (neither  layer  assumed  cracked 
for  calculation) 
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Substituting  Equation  (61)  in  (60),  gives 


ahi  (hi  + h2  )o  t 

ab  I ^max  (62) 

4.2.1  Load  Transfer  Factor 

Define  the  term  load  transfer  factor  with  symbol  M as  the  ratio  of  the 
load  transferred  (to  the  sound  layer)  to  the  load  released  due  to  a crack. 

2ah*at 

M ^ S = — 

2ah^<J  a (63) 


Equation  (62)  reduces  to 


ahi  (hi  + hg  ) 


ct  = 


ao 

May 


max 


03 

a 


ahx  (hi  + hg  ) 
I 


My 


max 


= Be,  Bending  correction  factor  (64) 

The  effective  stress  in  computing  the  stress  intensity  factors  is  increased 

bv  an  amount  c.  . For  the  case  hi  = hp  = h and  - Ep  , y “ h and 
J b * * x max 

I = 1/12  W(2h)3 


00 


a 


3a 

w 


M 


(65) 


4,2.2  The  Concept  of  Effective  Stress 

In  a cracked,  finite  width,  single- layer  structure  (Figure  24a),  the 
stress  intensity  factor  varies  as  a function  of  crack  length  according  to 
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the  equation 


K - o /na  f (a/W) 
s 

where  f(a/W)  is  the  finite  width  correction 


or 


f(a/W) 


K 

s 


Let  Ka  he  the  stress  intensity  factor  in  a crack*  adhesively  bonded 
A 

structure  with  half-crack  length  a and  subjected  to  an  applied  stress  of 

30 

a {Figure  25a).  The  equivalent  stress  acting  remotely  on  a single- 

layer  structure  such  as  shown  in  Figure  24a,  and  giving  the  stress  inten 

sity  factor  K at  half-crack  length  a,  is  given  by 
A 


K = a / na  f (a/W) 
A e 

using  Equation  (66b),  can  be  expressed  as 


or 


a = 


e 


K 


5 


oo 

a 


4.2.3  Interrelation  Between  Load  Transfer  Factor  and  Stress  Intensity 

Factors 

The  stress  transferred  to  the  sound  layer  is  the  difference  between 
the  remotely  applied  stress  and  the  uniform  inplane  stress  that  matches 
the  required  singularity  effect  in  the  cracked  layer,  i.e,. 


a - <J  - a 

t e 


( 66a ) 


(66b) 


(67) 


(68) 
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Using  Equation  (68) 
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GC 

a 


K 
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The  load  transfer  factor  is 


M ~ 


(69) 


(70) 


is  the  stress  intensity  factor  obtained  from  the  finite  element  or  mathe- 
matical analysis  of  an  adhesively  bonded  structure  with  a crack  in  one  layer , 
and  is  the  stress  intensity  factor  for  a single  layer  with  a crack  the 
same  size  as  the  adhesively  bonded  structure,  K is  available  in  the  form 
of  Equation  (66a)  (References  11  and  12)*  The  value  of  M obtained  here, 
is  used  in  computing  bending  correction  in  Equation  (64), 

4,3  EXAMPLE  OF  BENDING  CORRECTION 

As  a bending  correction  example,  consider  a two-ply,  adhesively  bonded 
panel  with  a center  crack*  The  panel  width  is  taken  as  six  inches,  and 
each  layer  has  the  same  modulus  with  a thickness  of  0*063  inch*  The  stress 
intensity  factors  for  this  structure,  assuming  an  elliptical  debond  with 
b/a  - 0*1  (semiminor  to  semimajor  axis  ratio  equal  to  0,1),  are  determined 
by  finite  element  analysis  and  are  show^i  in  Figure  28,  and  Column  2 of 
Table  1*  For  the  structure  under  consideration,  Equation  (65)  reduces  to 


CO 


a 


for  W = 6 inches, 


h 

— = 0 , 5aM 

CO 

u 


Various  steps  involved  in  the  application  of  the  bending 
shown  in  Table  1.  The  value  of  a/W  is  computed  for  each 


correction  are 
crack  length 
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Figure  lb.  Influence  of  Bending  on  Stress  Intensity  Factors 
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and  is  shown  in  Column  3,  For  these  values  of  a/W,  finite  width  correction 
factors  can  be  obtained  from  any  fracture  mechanics  handbook  (References  11 
and  12}  , and  are  shown  in  Column  4,  Load  transfer  factor  M,  given  by  I - — , 
is  computed  and  is  shown  in  Column  5*  Bending  correction  factor  Be  is 
given  by  0,5aM  and  is  shown  in  Column  6,  The  stress  intensity  factors 
corrected  for  bending  are  given  by  (1  + Be)  times  those  obtained  from  the 
two-dimensional  analysis,  and  are  shown  in  CoLumn  7.  The  plot  of  bending 
corrected  stress  intensity  factors  is  shown  in  Figure  28,  It  is  seen  that 
for  a half-crack  length  of  1,5  inches  in  a six- inch  wide  panel,  the  stress 
intensity  factors  may  increase  by  as  much  as  53  percent  due  to  bending. 

These  stress  intensity  factors  corrected  for  bending  have  shown  excellent 
agreement  with  those  obtained  from  experiments.  Also,  using  these  bending 
corrected  stress  intensity  factors,  it  was  possible  to  predict  crack  growth 
life  within  ten  percent  of  actual  life. 


It  may  be  noted  that  the  bending  correction  discussed  above  is  appli- 
cable to  thin  plates.  Also,  it  is  assumed  that  entire  width  of  the  plate  is 
effective  in  resisting  bending.  If  the  panel  width  is  large  in  adhesively 
bonded  structure,  the  influence  of  bending  will  be  localized  in  some  width 
on  either  side  of  the  crack  tips.  Hence,  the  effective  width  resisting 
bending  will  be  smaller  than  the  width  of  the  plate. 


4.4  BENDING  CORRECTION  FOR  CRACKS  AT  STRESS  CONCENTRATIONS 

In  Section  4,2  bending  correction  for  center  crack  panels  was  discussed. 
In  following  subsections  the  application  of  the  method  to  cracks  at  stress 
concentrations  is  discussed. 


4.4,1  Application  of  Analysis  to  Cracks  at  Holes 

The  method  of  bending  correction,  based  on  a load  transfer  factor  M, 
has  a general  application  even  to  cases  where  cracks  are  emanating  from  holes, 
etc.  For  the  case  of  a center-cracked  panel,  the  force  release  was  taxon  as 
ZaCJ^h^  bu  t f or  c ases  where  cracks  are  emanating  from  holes,  this  stress  release 
will  be  different.  For  example,  the  adhesively  bonded  structure  with  a crack 
emanating  from  a hole  as  shown  in  Figure  29,  will  have  a stress  concentration 
around  the  hole,  hence  the  stress  release  will  be  Larger  than  the  center  cracked 
panel.  The  stress  distribution  around  a hole  in  an  uncracked  plate  is  shown  in 
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SECTION  A.- A 


Figure  29.  Adhesively  Bonded  Structure  with  a Crack  at  Hole 
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Figure  30,  the  stress  being  ^ (ot  > cr  ) at  distance  a from  the  edge  of 

the  hole.  In  the  elastic  range,  the  stresses  cr  and  criwil!  be  directly 

max 

proportionaL  to  the  remote  appLied  stress.  The  shaded  area  A shown  in 
the  figure,  multiplied  by  the  thickness  of  the  plate  will  be  the  force  re 
release  due  to  the  presence  of  a crack  (Figure  SOb),  Let  G be  the  average 
stress  concentration  factor  over  the  entire  crack  length.  The  area  of  the 

03 

stress  diagram  in  Figure  30a  will  be  Gaa  . The  load  transferred  to  the 

CO 

sound  Layer  = ah^Gat  ~ ah^  GMct  , The  bending  moment  (couple)  caused  by 
load  transfer  is 


* ( hi  + hs  ) 

C — aGMo  ~ 


(72) 


a 


(a)  Stress  concentration 
at  a hole 


(b)  Single  layer  with  a 
crack  ai  a hole 


Figure  JO.  Crack  at  Stress  Concentration 
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The  stress  due  to  bending  is 


a 


b 


Cy 


max 


The  value  of  M is  given  by  Equation  (70) 


(73) 


M = 


where  is  the  stress  intensity  factor  in  the  adhesively  bonded  structure 

(Figure  29),  and  K is  now  the  stress  intensity  factor  for  a single  layer 
s 

with  a radial  crack  at  a hole  (Figure  30b) 

4,4.2  Example  of  Bending  Correction  for  a Crack  at  Stress  Concentrations 

Consider  the  two-ply  bonded  structure  with  a crack  emanating  from  a 
hole,  as  shown  in  Figure  29,  The  panel  has  a width  of  12  inches  and  a hole 
diameter  of  3/16-inch.  Each  layer  is  7075-T73  aluminum  with  a thickness  of 
0,063-inch,  The  bending  moment  is  given  by  Equation  (72)  as 

. M 

C — a GMa  h^ 


and 


cw 

ab = -r- 


For  the  structure  under  consideration 


hx  = h,,  = h 


y - h 
max 


and 


I = — W(2h)3 
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using  these  values 


3affGM 
b 2W 


for 


W = 12" 


0,  = 0.125  aGM  (74) 

D 

Various  steps  in  applying  bending  correction  are  shown  in  Table  2,  The 

stress  intensity  factor  for  crack  at  a hole  in  a single  layer  (K  } is  obtained 

s 

by  multiplying  (stress  intensity  factor  for  crack  at  a hole  in  a single 
layer  infinite  width)  by  the  finite  width  effects  for  center  crack  paneLs. 

This  finite  width  correction  does  not  account  for  the  presence  of  a hole. 

The  values  of  are  obtained  by  finite  element  analysis.  The  values  of 
G are  computed  from  unknown  solutions  of  stress  distributions  around  holes 
(e.g.  Reference  13). 

Table  1 shows  that  the  influence  of  bending  correction  on  stress 
intensity  factors  is  significant  for  large  crack  lengths.  The  influence 
of  stress  concentration  (G)  on  bending  correction  is  significant  for  small 
crack  lengths  only.  The  total  bending  correction  at  the  small  crack  lengths 
may  be  negligible. 
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TABLE  2.  BENDING  CORRECTION  FOR  A PANEL  WITH  A CRACK  AT  A HOLE 
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3/16-inch  hole 


SECTION  5 


CRITERION  FOR  THE  CRACKING  OF  A SOUND  LAYER 


Under  cyclic  loading  a crack  may  initiate  in  a sound  layer  that  is 
directly  attached  to  a cracked  layer  because  there  is  a large  amount  of 
load  transfer  between  cracked  and  initially  sound  layers.  The  load  trans- 
ferred to  the  sound  layer  causes  high  stress  concentrations  in  the  load 
transfer  regions,  and  subsequent  cycling  initiates  a crack*  The  cracking 
of  an  initially  sound  layer  will  depend  on  the  level  of  the  load  trans- 
ferred and  the  number  of  fatigue  cycles  applied.  As  shown  in  Section  4, 
the  Load  transferred  to  the  sound  layer  will  depend  on  the  crack  length, 

A general  criterion  for  the  cracking  of  a sound  layer,  based  on  load 
transfer  to  the  sound  layer,  is  developed  here*  It  is  assumed  that  the 
initially  cracked  layer  has  a small  flaw  that  propagates  under  fatigue 
loads.  During  the  crack  propagation,  the  sound  layer  is  subjected  to  in- 
creasing fatigue  Loads  in  the  load  transfer  region  as  the  load  transfer 
increases  with  crack  length,  When  the  load  transferred  to  the  sound 
layer  reaches  a critical  value,  the  crack  will  initiate  in  the  sound 
layer.  The  application  of  this  criterion  to  various  structures  is  dis- 
cussed in  the  following  paragraphs, 

5.1  CRACKING  OF  A SOUND  LAYER  IN  CENTER- CRACKED  BONDED  STRUCTURES 

The  concept  of  load  transfer  factor  H was  introduced  in  Section  4 and 
was  defined  as  the  ratio  of  the  load  transferred  to  the  cracked  layer,  to 
the  load  released  due  to  the  presence  of  a crack  (Equation  63),  It  is 
assumed  that  a crack  will  initiate  in  the  sound  layer  when  the  load  trans- 
fer factor  reaches  the  critical  value  M , defined  as  the  ratio  of  the  load 

c. 

transferred  to  the  sound  layer,  to  the  remote  load  in  this  critical  portion 
under  consideration.  Consider  the  center-cracked,  two-ply,  adhesively 
bonded  structure  shown  in  Figure  25a.  In  this  structure,  force  release 
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30 

clue  to  the  presence  of  a crack  is  2 ar  hi  and  the  load  transfer  factor 
associated  with  bending  for  this  structure  was  defined  in  Section  4 as 

2ahia  o 

M = - = -r 

2ah xCT  cr 

In  Section  4,  the  discussion  on  interaction  between  cracked  and  sound  layers 
was  restricted  to  load  transfer  without  any  special  consideration  to  the 
local  distribution  of  stresses.  The  sound  layer  will  be  subjected  to 

X 

stresses  higher  than  a (remote  applied  stress),  on  the  crack  plane,  over 
a width  approximately  equal  to  2a,  with  maximum  stress  expected  to  occur 
in  the  center  of  the  plate.  Thus,  a crack  will  initiate  at  the  center  of 
the  crack.  The  method  of  computing  the  value  of  load  transfer  factor  was 
discussed  in  Section  4 and  calculations  for  a 6-inch  wide  panel  were  shown 
in  Table  1.  A plot  of  load  transfer  factor  M versus  half-crack  length  a for 
this  panel  is  shown  in  Figure  31.  A similar  plot  fox  a 12-inch  wide  center 
cracked  bonded  panel  is  also  shown  in  the  figure. 

The  experimental  data  on  various  center  crack  two-ply  bonded  panels 
showed  crack  initiation  in  the  sound  layer.  The  experimental  values  of  the 
crack  lengths  at  which  the  sound  layer  cracked  for  various  panels  are  shown 
in  Table  3,  The  load  transfer  factor  at  these  crack  lengths  is  obtained 
from  Figure  31  and  is  shown  in  Table  3.  It  is  seen  that  the  value  of  the 
critical  load  transfer  factors  are  fairly  constant. 

Assuming  that  the  crack  initiated  in  the  sound  layer  when  the  half- 
crack length  a in  the  cracked  layer  was  about  0.8  inch,  the  number  of 
cycles  required  to  propagate  a crack  from  a haLf-crack  length  of  0.8  inch 
to  1.4  inch  are  shown  in  Table  3.  These  cycles  may  be  considered  as  the 
cycles  for  crack  transfer  to  the  sound  Layer. 

Also  shown  in  Table  3,  are  the  stress  Intensity  factors  in  the  cracked 
plate  at  which  the  sound  layer  cracked.  It  is  seen  that  except  for  Panel 
11-16.  the  values  of  K are  fairly  constant.  Hence,  it  would  appear  that 
a crack  will  initiate  in  a sound  layer  when  the  crack  length  in  an  initially 
cracked  layer  reaches  a critical  value,  or  at  a critical  stress  intensity 
factor.  This  criterion  for  the  cracking  of  a sound  layer,  based  on  the 
critical  stress  intensity  factor  in  the  crack  layer,  may  be  valid  for  a 
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Load  Transfer  Factor  M versus  Half-Crack  Length 


I’ ABLE  3.  CRITICAL  LOAD  TRANSFER  FACTORS  FOR  THE  CRACKING  OF  A SOUND  LAYER 


SPEC I MEN 

AND  TYPE 

OF  FLAW 

SECOND  PLY 

CRACKING  AT 

a = AND 
LOCATION  ( ) 

( inches ) 

CRITICAL  LOAD 
TRANSFER  FACTOR 

MC 

( inches ) 

N 

CYCLES 

FOR 

CRACK 

TRANSFER 

Kmax** 

Ksi  k'i  n 

Center  Crack 

1-1 

12"  Wide 

1 . 363 

(A) 

0.640 

33,400 

15.870 

Ua  = 0.0125" 

1.433 

(B) 

0.64  5 

37,400 

16.220 

1.339 

CO 

0.635 

34 , 000 

15.730 

1 1-2 

12"  Wide 

1.357 

(A) 

0.645 

33,800 

15.620 

ha  = 0.0125" 

1.295 

(B) 

0.631 

37,500 

15.230 

11-16 

12"  Wide 

1.114 

(A) 

0.605 

12,400 

19.480 

ha  = 0.008" 

1 .270 

(C) 

0.630 

14,100 

20. 560 

II- 1 0 

6n  Wide 

1.353 

(B) 

0.6  57 

36,500 

1 5.980 

1-11 

6M  Wide 

1.382 

(A) 

0.680 

35,000 

1 

1 5 . 890  ' 

ha  = 0 ,01 5" 

1.460 

(C) 

0.685 

34, 500 

16. 180 

Average  M(;  - 

^Stress  intensity  in  the  cracked 


0,645 

layer  when  the  sound  layer  cracked 


particular  geometry,  however,  it  will  not  be  generally  applicable  due  to 
the  fact  that  for  c racks  emanating  from  holes,  there  is  a large  load  trans- 
fer to  the  sound  Layer  near  the  holes  even  though  the  corresponding  stress 
intensity  factors  may  not  be  high, 

5. 2 CRACKING  OF  A SOUND  LAVER  IN  AN  ADHESIVELY  BONDED  STRUCTURE  WITH 

A CRACK  AT  A HOLE 

Consider  the  case  of  a crack  emanating  from  a hole  in  a two-ply  bonded 
structure  as  shown  in  Figure  2*),  In  such  a structure,  the  presence  of  a hole 
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causes  stress  concentration  with  a maximum  value  at  the  edge  of  the  hole,  as 
shown  in  the  single  layer  of  Figure  30a*  Hence,  in  a bonded  structure  with 
a crack  at  a hole,  maximum  load  transfer  to  the  sound  layer  will  be  at  the 
edge  of  the  hole.  Thus,  a crack  in  the  sound  layer  will  initiate  at  the 
edge  of  the  hole  when  the  load  transfer  factor  there  reaches  a critical 
value*  In  computing  the  critical  load  transfer  factor,  some  arbitrary  width 
of  plate  € around  the  hole  has  to  be  defined* 

For  computing  the  critical  load  transfer  factor,  it  is  assumed  that  a 
plate  width  of  e = 0,1  inch  near  the  hole  is  effective,  as  shown  in  Figure  32* 
For  an  uncracked  plate  with  a crack  at  a hole,  let  the  shaded  area  (over  0,1 
inch  width,  Figure  32,  of  the  plate)  be  given  by  0.1  Her  t h being  the  average 


CO  00 

a a 


cr 


Area  = cr  (0. 1H) 

H = EFFECTIVE  STRESS 
CONCENTRATION 
FACTOR 


Figure  32,  Effective  Width  for  the  Critical  Load 
Transfer  Factor 
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stress  concentration  factor  over  0.1  inch  width.  The  load  released  is 

oc? 

equal  to  O.Lh1Ha  and  the  load  transferred  to  the  sound  layer  equals 

CC 

hat(0,l)h1  = 0.1h1HMa  . The  local  load  transfer  factor  associated  with 
crack  transfer  is  M* 

0.1HMo\ 

M*  = 

0 . iaC°h1 

= HM 

- Average  stress  concentration 
factor  x load  transfer  factor 

The  average  stress  concentration  factor  H for  a plate  with  a 3/16-inch 
diameter  hole  over  0,1-inch  wide,  is  1.693 

5.3  EXAMPLE  OF  THE  CRACKING  OF  A SOUND  LAYER  IN  A PANEL  WITH  A CRACK 

AT  A HOLE 

The  load  transfer  factor  M for  a two-ply,  12-inch  bonded  structure 
with  a crack  at  a central  hole,  was  computed  in  Section  3 and  shown  in 
Table  2,  Using  these  values,  local  load  transfer  factor  M*t  correspond- 
ing to  e = 0.1  (H  “ 1.693),  was  computed.  The  plot  of  the  local  load 
transfer  factor  M“  as  a function  of  crack  length  a,  is  shown  in  Figure  33, 

The  figure  also  shows  the  plot  of  local  load  transfer  factor  for  a 
12-inch  wide  panel  with  a crack  at  an  eccentric  hole. 

The  value  of  critical  load  transfer  factor  at  the  cracking  of  the 
sound  layer  was  observed  to  be  0.645,  from  Table  3 for  center-cracked  panels. 
Using  this  value  of  the  crack  length  at  which  a crack  will  initiate  in 

the  sound  layer  for  the  panel  with  a crack  at  a hole,  is  obtained  from  Fig- 
ure 33.  These  crack  lengths  are  0.415  inch  for  the  craek-a t-a-centra  1-hole 
panel  and  0,385  inch  for  the  crack-at-an-eccentr ic-hole  panel.  Table  4 
shows  the  comparison  of  experimental  and  predicted  crack  lengths  at  which 
the  cracking  of  the  sound  layer  takes  place  in  various  crack-at-a-hole 
panels.  It  is  seen  that  the  correlation  between  predicted  and  experimental 
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TABLE  4.  PREDICTED  CRACK  LENGTHS  FOR  CRACKING  OF  A SOUND  LAYER 


SPECIMEN 

AND  TYPE 

OF  FLAW 

SECOND  PLY 
CRACKING  AT 
a = AND 
LOCATION  ( ) 

( inches) 

PREDICTED  ’a' 

FOR  SECOND  PLY 
CRACKING 

(inches) 

K ** 
max 

Ks  i /Tn 

Crack  at  a 
central  hole 

1-5 

12n  wide 

0.370 

(A) 

0.415 

9.924 

h = 0.007" 

0,  395 

(B) 

0.415 

9.900 

a 

0.360 

(c) 

0.415 

9.957 

11-22 

12"  wide 

0.528 

(A) 

0.415 

9.198 

h = 0,007" 

0.420 

(B) 

0.415 

9.935 

a 

0.403 

(C) 

0.415 

9.910 

1-24 

6"  wide 

0.495 

(A) 

0.415 

10.488 

h = 0.01" 

0,395 

(B) 

0.415 

10.070 

a 

0.423 

(C) 

0.415 

10.280 

Average 

0.421 

0.415 

9.962 

Crack  at  an 
eccentric  hole 

11-12 

12"  wide 

0.364 

(A) 

0.385 

9.843 

h = 0.009" 

0.328 

(B) 

0.385 

9.827 

a 

0.384 

(C) 

0.  385 

9.935 

Average 

0.359 

0.385 

9.868 

**S tress  intensity  in  cracked  Layer  at  sound  layer  cracking. 
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eccentric  hole 


w 

o 

x: 

u 

n 


-c 

Ai 

&0 

c 

a. 

nJ 

jt* 

c 

cfl 

c. 

u 


C4} 

*rH 

Lu 


78 


Local  Load  Transfer  Factor  M*'  versus  Crack  Length 


vaLues  is  good.  The  table  also  shows  that  the  stress  intensity  factors 
in  the  cracked  layer  at  the  cracking  of  the  sound  layer  arc  fairly  con- 
stant, However*  these  values  differ  significantly  from  those  for  the 
case  of  the  center-cracked  panel*  Hence*  the  cracking  of  a sound  layer 
could  not  be  based  on  the  critical  load  transfer  factor  in  a cracked  layer* 

The  predicted  values  of  crack  lengths  at  the  cracking  of  a sound  layer 
will  be  affected  by  the  assumed  plate  width  € for  the  local  load  transfer 
factor.  The  value  of  average  stress  concentration  factor  H will  depend  on 
the  value  of  e,  For  e = 0.05  and  0.15,  the  values  of  H are  2.047  and  1,533, 
respectively.  The  values  of  local  load  transfer  factor  M for  e = 0,05, 

0*1,  and  0.15  are  shown  in  Table  5 for  a 12- inch  wide  panel  with  a crack 
at  a central  hole.  The  plot  of  local  load  transfer  factor  H'  as  a function 
of  crack  length  a for  the  three  values  of  € is  shown  in  Figure  34,  Using 
the  value  of  critical  local  transfer  factor  as  0*645,  the  predicted  crack 
lengths  for  the  cracking  of  a sound  layer  are  shown  in  the  Figure  33.  The 
predicted  crack  lengths  may  vary  between  0.255  inch  and  0,485  inch  for  e 
between  0,05  inch  and  0,15  inch. 


TABLE  5.  LOCAL  LOAD  TRANSFER  FACTOR  FOR  A 12- INCH  WIDE  PANEL  WITH  A 
CRACK  AT  A CENTRAL  HOLE 


a+r^ 

[ inch) 

a 

{ inch) 

a 

r 

M=l-  ^ 

M*  = MH 
= 2.047M 
(e  = 0.05") 

M*  = MH 
= 1.693M 

(E  = 0.1") 

M*  = MH 
= 1.533M 
(e  = 0.15") 

0.2 

0.1062 

1.13 

0.2575 

0.5271 

0.4360 

0.3947 

0.5 

0.4062 

4.33 

0.3656 

0.7484 

0.6190 

0.5605 

1,0 

0.9062 

9.66 

0.5204 

1.0653 

0.8810 

0. 7978 

1.5 

1.4062 

14.99 

0.6188 

1.2667 

1.0476 

0.9486 

2.0 

1.9062 

20.30 

0.6765 

1.3848 

1.1454 

1.0371 

3.0 

2.9062 

31.00 

0.7412 

1.5172 

1.2550 

1.1363 

*3/ 16- inch  hole 
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Figure  34.  Local  Load  Transfer  Factor  M"  versus  Crack  Length  a for  Various  Values  of 


5. 4 LIMITATIONS  OF  CRACK  TRANSFER  CRITERIA 


The  crack  transfer  criteria  based  on  critical  load  transfer  factor  Mf 
are  applicable  to  cases  where  the  cracked  layer  and  sound  layer  are  approxi- 
mately the  same  thickness.  The  value  of  may  depend  on  the  material  and 
thickness  of  the  sound  layer.  The  thickness  and  the  modulus  of  the  sound 
layer  material  are  taken  into  consideration  in  the  analysis*  However , the 
metallurgical  factors  will  influence  the  initiation  of  a crack  in  the  sound 
layer,  hence  the  value  of  M may  be  different  for  different  materials  and 
thicknesses. 

All  but  one  test  was  conducted  at  same  maximum  stress.  In  Table  3, 
panel  11-16  had  a maximum  stress  about  30  percent  higher  than  that  for  the 
other  panels.  The  value  of  for  this  panel  is  still  within  the  experi- 
mental data  scatter.  The  value  of  maximum  stress  may  have  some  effect  on 
the  value  of  critical  load  transfer  factor  M^, 
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SECTION  6 


ANALYSIS  OF  DEBQND  PROP  AG  AT  ION 


Cracked,  adhesively  bonded  structures  subjected  to  cyclic  loading 
exhibit  debonding  in  the  adhesive  layer  around  the  crack.  The  debond  propa- 
gates as  the  crack  in  the  metallic  layer  propagates.  The  rate  of  debond 
propagation  will  depend  on  the  type  of  adhesive  and  load  transferred  to  the 
sound  layer  through  the  adhesive.  The  load  transferred  through  the  adhe- 
sive will  depend  on  the  properties  and  thicknesses  of  the  adhesive  and 
adherend  materials  and  the  applied  loads.  The  propagation  of  the  debond 
will  be  governed  by  the  criterion  connected  with  the  failure  of  the  adhe- 
sive, For  materials  exhibiting  plastic  deformation,  a simple  criterion 
for  failure  of  the  adhesive  could  be  the  critical  strain.  It  may  be  assumed 
that  the  adhesive  has  failed  if  the  strain  in  the  adhesive  reaches  a criti- 
cal value.  This  critical  value  of  strain  may  be  taken  as  the  strain  at  fail 
ure  in  a static  test  on  the  adhesive.  Based  on  these  assumptions,  a crite- 
rion for  the  propagation  of  a debond  was  developed  in  this  research  program. 
The  details  of  this  criterion  for  debond  propagation  and  its  application  to 
various  crack  geometries  is  discussed  in  the  following  paragraphs. 

6,1  CRITERION  FOR  THE  PROPAGATION  OF  A DEBOND 

Consider  a single-layer  plate  with  a crack,  as  shown  in  Figure  35a. 
Under  applied  loading  the  crack  will  open,  as  shown  in  Figure  35b-  Next, 
consider  a case  where  this  plate  is  the  cracked  layer  in  a two- Layer,  adhe- 
sively bonded  structure  (Figure  1),  In  such  a structure,  the  crack  open- 
ing will  be  smaller  than  the  single- layer  opening,  as  shown  in  Figure  35c, 
The  reduction  in  the  crack  opening  will  depend  on  the  load  transferred  to 
the  uncracked  layer,  which  is  a function  of  the  thickness  and  material 
properties  of  the  adhesive  and  adherends,  Also,  the  presence  of  a crack 
in  one  layer  of  a bonded  structure  gives  rise  to  out-of-plane  bending. 

This  load  transfer  to  the  sound  layer  and  the  out-of-plane  bending  will 
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(a)  Single-layer  Plate 
with  a Crack 


(b)  Crack  Opening  in  a 
Single- layer  Plate 


(c)  Crack  Opening  in  a 
Two-ply , Adhesively 
Bonded  Structure 


ADHESIVE 


(d)  Cross-section  of  Adhesive 
in  a Two-ply , Cracked, 
Bonded  Structure 


Figure  35.  Crack  Openings  in  Cracked  Structures 
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create  normal  and  shear  stresses  in  the  adhesive.  If  both  layers  had 
cracks  of  the  same  size,  the  crack  openings  in  both  layers  would  be  the 
same,  and  there  would  be  no  shear  or  normal  stresses  in  the  adhesive. 

Consider  the  cross-section  of  the  two-ply,  bonded  structure  shown  in 
Figure  35d*  The  displacements  in  the  cracked  and  sound  layer  far  from  the 
crack  plane  will  be  the  same,  hence  no  load  transfer  through  the  adhesive. 

In  the  vicinity  of  the  crack  plane,  the  cracked  sheet  will  undergo  large 
deformation  compared  to  the  sound  layer,  with  the  result  that  the  adhesive 
on  the  surface  of  the  cracked  plate  will  deform  more  in  the  y direction 
than  the  adhesive  on  the  surface  of  the  sound  layer.  The  difference  in  the 
deformation  of  the  two  plates  will  depend  on  the  distance  from  the  crack 
plane.  The  adhesive  in  the  bonded  structure  may  be  looked  on  as  a discrete 
element  (both  in  the  x and  y directions)  connecting  two  plates,  as  shown  in 
Figure  36a,  These  elements  are  denoted  by  Al  , Ag  Bg  , A^  B^  , etc*  (at 
x = 0),  and  Eg  . . . , , . at  other  x locations,  as  shown  in  Fig- 

ure 36c*  Under  applied  loads,  the  elements  will  undergo  deformation  (Fig- 
ure 36b).  The  points  Aa,B1  , Ag,Bg  . , , will  displace  to  A\  , , A*s  t 

B 1 g , * • , respectively,  as  shown  in  Figure  36b,  which  also  shows  that  the 
displacements  are  larger  in  a cracked  plate.  Point  A1  displaces  VCi 
% is  displaced  by  Vsi*  The  relative  displacement  of  the  two  points  will 
be  VC1  - VS1  . 

The  strain  Gj  in  the  adhesive  (element  Ai  at  x = 0),  is  given  by 

VC1  - vsl 

Gj  = (at  y = 0)  (76) 

a 


Similarly,  strain  eg,  g3  , £,  , in  elements  Ag  Bg  , A3B3,  and  A4B4  (at  x = 0) 
is  given  by 


e2 


H 


vcs  - vS3 


Vc3  - VS3 


VC*  - VS4 


h 

a 


(y  = yi  ) 
(y  = ys  ) 
(y  - ya ) 


(77) 
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a>  m 
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(a)  Cross-section  showing  (b)  Displacements  in  plates  (upper  HaLf  only) 

representation  of  adhe- 
sive at  x = 0 


h Pi  h A Fi 


BY  SHEAR  ELEMENTS 

(c)  Sectional  plan  at  y = 0 showing 
representation  of  adhesive 


h = Adhesive  thickness 
a 


Figure  36.  Discrete  Element  Representation  of  Adhesive 
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In  general,  the  strain  e,  at  any  location  is  given  by 


V - V 


a 


The  value  of  the  strain  will  decrease  as  the  distance  from  the  crack 
plane  increases.  The  strains  at  other  x locations  (elements  Eh  Ex  , ,,, Fig- 
ure 36c)  can  be  obtained  in  a similar  manner.  The  adhesive  will  fail  if 
the  strain  in  the  adhesive  reaches  a critical  value  which  is  defined 
as  the  resistance  of  the  adhesive  to  fracture,  is  taken  as  the  failure 

strain  of  the  adhesive  from  the  tensile  test  on  the  adhesive. 

The  displacements  and  in  cracked  and  sound  layers  are  obtained 
from  the  mathematical  or  finite  element  analysis  of  the  cracked w adhesively 
bonded  structures  at  a particular  crack  length  a^.  The  strain  in  the  ad- 
hesive is  obtained  from  these  displacements  using  Equation  78,  The  strain 
in  the  adhesive  is  obtained  from  the  displacements  at  various  (x,  y)  loca- 
tions and  from  this,  lines  of  constant  strain  £ in  the  adhesive  can  be 
plotted.  These  lines  are  schematically  shown  in  Figure  37,  The  debond 
shape  and  size  will  be  given  by  the  curve  corresponding  to  critical  strain 
e , The  predicted  shape  of  the  debond  in  an  adhesively  bonded  structure 
corresponds  to  a particular  crack  length  in  the  cracked  layer. 


Figure  37-  Lines  of  Constant  Strain  in  Adhesive  for  a 
Fixed  Applied  Stress  in  a Bonded  Structure 
(one-quarter  of  contours  shown) 
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As  the  crack  propagates  under  fatigue  loads,  the  crack  opening  and  load 
transfer  to  the  sound  layer  will  increase,  hence  the  relative  displacements 
between  the  two  layers  will  increase,  or  the  strain  in  the  adhesive  at  vari- 
ous x,  y locations  will  continue  to  increase.  As  soon  as  the  strain  at  any 
(x,  y)  location  in  the  bonded  region  reaches  the  critical  strain  e ^ , the 
adhesive  will  fail,  hence  the  size  of  the  debond  will  increase.  Thus,  the 
size  of  debond  can  be  determined  for  every  crack  length  as  the  crack  in  L tie 
cracked  layer  propagates.  From  the  above  criterion  of  debond  propagation, 
the  following  points  are  observed. 

1.  The  size  of  the  debond  is  dependent  on  the  adhesive  thickness  and 
properties*  From  Equation  (78),  it  is  seen  that  if  the  thickness 
of  the  adhesive  is  large,  the  strain  in  the  adhesive  will  be  small, 
hence  the  size  of  the  debond  will  be  small*  If  the  adhesive  thick- 
ness is  very  large,  the  size  of  the  debond  may  be  zero. 

2.  The  size  of  the  debond  will  depend  on  the  applied  loads  (a 

r max 

for  fatigue  loads)*  The  large  applied  loads  will  cause  large 
crack  openings,  hence  large  shear  strain  in  the  adhesive. 

3.  The  size  of  the  debond  will  depend  on  the  relative  clastic 
properties  of  the  sound  and  cracked  layer,  as  these  will  influ- 
ence the  relative  displacement  between  the  two  layers. 

6.2  PROCEDURE  FOR  OBTAINING  CONSTANT  STRAIN  LINES 

The  size  of  debond  is  determined  by  constant  strain  lines.  The 
following  analytical  steps  are  observed  to  obtain  constant  strain  lines. 

1*  The  finite  element  or  mathematical  analysis  of  the  cracked  bonded 
structure  under  consideration  is  carried  out  assuming  no  debond 
(if  debond  size  not  known),  or  assuming  known  debond, 

2*  The  displacements  in  the  y direction,  at  various  locations  around 
the  crack,  in  the  cracked  and  sound  layer  are  obtained  at  the 
known  applied  load.  Using  these  displacements  and  Equal  ion  (78), 
the  strain  in  the  adhesive  at  various  locations  is  obtained* 
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3 * 


The  strain  in  the  adhesive  can  be  plotted  (corresponding  to  any 
applied  load)  as  a function  of  distance  x from  the  centerline  of 
the  crack  for  various  distances  from  the  y plane.  These  plots  are 
schematically  shown  in  Figure  38.  The  figure  also  shows  the  line 
corresponding  to  strain  €1  . It  is  seen  that  Line  y = 0 cuts  line 
Ci  at  x = Xq  , hence  on  the  y = 0 plane,  the  strain  in  the  adhesive 
is  larger  than  e:  , up  to  x = x0.  For  x greater  than  x0,  the 
strain  in  the  adhesive  at  y - 0 is  smaller  than  €i  * Thus,  one 
end  of  the  constant  strain  line  will  be  at  y = 0 and  x = x0,  At 
x - 0,  the  line  Ci  , cuts  the  y - y3  line,  hence  for  y greater  than 
y3  , the  strain  in  the  adhesive  is  less  than  a e*  , thus  the  other 
end  of  the  constant  strain  line  is  at  x = 0 and  y = y3  . The  other 
boundaries  of  the  constant  strain  line  (Figure  38)  are,  at  x = 
y - y1  , and  at  x = ife  , y - yg  . Thus  the  shape  of  the  constant 
strain  line  corresponding  to  ei  , is  obtained.  Similar  lines  can 
be  obtained  for  other  € values.  The  lines  of  constant  strain  are 
shown  in  Figure  37, 

4.  The  size  of  debond  is  given  by  constant  strain  line  having  a 
strain  equal  to  critical  strain  £^. 

5,  If  the  analysis  was  carried  out  assuming  no  debond  in  the  adhesive 
and  the  predicted  debond  size  is  large,  the  finite  element  or 
mathematical  analysis  is  carried  out  again  assuming  the  predicted 
debond  size  of  Step  4 in  the  adhesive*  Using  this  new  analysis, 
the  size  of  the  debond  is  predicted  again.  Thus,  an  iteration 
process  will  be  necessary  to  compute  the  debond  size.  However, 
if  the  predicted  debond  size  in  Step  4 is  small,  the  iteration 
process  need  not  be  done  and  the  predicted  debond  size  (based  on 
assuming  no  debond  in  the  adhesive)  will  be  close  to  actual  debond 
size , 

6,3  DEBOND  PREDICTION  IN  A TWO- LAYER,  ADHESIVELY  BONDED  STRUCTURE  WITH  A 

CENTER  CRACK 

A two-ply,  adhesively  bonded  structure  12  inches  wide,  with  a center 
crack  (Figure  l)  is  considered.  The  structure  has  a half-crack  length  of 
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Shear  Strain  in  Adhesive 


Figure  38*  Schematic  Variation  of  Adhesive  Strain  with  Distance  from 
Centerline  of  Crack  for  Various  Values  of  y {one-quarter 
of  the  panel  considered) 
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one  inch  and  the  maximum  applied  stress  is  15,500  psi.  From  the  finite 
element  analysis  (assuming  no  initial  debond  in  the  adhesive),  the  dis- 
placements in  the  y direction  at  various  grid  points  (x,  y locations)  are 
known.  The  shear  strain  in  the  adhesive  will  be  given  by  Equation  (78). 

The  thickness  of  the  adhesive  layer  for  the  panel  is  taken  as  0.008  inch. 

The  values  of  strain  are  computed  at  various  (x,  y)  locations  around  the 
crack  and  plotted  as  a function  of  x for  various  values  of  y,  as  shown  in 
Figure  39,  The  value  of  critical  strain  for  adhesive  FM-73  is  approxi- 
mately 0,04  (Reference  14)  and  is  also  shown  in  Figure  39.  From  the  figure, 
at  x 3 0,  all  points  greater  than  y = 0,11  (by  interpolation)  will  have  a 
strain  less  than  a critical  strain,  thus  the  boundary  of  debond  at  x = 0 is 
y = 0.11.  Similarly,  other  points  on  the  boundary  of  the  debond  are  given 
by 

y = 0.10  at  x - 0,60 

y = 0,05  at  x = 0,82 

y = 0 at  x - 0.95 

The  predicted  shape  of  the  debond  (one— quarter)  is  shown  in  Figure  40a 

by  the  area  between  Curve  A and  the  x and  y axes.  It  may  be  noted  that  the 
shape  of  the  debond  is  approximately  elliptical  and  similar  to  the  shape  of 
the  crack  openings, 

6.4  DEBOND  PREDICTION  IN  A CRACKED  PLATE  WITH  AN  ADHESIVELY  BONDED  STRINGER 

Consider  the  cracked  plate  with  an  adhesively  bonded  stringer  (Fig- 
ure 14).  For  this  structure,  the  crack  surface  openings,  as  a function  of 
the  distance  from  the  centerline  of  the  crack  obtained  from  finite  element 
analysis  assuming  no  debond,  are  plotted  in  Figure  41  for  various  crack 
lengths.  For  crack  lengths  smaller  than  the  width  of  the  stiffener,  the 
shape  of  the  crack  surface  displacement  curve  is  elliptical,  hence  the  de- 
bond shape  will  be  similar  to  the  case  of  the  center  crack  in  a two-ply, 
adhesively  bonded  panel.  The  predicted  shape  for  a - 0,75  is  shown  in 
Figure  40b,  For  crack  lengths  larger  than  the  width  of  the  stiffener,  the 
crack  surface  openings  show  an  increasing  trend  toward  the  edge  of  the 
stiffener,  as  shown  in  Figure  41,  The  shape  of  the  debond  will  be  different 
in  this  case.  The  relative  displacement  at  various  (x,  y)  locations  are 
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Shear  Strain  in  Adhesive  (in/in) 


Half-Crack  Length  = 1.0  inch 


00 


Figure  39.  Variation  of  the  Shear  Strain  in  the  Adhesive  with  the 
Distance  from  the  Centerline  of  the  Crack,  for  a Two- 
Ply,  Bonded  Structure,  with  a Center  Crack  (one-quarter 
of  the  panel  considered) 
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Figure  41.  Variation  of  Crack  Opening  with  Distance  from  Center  of  Crack 
(one-quarter  of  panel  considered) 


computed  from  the  finite  element  analysis  and  from  these  locations,  the 
strain  in  the  adhesive  is  computed  from  Equation  (78),  The  plot  of  strain 
as  a function  of  the  distance  from  the  centerline  of  the  crack  for  various 
values  of  y and  a half-crack  length  equal  to  2.0  inches  is  shown  in  Fig- 
ure 42.  The  critical  strain  fi  = 0,04  (FM-73)  is  also  shown.  At  x = 0, 
for  the  values  of  y larger  than  Q.l5t  the  strains  in  the  adhesive  are 
smaller  than  the  critical  strain.  Similarly,  at  x = 1,0,  the  strains  at 
points  farther  than  y “ 0.4  are  smaller  than  the  critical  strain.  Other 
points  on  the  boundary  of  the  debond  are  obtained  by  interpolating  between 
intermediate  y values.  The  shape  of  the  debond  is  plotted  in  Figure  39c. 
The  actual  debond  for  a one-quarter  panel,  will  be  the  area  enclosed  be- 
tween the  x and  y axes,  Curve  B,  and  the  edge  of  the  stiffener  (x  - 1*0). 

For  adhesive  FM-400,  the  critical  strain  is  about  0,02  (Reference  14), 
and  is  plotted  in  Figure  42.  At  x - 0,  the  points  farther  than  y = 0.3 
(approximately)  will  have  strain  smaller  than  0.02.  Similarly,  the  boun- 
dary of  the  debond  is  determined  for  other  x locations  and  is  shown  in 
Figure  40c,  It  is  seen  that  the  size  of  the  debond  is  considerably  larger 
than  that  for  the  FM-73  adhesive. 

6.5  DEBOND  PREDICTION  IN  A TWO- FLY,  ADHESIVELY  BONDED  STRUCTURE  WITH 

A CRACK  AT  A CENTRAL  HOLE 

Consider  the  two-ply,  bonded  structure  with  a crack  at  a central  hole 
shown  in  Figure  43.  This  structure  has  a width  of  six  inches  and  an  adhes- 
ive thickness  of  0,008  inch.  The  analysis  of  the  structure  was  carried  out 
assuming  no  debond  in  the  adhesive.  The  strains  in  the  adhesive  are  com- 
puted from  the  displacement  in  the  cracked  and  uncracked  layers,  and  the 
plot  of  shear  strain  as  a function  of  the  distance  from  the  centerline  of 
the  crack  is  shown  in  Figure  44*  It  is  seen  that  the  curves  for  y larger 
than  0.0  show  an  upward  trend  due  to  the  presence  of  a hole.  The  predicted 
shape  of  debond  size,  corresponding  to  FM-73  (e  = 0,04)  is  denoted  in  Fig- 
ure  40d.  It  is  seen  that  the  debond  is  almost  in  the  form  of  a strip. 
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Strain  in  the  adhesive  (in/in) 


Distance  from  centerline  of  crack  (inches) 


Figure  42.  Variation  of  Strain  in  Adhesive  with  Distance  from  Centerline 
of  Crack  for  a Cracked  Plate  with  Bonded  Stiffener 
(one-quarter  of  panel  considered) 
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Figure  43,  Two-Plyf  Adhesively  Bonded  Structure 
with  a Crack  Emanating  from  a Hole 
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Strain  in  adhesive  (in/ in) 


Figure  44.  Variation  of  Strain  in  Adhesive  with  Distance  from  Edge 
of  the  Hole  for  a Two-Ply,  Bonded  Panel,  with  a Crack 
Emanating  from  a Hole  (one-quarter  of  panel  considered) 
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It  may  be  noted  that  for  predicting  debond  sizes  here,  the  finite 
element  analysis  for  various  structures  was  carried  out  assuming  no  initial 
debond  in  the  adhesive-  The  displacements  in  the  cracked  layer  and  sound 
layer  will  be  slightly  different  due  to  the  presence  of  a debond,  and  these 
will  influence  the  strain  in  the  adhesive-  Hence,  the  size  of  the  predicted 
debond  will  be  slightly  different  if  the  initial  debond  present  in  the  ad- 
hesive is  taken  into  consideration.  However,  this  influence  will  be  neg- 
ligible for  ductile  adhesives  where  the  size  of  the  debond  is  small.  For 

brittle  adhesives  having  a low  value  of  e , the  size  of  the  initial  debond 

K 

may  have  a significant  influence  on  the  debond  propagation.  This  can  be 
taken  into  consideration  in  the  analysis,  by  assuming  an  initial  debond  if 
known  or  using  the  iteration  procedure  discussed  in  Paragraph  6.2, 
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SECTION  7 


GUIDELINES  FOR  APPLYING  METHODOLOGY 


In  Sections  2 through  6,  analytical  techniques  were  developed  to  study 
crack  problems  in  adhesively  bonded  structures.  In  this  section  guidelines 
are  provided  for  applying  the  methodology.  The  steps  listed  below  are  to 
be  followed  in  using  the  techniques  developed  here. 

1,  Decide  the  method  of  analysis  to  be  used.  Mathematical  methods 
(computer  programs  given  in  the  Appendices)  are  recommended  for 
parametric  studies.  Finite  element  analysis  is  recommended  for 
crack  growth  analysis. 

2,  The  analysis  of  the  structure  is  carried  out  assuming  no  debond  (it 
debond  size  not  known)  or  known  debond. 

3,  Using  the  displacements  in  cracked  and  sound  layers  and  the  failure 
strain  of  the  adhesive,  the  size  of  the  debond  is  determined  by  the 
procedure  outlined  in  Section  6. 

4,  The  analysis  of  the  bonded  structure  is  carried  out  once  again 
assuming  the  debond  size  determined  in  Step  3.  If  the  debond 
size,  obtained  in  Step  3,  is  small  (say  elliptical  with  a minor- 
to-major  axis  ratio  less  than  about  0,l)  the  debond  size  need  not 
be  recomputed.  If  the  debond  size  in  Step  3 is  large  (or  signifi- 
cantly larger  than  known  debond  size)  the  debond  size  is  recom- 
puted with  the  analysis  of  this  step.  The  bonded  structure  is 
analyzed  once  again  with  the  new  debond  size. 

5,  Compute  the  stress  intensity  factors  using  the  analysis  of  Step  4. 

6,  Apply  bending  correction  to  the  stress  intensity  factors  using  the 
method  outlined  in  Section  4* 

7,  Using  the  method  described  in  Section  3,  determine  the  load  transfer 
factor  or  local  load  transfer  factor  (for  cracks  at  stress  concen- 
tration), If  the  load  transfer  factor  is  less  than  M (the  critical 
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load  transfer  factor  of  the  sound  layer  material),  the  uncracked 
backup  layer  has  not  cracked*  The  stress  intensity  factors  ob- 
tained in  Step  6 are  correct* 

8.  The  stress  intensity  factors  are  determined  for  various  crack 
lengths  using  the  procedure  outlined  in  the  above  steps, 

9.  Obtain  the  crack  growth  data  on  a single  monolithic  layer.  Using 
these  data  and  a crack  growth  equation,  say  ^ = CAK11,  the  constants 
C and  n are  obtained, 

10.  Using  the  stress  intensity  factors  of  Steps  S and  crack  growth 
constants  C and  n,  the  crack  growth  life  of  a bonded  structure  is 
obtained. 

11.  The  crack  growth  Life  in  Step  10  is  predicted  only  up  to  crack 
length  a^,  at  which  the  sound  layer  cracks.  For  crack  lengths 
beyond  sound  layer  cracking,  the  analysis  can  be  carried  out 

by  assuming  crack  length  a^  in  the  initially  cracked  layer,  and 
a crack  in  the  initially  sound  layer. 


102 


REFERENCES 


1,  Ratwani,  M.M.  and  Wilhem,  D.P,,  "Development  and  Evaluation  of  Methods 
of  Plane  Stress  Fracture  Analysis  * A Technique  for  Predicting  Residual 
Strength  of  Structures,"  AFFDL-TR-73-3  Part  II,  Volume  I,  April  1975, 

2.  Erdogan,  F,  and  Arin  K. , "A  Sandwich  Plate  with  a Part-Through  and  a 
Debonding  Crack,"  Engineering  Fracture  Mechanics,  Volume  4,  1972  pp  449-458, 

3*  Arin,  K,  , ,rA  Plate  with  a Crack,  Stiffened  by  a Partially  Debonded 

Stringer,"  Engineering  Fracture  Mechanics,  Vol.  6,  1974,  pp  133-140, 

4.  Arin,  K. , "A  Note  on  the  Effect  of  Lateral  Bending  Stiffness  of  Stringers 
Attached  to  a Plate  with  a Crack,"  Engineering  Fracture  Mechanics, 

Volume  7,  1975,  pp  173-179, 

5,  Alturi,  $.,  Georgia  Institute  of  Technology,  Atlanta,  Georgia,  A Private 
Communication,  May  1976. 

6*  Erdogan,  F.,  and  Ratwani,  M,M,  , "Stress  Distribution  in  Bonded  Joints," 
Journal  of  Composite  Materials,  July  1971,  pp,  378-393, 

7,  Tong,  P, , Fian,  T.H.H.,  and  Lasry,  S,J,,  "A  Hybrid  Element  Approach  to 
Crack  Problems  in  Plane  Elasticity,"  International  Journal  of  Numerical 
Methods  in  Engineering , Volume  7,  1973,  pp,  297-308. 

8,  Erdogan,  F.,  "On  the  Stress  Distribution  in  Plates  with  CoIIinear  Cuts 
Under  Arbitrary  Loading,"  Proceedings  of  4th  U,S,  Nstional  Congress  of 
Applied  Mechanics,  ASME,  1962,  pp,  547-553, 

9,  Muskhellishvilli , N.I,,"  Some  Basic  Problems  of  the  Mathematical  Theory 
of  Elasticity,"  Nordhoff  Publication,  1953, 

10,  Higdon,  A,,  Ohlsen,  E.H. , and  Stiles,  W.B,,  "Mechanics  of  Materials," 

John  Wiley  & Sons  Inc.,  New  York,  1967. 

11,  Sih,  G.C, , "Handbook  of  Stress  Intensity  Factors."  Institute  of 
Fracture  and  Solid  Mechanics,  Lehigh  University,  Bethlehem,  Pa,  1973, 

12,  Tada,  H.,  Paris,  P,  and  Irwin,  G.  "The  Stress  Analysis  of  Cracks  Hand- 
book," Del  Research  Corporation,  Hellertown,  Pa.,  1973, 

13,  Timoshanko,  S.P.,  and  Goodier,  J.N.,  "Theory  of  Elasticity,"  McGraw- 
Hill,  1951,  pp,  78-81. 

14,  "Fatigue  Behavior  of  Adhesively  Bonded  Joints,"  Quarterly  Progress 
Report  No,  2,  Air  Force  Contract  No,  F33615-76-C-5220,  prepared  by 
General  Dynamics,  Fort  Worth,  Texas,  January  1977. 


103 


APPENDICES 


The  computer  programs  for  the  mathematical  analysis  of 
adhesively  bonded  structures  are  outlined  here.  The  computer 
programs  for  the  solution  of  three  different  problems,  as  dis- 
cussed in  Section  2,  are  outlined  in  the  following  pages,  AIL 
the  computer  programs  are  written  in  FORTRAN  IV  language.  The 
programs  are  operational  on  IBM  and  CDC  computers.  The  essential 
features  of  the  programs  for  each  problem  are  outlined  in  the 
fo 1 lowing  appendices , 


APPENDIX  A 

A CRACKED  SHEET  WITH  AN  ADHESIVELY  BONDED  STRINGER 
AT  AN  ARBITRARY  LOCATION 
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INPUT  OF  DEFINITIONS 


The  computer  program  has  the  fol  lowing  inputs : 

I*  modulus  , Poisson's  ratio,  and  thickness  of  the  cracked  sheet 
(EP,  NUP,  HP) 

2.  modulus,  Poisson's  ratio,  width  and  thickness  of  the  stiffener 
(ES,  NUS,  DS,  TS) 

3.  Shear  modulus  and  thickness  of  the  adhesive  layer  {MUA,  HA) 

4*  length  of  the  debond  (B) 

5.  location  of  the  stringer,  i .e . , distance  from  centerline  of  sheet  (D) , 
The  distance  D is  always  positive,  or  zero 

6.  ha  If -crack  length  a (A) 

7.  number  of  increments  in  ha If -crack  length  (NAA) 

8.  increment  in  half -crack  length  (DA) 

DESCRIPTION  OF  OUTPUT 

The  outputs  of  the  program  are  the  shear  stress  in  the  adhesive,  and 
the  stress  intensity  factors  at  the  right"- (+A)  and  left-hand  (-A)  crack 
tips  (K/a/na), 
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LISTING  OF  PROGRAM 


C PROGRAM  DEVELOPED  BY  M.M.RATWANI  OF  NORTHROP  CORPORA! ION.AIRCRAFT  DIVISION 
C HAWTHORNE , CALI FORN1 A.  Tt  L6  PHONE  <2131  970-5285 

C PROGRAM  COMPUTES  STRESS  INTENSITY  FACTORS  AND  SHEAR  STRESSES  IN  A 
C CRACKED  SHEET  WITH  ADHESIVELY  BONDED  STRINGER 

C THE  PROGRAM  HAS  OPTION  TO  PUT  DtRUND  OR  NO  QE BUND  IN  THE  ADHESIVE 
C TFLESCUP1C  GRID  IS  USED  IN  THE  INTEGRATION  UF  KERNELS  WHICH  HAVE 
C LUGAR I ThIMIC  SINGULARITIES 

C I HE  PROBLEM  IS  REUUCEO  TO  SOLVING  AN  INTEGRAL  EQUAIION  IN  WHICH 
C SHEAR  STRESSES  ARE  IAKEN  AS  UNKNOWN  QUANTITIES 
C 

C INPUTS  ARE  - EP.NUP ,HP ,fc S,NUS, TS.OS , D , MU A ,H A , 8 , A,  DA,  NAA 
C N IS  THE  NUMBER  OF  COLLOCATION  POINTS 

C S.i.F  IS  THE  VALUE  UF  S1RESS  INTENSITY  FACTOR  IN  THE  OUTPUT 
C 

REAL  NUP iNUS  fMUPtMUSt  MU A 

COUPLE  X Z,Z2,ZR,ZB2,ZS,ZBS,ZZ,ZZB,ZZ2,ZZStZZB2,ZZ6S,THLl,TH12,TH13 
3, THlA, THAI , THA2,TH3L, TH32, TH33, TH3^,TH51 ,TH52»  THfal,TH62f  ZP, AZZ. BZZ 
3, AZZB, BZZBtZOM,ZPM 

OIMINSIUN  SB! 3faOOI  ,AB(60),ALI60l,CD(fa2l,ATl60I 
SO  FORMAT ( E 10. 3,AE 1 0. A) 
bl  F0RMA1 ( 2F 10. A ,bX,  I 21 
7 fa  FORMAT ( 2X, 8E 1 3. A ) 

BS  FORMA  T I X 1 1 7H  S.I.F  AT  + A = , E 1 5 . 7 , 5 X , l 7H  S.I.F  AT  - A = .E15.7I 

bfa  FORMAT I2X*5H  Y - ,8X,7H  TUY  = ,10X,2LH  HALF-CRACK-LENGTH  = ,F7.A> 

d 7 FORMAT! 2X  * 7H  MUA  = ,EtO.A,5X,bH  EP  = , E 1 0 . A * 5X , 6H  ES  =■  .E10.A) 

88  FURMAT12X, 19H  PL  AT  E - TH  I CKNESS  = , Ffa  . 3 . 5X  , l 8H  S T I FF  ENER  - AR  E A * ,F6. 

331 

d9  FORMA  1 1 2 X . 2 2 H ADHE S I VE- THI CKNESS  = tF6.3,5X,17H  DE  BOND-L  ENGT  H ^ ,F 

3 fa . J» 

VO  FORMAT  I X t b BH  DISTANCE  OF  CENTERLINE  OF  STRINGER  FROM  CENTER  OF  PLA 
3 TL  = , Ffa.  31 


C 

C 

C PL  A 1 1 PaRAMETFKS 


C F P = PLASTIC  MODULUS  OF  THE  PLATE 
C NUP  * PUIS StJN  RAIIU  OF  PLA1L  MATERIAL 
C HP  = THICKNESS  OF  THE  PLATE 
C 

RLAOlb.bO)  FPtNUP,HP 


C 

c 

c 


STIFFENEP.  PARAMETERS 


C ES  = ELASTIC  MODULUS  OF  THE  STRINGER 
C NUS  = POISSON  RATIU  uF  STRINGER  MATERIAL 
C TS  = THICKNESS  OF  STIFFENER 
C OS  = WIDTH  OF  THE  SIRINGER 

C D = DISTANCE  OF  THE  CtNTEKLINE  OF  THE  STRINGER  FROM  CENTER  OF  THE  PLATE 


C 
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P E AD  I 5 , 50)  ES.NUS,  TS,US  .0 
C AS  = AktA  Of  THE  STRINGER 
A S=D  S*  T S 
C 

C A*****************************  *******  ********************************* 

C 

C ADHESIVE  PARAMETERS 

C M UA  = SHEAR  MODULUS  Of  ADHESIVE 
C HA  = THICKNESS  OF  ADHESIVE  LAYER 
C 3 = LENGTH  OF  OE  BUND 

PEADI5.50)  MUA,HA,S 
C 

C #***********+********************************************************* 

C 

C CRACK  LENGTH  PARAMETERS 

C A = HALF-CRACK  LENGTH 
C DA  = INCREMENT  IN  HAlF-CRACK  LENGTH 
C NAA  = NUMBER  OF  INCREMENTS  IN  HAJ^F -CRACK  LENGTH 
R£A0|5,5l)  A.OA.NAA 
C 

c ******** ** ******* * **  ** ************ ******** ******** **************** 

P I =3. 1415V2 
N=  50 

NN3-  N*  2 
N 1 = 20 
N S=  35 
N2=N-N1 
N J = NS-N1 
NB=N-NS 
NR  1 = NQ* I 
NN  =N*  1 
NL=N*N 

UNP=  2 • • ( l.+NUP) 

UN  S*2 . * f l.+NUS) 

M UP=E  P/UNP 
M uS  = E S/UN5 
Pk= l 3.-NUP) /I l.+NUP) 

PK1=1.+PK 
PK  12=PKl/2. 

PK2=PK  L*PK 1/2. 

PKK=1.-PK 
A S S=A S*E  S 
PKR=1.-PK 

ZP 1 = 1./ t PI + 4. *MUP*HP*PK1) 

Z P 2= 1 . / I 4.*MUP) 

ZP3=1./I2.*P1*HP*PK1) 

YD=B 

YL-YD+O. A 
YL I* YD  + 1 . 

YL  2=  YD  + 4. 

DEL  1=1 YL-fOl/Nl 
DtL2=I YL1-YL) /N3 
D£l3*I YL2-YL1) /NB 
DEL4=IOELl+DEL2)/2. 

DtL5=(DEL2+DEL3)/2. 

CDI  1 ) = YD 

***** * * * ********  * * ** ** ** ** ** ** **** ***** **************************** *** 

DO  LOOPS  13V, 140,141  COMPUTE  Y CO-ORDINATES  OF  POINTS  AT  WHICH 
STRESSES  ARE  COMPUTED.  THESE  CO-ORDINATES  ARE  THE  MID-POINTS  OF 
THE  TELESCOPIC  GRID  USED  IN  SOLVING  INTEGRAL  EQUATIONS 
DO  13V  1=1, N1 
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139  GDI  l *1 ) =C0 I I 1*0EH 
DO  UO  1 = 1 , N 3 

1*0  COIN  1 + 1 ♦ I 1 =Cni N1 ♦ t ) +0612 

OH  1*1  1=1»NG1 

[ <t  i COINS*  l*-n*CU(NSM)*DEL3 
A SO  = DS  *MUA /HA 
wRITE(b,90)  D 
WRIIE16.87)  MUA.EP.ES 
HR  I lEtG.BB)  HP, AS 
WR1TEIG,89>  HA, YD 
WRITMG,86)  A 

NAf=NAA* 1 

C LOOP  ON  CRACK  LENGTH 

DO  800  NAP  = 1 ,NA  T 
A 2 = A *A 


THIS  LOUP  COMPUTES  THE  RIGHT  HAND  SIDE  OF  INTEGRAL  EQUAT IONS l EOUAT 
ION  20  UP  SECT  I UN  2.2.11  IN  NUMERICAL  ANALYSIS.  IT  IS  ALSO  LOOP 
FUR  EVALUATION  UF  KERNELSIEQUATION  20  OF  SECTION  2.2.11 
DU  69  1*1, N 
VI =C0 I 1*1) 

2 = CMPL  X ( D i Y l 1 
Z b = CON JG 1 Z ) 

Z2=Z*Z 
Z B2  = ZB*Z8 
Z S=CSQRTl Z2-A21 
ZBS=CSURT i Z B2-A2 1 

BS»PK1*A I MAG  I ZSI*PKK*Y1-2.*YI *REAL IZ/ZS) 

ABl 1 l=ZP2*ASD*B5 
ZDM=ZB-Z 

ZP-U*ZBS-ZB1  /(  A-ZB) 

ZPM=(-A*ZBS-ZB) /I-A-ZB1 

AL  < 1 ) =A l MAG ( ZDM* ( IZB/ZB$-1.1*ZP1/<A-ZB|-PK*ZP*(A*ZS-Z)/(A-ZH 
A HI  I = A l MAG  I ZDM*(  I ZtJ/ZBS-l  .1  *ZPM)  / t-A-ZBI-PK*ZPM+(-A*ZS-Z  1/  I-A-Z  1 1 
ALII) =“ AL I 1 t 
ATI  I ) =-  AT  ( I ) 


THIS  DO  LOOP  EVALUATES  KERNELS  WHICH  HAVE  LOGARITHMIC  SINGULARITI 

ES  AND  KERNELS  WHICH  ARE  NON-S I NGUL AR I EQUAU ON  20  OF  SECTION  2.2.1) 

DO  f,<V  J = l,N 

YJ  = CD I J+l 1 

ZZ-LMPLXIO, YJ1 

Z Z 8=C0U JG I Z Z I 

ZZ2=ZZ*ZZ 

ZZS=CSURT I ZZ2-A2) 

ZZl!2  = ZZH*ZZB 
/ZHS-LSUR I IZZD2-A2) 
lHU  = LLUGIZb*ZZB-AZ*ZBS*ZZaSl 
TH 1 2-C  LUG ( Zfi*ZZ-A2+ZBS*2ZSI 
TH1 J=CLUO( Z*ZZ-A2*  ZZS*ZS) 

MU**CLUGt  Z*/Zb-A2*  ZS*ZZUS) 

TH*2=0. 5*( ZBS-ZB-ZZS+ZZ1*! Z-ZB1/ I ( ZB-ZZ)*ZBSt 

AZ  Z = A2- Z Z 2 

AZZB=A2-ZZB2 

BZZ  = I ZZ-ZZB) /I ZB-ZZ) 

TH 32=0.  5*U.*ZS*ZZB S/A ZZBI  MZZB-ZZ1/  IZ-ZZ81 

TH  34=  0*  5 * ( l.*Z6S*Z2S/AZZ)*BZZ 

TH51=BZZ*( TH*2-0.5*I Z-ZBI * I 11/ 1 IS- 1 . 1 / Z 8S I 
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160 


161 


162 


30 

31 

a 

201 

202 

?oi 

205 

206 

207 

204 

30 

39 

69 


TH62-C L(JG<  Z-ZZb)*CL0G( ZB-ZZ) 

IF  IYI-YJ)  160,U1,16(J 
TH6l=CLOUt  Z-ZZ)*CLGG(ZB-ZZB) 
bzzb=i zzo-zz i n za-zzBi 

SP 1 =R  E A L ( TH61-TH62J 

TH31-0. 5*1  l.*ZS*ZZS/AZZ> *1 ZZ-ZZB1/J  Z-ZZ) 

TH33  = 0. 5*1  l. *ZbS*ZZ6S/AZZb) *BZ  ZB 

TH4  1 = 0. 5*1 ZBS-ZB-ZZBS*ZZB) •1Z-ZB1Z l(ZB-ZZ8l*ZBS) 

Th52=BZZB*l TH41-0. 5*1Z-ZB) *1 ZZB/ZZBS-l. 1 / Z BS » 

GO  TO  162 
EP 1=CD( J ) 

E Pk=CU ( J + l ( 

E P 2=CD I J *■  2 I 
FK  1 = 6 P2-EPR 
ER2=EP1-EPK 
E R 3=t  P 2-EP  L 
ER4=EK  3/2. 

YLPrt=2.*( ER 1*AL0GIABSIER1) )— ER2*ALQG(AHSI 6R2) J-ER3) 

YLPM=YLPM/ER4 
SP  !■=  YLPM-RE  AL(  TH62  > 

TH31=-0.5*Z*(ZZ-ZZG)/(ZS*ZS) 

T rl 33=-Q.  5*ZB*  I ZZtJ-ZZl  / ( ZBS*ZBS  ) 

Th4L=Q.5*I Z-ZB)*( ZB/ZBS-1. 1/ZBS 
THS2=-A2*IZZB-ZZ»*IZ-ZBJ/l4.*ZS*ZS*ZS*ZSI 
SP  2-RE  AL  ( PK1 2 *(  TH1  1-TH12*PK*(TH13-TH14)  ) *PK1*  ( TH41-TH42  1 *PK*(  TH3  l- 
3TH321+TH33-TH34+TH51-TH52) 

SP=Pk*SP1-SP2 
IF l Y i — Y J > 30,30,31 
SD=YI /ASS 
GO  III  32 
SD*YJ/ASS 
I L - 4 J- 1 1 *N* I 
IF  IJ-N11  201,202,203 
B 1 S= A $0  * ( SD- SP*ZP  1 ) *0ELl 

r.n  to  204 


U 1 Si A SO* I 5D-SP*ZP1 ( *0LL4 
GO  TO  204 
IF  r J— NS  J . 

B l S*A  SO  * ( S0- 
GU  TO  204 
B I S1 A 50* I SD- 
GO  TO  204 
b I 5= A SD* I SD- 


205,206 ,207 
D- SP*ZP 1 ) *0EL2 


3-$P*ZPl>  *Ofc  L 5 
D-SP*ZP1 ) *DEL3 


IFU-J)  38,3y,3b 
Sal IL 1=81 S 
GO  TO  69 
SB { I L 1 = l. + B1  S 
CONTINUE 


K S=  0 


The  SIMQ  SOLVES  FOR  UNKNOWN  SHEAR  STRESS  IN  THE  ADHESIVE 
CALL  SI MQI SB,A6,N, KS) 


795 


THIS  00  LOOP  PRINTS  OUT  SHEAR  STRESSES 
DO  795  I*1,N 

HR  1 TE  16,76)  CDI  1*1 ) , AB(I) 

SR-=0. 


C 


SR  2 = 0. 
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C THIS  DO  LOOP  CUMPUItS  NORMALIZED  STRESS  INTENSITY  FACTORS 

DU  SOS  1=1, N 
IF  { I — N 1 > 401  ,<,02,403 

401  Sk=SR»2. *A8( I ) *AL1  1 1 *DEL1/A 
SR2=SR2-2.*A6( I 1 * A T ( I 1 *DE  L L / A 
GO  TD  SOS 

402  SK»SKf2**AR( I )*AL( I 1 *DEL4/A 
SR2=SR2-2. *AB( 1 I *AT ( I > *DE  L4 / A 
GO  TO  SOS 

403  IF  I l-Nb)  406, 407,40b 

406  SK=SK*2.*AB( I 1 * AL  l 1 I *DFL2/A 
SR2=SK2-2. *ABI1 ) *AT( 1) *DEL2/A 
GO  Tl)  SOS 

407  SK=SK*2. *AGl l ) *AL II) *DEL52 A 
SR  2= SR 2-2. *AHl l ) *ATI I ) *DEL5/A 
GU  TO  SOS 

4 UB  SR  = SK+2.*AHU)  *AL(  1 1 *L)t  L3/A 

SR  2=  SR  2-2.  *ABIII  HH  l 1 *DFL3/ A 
SOS  LUNIINUb 

$Z  = SP*ZP3M. 

SZ2=5K2*ZP3+1. 

V-R1TE  (6,851  SZ.SZ2 
800  A = A *DA 
STOP 
PND 

SU6RQUT  t NE  S l RQ I A , B , N , KS 1 

D 1 ME  NS  1 UN  At  t)  ,6(11 

T(H_  = 0.0 

KS-0 

JJ=-N 

DO  6S  0=1, N 
JY=J+1 
JJ  = J J + N* 1 
6 I G A = 0 
1 T = J J- J 
DO  30  I = J , N 
1 J=I T , I 

IF (ABS(BIGA)— A6S(A( IJ) ) ) 20,30,30 
20  B l GA  =A  l [ J) 

I M A X—  I 

30  CONTINUE 

IFl AbStflIGAl— TOLJ  35,35,40 
35  KS=L 

Rfc  TURN 

40  I1=J+N*1J-21 

I T = IMAX-J 
DD  SO  K= J » N 
1 1= [ 1+N 
12=1  1 + 1 T 
SAVE=At 111 
A ( 1 1 1 = A ( 121 
At  I 2 1 = S A VE 

SO  At  lll=A( I ll/BIGA 
SA  Vt  = 0 ( I MAX) 

8(  I MA X > =0 l 01 
B(  J 1 =SA Vt  /U IGA 
lFU-N)  SS, 70,55 
55  I C $ = N*  t J- 1 1 

DO  65  1 X= J Y , N 
IXJ=IUS*1 X 
II =0-1 X 

00  oO  JX=JY,N 

1 X J X = N * ( JX-  1 ) + I X 
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J JX=  IX JX+ I T 

60  At  I XJX ) * A t 1 X JX ) - ( A ( 1 X J ) *At JJX)  ) 
6 b 8 ( i X ) = 0 ( IX)-lBtJ)*A(IXJ)) 

70  NV=N-1 
I T=N*N 

00  ao  j * 1 1 n y 

! AM  T- J 

1 8 *N- J 
IC»N 

DO  80  KM.J 

8 ( IB )=Bt I B l-Al 1 A) *B  t 1C  I 
I AMA-N 

eo  ioic-1 

At  TURN 
END 
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SAMPLE  INPUT 


5INQ.E  STRINGER 

rn^inFFH 

; 2L- 

*’**  iCBiAt  -*Q  NQP  -,o. 

mu  0HGN-  NO, 

S*  T t 

: 1 

1 ■ ™ - 1 1 1 l L_J £ - 1 

o u o a o g o o n t i 

■ _7  S * 5 * T fc  1 0 1 

1 2 2 * i i 7 1 *:  0 1 . J 

m l 2 i i ) t 

I * S 4 t 1 t 0 1 

a'VA  1 4 ’ * *01  7 2 t!l_T  ’ » T <0  I t T Lit  > A ) . P 1 1 . j , 1 . *_  , 9 * ^ ■ j— L_!_3— l-*-i 

10.  1 t:  + ll 

. 1 1 

. .0  fi  3 

1 

1(1,  1 F.  + fi 

. 3 3 

,12  5 

r j 

o 

0 . U E + 5 

. OOB 

0.0 

0 . 6 

. o..  o 

0 . 

- 

i ::::: j 

OUTPUT  FOR  SAMPLE  INPUT 


DISTANCE  of  cfntfrune  of  stringer  from  center  of  PLATE  = 0.100 

MUA  = 0. 6000E  *05  EP  1 0.1030E+08  ES  = 0.1030E+08 

PLATE-THICKNESS  = 0.063  ST  I FF ENE R- ARE  A = 0.250 

ADHESIVE-THICKNESS  = 0.008  DE BOND-LE NG TH  = 0.0 

y = Tiiy  = HALF-CRACK-LENGTH  = 0.6000 


0.2000E-01 

0. L325E+00 

0. 900QE- 0 1 

0.1 150E  +00 

0.6000E-01 

0. L092E  +00 

0.8000E-0L 

0.96L8E-01 

0.  1000E+00 

0.  89  75E-01 

0. 1200E+00 

0. 8928E-01 

0.  1 9 OOF  +00 

0. 7997E-01 

0. 1600E+00 

0.  7 5 1 9E  -0  L 

o.  laooE  + oo 

0.7117E-01 

0.20  00E  +00 

0.  fa  75  IE-  0 1 

0. 2200E  +00 

0.fa908F-01 

0.29Q0E  +00 

0.608  7E- 0 1 

0.2600E+00 

0.57S5E-01 

0.2800F+0Q 

0. 5501E-01 

0. 3000E+00 

0. 5233E-01 

0.3200E+00 

0. 9982E- 0 1 

0.3900E+0Q 

0.9798E-01 

0.3600E+00 

0. 9529E-01 

0.3800E+00 

0.932  IE-01 

0.9000F+00 

0. 3820E-01 

0.99006+00 

0.3276E-U1 

0.98006+00 

0. 295  IE- 0 1 

0.5200E+00 

0.27LfaE-01 

0.5600F+00 

0 . 246 1 F- 0 1 

0.6000E+00 

0.222faE-01 

0.fa900E+00 

0. 2013E-01 

0.6800F  +00 

C.  LB20E-01 

0. /200E+00 

0.  1696E-01 

0.7600E+00 

0.  1990E-01 

0.8000E+00 

0.  1351E-01 

0.8900F  +00 

0.  1226E-01 

0. 8BG0E+00 

0.  1 l 19E-U1 

0.9200F  +00 

0.  1010E-OI 

0. 06001+00 

0.89  77E-02 

0.  1000E+01 

0.6779E-02 
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o.  mot  *oi 
0.  UOOE+OL 
0.  1CjJJF*U1 
o. laoor »oi 
o.i*ooor  *oi 

0.2JOOE  * 0 i 

0.2600E  *01 
0.280QE  *01 
0.3000E  *01 
0 . 320QE  *0 1 
0 . 3AOOE  *0  1 
0.3600E  *01 
0. 3303F  *0  1 
0. A OOOE *01 
S. 1 .f  AT  * A 


O.<«Ot30L-O2 
0.  2890E-U2 
G.  1719002 
0.  UJ')6O02 

0 . 6l)‘t  5F  - 03 

0. 2897E-0  V 
-0.  W30E-U<* 
-0.2246E-03 
-0. 3286F-03 
-0.«t271E-03 
-0. 5323E-03 
-0.66731-03 
-0.8709E-03 
-0. 1316E-02 
= U.OO3SM0E  *00 


S.l.f  AT  - A 


0.  T077984E*00 
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APPENDIX  B 

CRACKED  SHEET  WITH  TWO  ADHESIVELY  BONDED  STRINGERS 
SYMMETRICALLY  LOCATED  ABOUT  THE  CENTERLINE 
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The  inputs  to  this  computer  program  are  the  same  as  for  the  single 
stringer  case  of  Appendix  A.  In  this  case,  distance  D is  the  distance  from 
the  centerline  of  each  stringer  to  the  centerline  of  the  crack.  The  shear 
stresses  in  the  adhesive  bonding  each  stringer  to  the  plate  are  the  same, 
also  the  stress  intensity  factors  at  both  crack  tips  are  equal. 

INPUT  OF  DEFINITIONS 

The  computer  program  has  the  following  inputs: 

1.  modulus,  Poisson's  ratio,  and  thickness  of  the  cracked  sheet 
(EP,  NUP,  HP) 

2.  modulus,  Poisson's  ratio,  width  and  thickness  of  each  stiffener 
(ES,  NUS , DS,  TS) 

3.  shear  modulus  and  thickness  of  the  adhesive  layer  (MUA,  HA) 

4.  length  of  the  debond  (B) 

5.  location  of  each  stringer,  i,e.,  distance  of  each  stringer  from 
centerline  of  sheet  (D).  The  distance  D is  always  positive,  or 
zero 

6.  half -crack  length  a (A) 

7.  number  of  increments  in  half-crack  length  (NAA) 

8.  increment  in  half  crack  length  (DA) 

DESCRIPTION  OF  OUTPUT 

The  outputs  of  the  program  are  the  shear  stress  in  the  adhesive, 
and  the  stress  intensity  factors  at  each  crack  tip  (K/a/Vra). 


118 


nrvonnno 


LISTING  OF  PROGRAM 


Q * **  * * *******  *************************************  **  ***  *****  *********  ** 

C PROGRAM  DEVELOPED  ti  Y M.M.RATWANl  OF  NORTHROP  CORPORA!  I ON , A I RCR  AFT  DIVISION 
C HAWTHORNE  , CAL IFORNI A.  TELEPHONE  12L3I  970-5235 

C PROGRAM  COMPUTES  STRESS  INTENSITY  FACTORS  AND  SHEAR  STRESSES  IN  A 
C CRACKED  SHEET  WITH  TWO  BONDED  STRINGERS  SYMMETRICALLY  LOCATED  ABOUT 
C CENTERL INE 

C THE  PROGRAM  HAS  OPTION  TO  PUT  DEBOND  OR  NO  OEBONO  IN  THE  ADHESIVE 
C TELESCOPIC  GRID  IS  USED  IN  THE  INTEGRATION  OF  KERNELS  WHICH  HAVE 
C LUGAR ITHIMIC  SINGULARITIES 

C THE  PROBLEM  IS  REDUCED  TO  SOLVING  AN  INTEGRAL  EOUATIUN  IN  WHICH 
C SHEAR  STRESSES  ARE  TAKEN  AS  UNKNOWN  QUANTITIES 
C 

C INPUTS  APE  - EP,NUPtHP,ES,NUS,TS,DSiD,MUA,HA,B(AtDA,NAA 
C N IS  THE  NUMBER  OF  COLLOCATION  POINTS 

C S.I.F  IS  THE  VALUE  OF  STRESS  INTENSITY  FACTUR  IN  THE  OUTOUT 
C 

C ****************  ************************************************ 

REAL  NUP .NUStMUP ,MUS,MUA 

COMPLEX  2, Z2.ZB, Z B2  , Z S » Z B$  , ZZ  » ZZ  B , 2 Z2  , ZZS  . Z ZB2  , ZZBS , T HI  1 , TH 12 , TH1 3 
3, TH  l4,TH4ltTH42. TH31( TH 32 , TH33 t TH34 t TH5 L ,TH52*  TH6 L , TH62»  ZP  , AZ  Z,  BZZ 
3, AZZB,BZZB,ZDM,ZPM 

3iZIiZBItZMI,ZBMI  t ZJtZBJ.ZMJ.Z&WJ,  ZPMM  , ZMPM , TMl  1 ,TML2,  TM13,TM14» 
3TM4l,TH42.TM31,TM32.TM33,TM34,TM5l.TH52*ZZMItZZMIB,ZZJM, ZZJBM 
3, TM6L.TH62 
3,  AL  f 

D I MENS  ION  SB l 360  0 I , ABI 6 0) » AL ( 60 \ , CD  I 6 2 I , AT ( 60) 

50  FURMAT (E 10. 3,4F 10.41 

51  FORMAT I2F 1 0 . 4 ( 5 X i 12 J 
76  FORMATI2X,BEI3.4I 

87  FOkMATl2X,7H  MOA  * , E 1 0 . 4 , 5X * 6H  £P  = ,ElO,4p5Xi6H  ES  * »E10.4J 

88  FORMATl 2X , 19H  P L A TF- I H I CK NC SS  * tF6.3*5X*l8H  STIFFENER-AREA  ■ tF6. 

331 

89  FORMA T I 2X t 22H  AUHt S I VE- TH I CK NE SS  1 tF6.3t5XtL7H  OEBUNO-L  ENGTH  * (F 
36.31 

86  FORMAT (2X (5H  Y 1 .8X(7H  TUY  = »10X*21H  HA LF-CR ACK-L ENGTH  * ,F7.4! 

85  FORMAT  IX»9H  S.I.F  - ,E15.7) 

93  F DR  M A T i X i 6 AH  DISTANCE  OF  CENTERLINE  OF  EACH  STRINGER  FRUM  CENTER  0 
3F  PLATE  = ,F6.3) 

*****+*************************#**********+*+************************+ 

plate  parameters 

EP=  ELASTIC  MOOULUS  OF  THE  PLATE 
NUP  = POISSON  RATIO  OF  PLATE  MATERIAL 
HP  = THICKNESS  OF  THE  PLATE 

REAO(5,50I  EP»NUP.hP 

£ **  **************  **********  *********************** 

C 

C STIFFENER  PARAMETERS 

C E5  = ELASTIC  MODULUS  OF  THE  STRINGER 
C NUS  - POISSON  RATIO  OF  STRINGER  MATERIAL 
C TS  = THICKNESS  OE  STIFFENER 
C DS  = WIDTH  OF  EACH  STRINGER 

C D = DISTANCE  OE  THE  CENTERLINE  OF  THE  STRINGER  FROM  CFNT E R OF  THE  PLATE 
C 

R E A D t 5t  50  I FS.NUS  , TS.DS  ,D 
C AS  = AREA  OE  EACH  STRINGER 
AS-DS*  T S 
C 
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c ********************************************** ************************ 

c 

C ADHESIVE  PARAMETERS 
C MUA  » SHEAR  MODULUS  OF  ADHESIVE 
C HA  - THICKNESS  OF  ADHESIVE  LAYER 
C B * LENGTH  OF  0E80ND 

READ! 5 i 50 ) MUA, HA,  B 

*****  ************  **********************************************  a****** 

CRACK  LENGTH  PARAMETERS 
A =>  HALF-CRACK  LENGTH 
DA  ■ INCREMENT  IN  HALF-CRACK  LENGTH 
NAA  * NUMBER  OF  INCREMENTS  IN  HALF-CRACK  LENGTH 
RE  ADI  5, 51 } A , DA | NAA 

*********************************** ****************  ******  *******  ****** 

Pi  = 3.  1AIS92 
N=  50 

NN  3=  N*2 
N L * 20 
NS*35 
N2-N-NL 
N3*NS— N 1 
NB*N-NS 
NB  1=NB* 1 
NN=N*  1 
NL*N*N 

UN  P = 2 » *(  l.*NUP) 

UN  S*  2 * M l.+NUS) 

MUP*  EP/UNP 
HUS-ES/UNS 
PK#( 3.-NUP)/! l.+NUP) 

PK  1=1  , + PK 
PK 12=PK 1 / 2, 

PK2  = PK  l*P  Kl/2. 

PKK  = 1 PK 
A SS=A  S*F  5 
PKR=  1 » — PK 

ZPl  = l ,/(P  I *A.*MUP*HP*PKl) 

ZP2=  1 ./ ( A,*MUP) 

ZP3=1  ./I2  ,*PI*HP*PKl» 

YD=  B 

YL=YD*0. A 
YL  1 = YD*  1 * 

YL2*=Yr*A. 

DEL  l = t YL-YDI/Nl 
DEL2-I Yl  l- YL ) /N3 
DE  L3=  I Yl  2-YL1I /NB 
DELA=(0ELWDEL2)  /2. 

DEL5=(DEL2*DEL3>/2. 

C D ( 1 } = YD 

*******  *************************************************************** 

DO  LOOPS  139, 1A0, LAI  COMPUTE  Y CO-ORDINATES  OF  POINTS  AT  WHICH 
STRFSSES  ARE  COMPLIED,  THESE  CO— ORD I NATES  ARE  THF  MID-POINTS  OF 
THE  TELESCOPIC  GRIO  USED  IN  SOLVING  INTEGRAL  EQUATIONS 
DO  139  I = 1 f N 1 
139  C0(  1*1 ) =CDl I ) *DE  L 1 

DO  LAO  1=1  ,N3 

l AO  CDINLH+t  l=CD(  NUI  ) *DEL2 
DO  LAI  t = 1 , NB 1 
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Ul  COINSM*  t)=CDl  NS*  I)*DE13 
A SD  = OS*MU A/HA 
W K I T F ( b , <5  0 > 0 
WR  I TE (6,871  MUA.EP.ES 
WRITE (6.88)  HP. AS 
WRITEI6.B91  HA, YD 
WRITE  (6,861  A 

C ******,♦*,,*,»*(■******»*+***************************#***************** 

NA  T=NAA+  1 

C LOOP  ON  CRACK  LENGTH 

DO  800  NAP  = 1 , NA  T 
A 2=  A* A 


THIS  LOOP  COMPUTES  THE  RIGHT  HAND  SIDE  OF  INTEGRAL  EOUAT IONS ( EQUAT 
ION  20  OF  SECTION  2.2.1)  IN  NUMERICAL  ANALYSIS.  IT  IS  ALSO  LOOP 
FOR  EVALUATION  OF  KERNE L S ( EQUA T ION  2D  OF  SECTION  2.2. L) 

DU  69  1=1, N 
Y l=CD( I *1 ) 

2 = CMPl  X( 0, Y! ) 
l B=CQNJG I Z ) 

Z 2=  Z*  Z 
Z B2  =Z  B*ZB 
l S=CSQRT( Z2-A2I 
Z8S*CONJG(2SI 

BS=PK 1*AIMAG(ZS >*PKK*Yl-2.*VI*REAL( Z/ZS) 

ABI  I t =Z  P2  * A SD*8  S 
21  = ZS-Z 
ZB  l=CONJG(Z ! ) 

ZM  l=-ZS*Z 
ZB M I =CDN JG I Z MI  I 
ZDM=Z-Z  Ft 
ZJ-Z/ZS-  l. 

I B J=CUN JG( Z J ) 

ZM J=ZJ 

ZBMJ=ZBJ 

A Z Z = A- Z 

AZ  Z 8 = A-Z  B 

8Z  Z = A + Z 

BZZB=A*ZB 

ZP  = (AtZI I /AZ  2 

Z PM3 | A ♦ Z fl I ) / AZ  Z B 

Z M PH= I At  ZM  I )/B ll 

ZPMM=(A*ZBMI IZBZZB 

AL I=PK1M-ZP*ZPM+ZMPM-ZPMH»*ZDH*< (ZP*ZJ)/AZZ*  ( ZMPM+ ZMJ  ) / 8ZZ  ♦ ( 
3ZPMM+ZBMJ  ) ZBZZBM  ZPM+ZBJI  / AZZB) 

AL  ( 1 > = A I MAGI AL I I 


C 

c 

c 

c 


THIS  DO  LOOP  EVALUATES  KERNELS  WHICH  HAVE  LOGARITHMIC  SINGULARITI 
ES  ANO  KERNELS  WHICH  ARE  NON-S I NGUL AP  I E OUAT I ON  20  OF  SECTION  2.2.1) 


DO  69  J*1.N 
Y J=  CD  l J* l ) 

Z Z = CMPL  XI D, YJ } 

ZZB=CONJGIZZ) 
l 12-11*12 
1 Z S = C SORT ( ZZ2-A2) 

ZZB2=ZZB*ZZB 

ZZB$=LONJGtZZS) 

TH1 l = CLOG l ZB*ZZ  B-A2*Z8S*ZZflS) 
T Hi 2 = CL0G ( ZH  + ZZ-A2+28S*ZZS) 
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TH13=CUNJG(TH1  1 ) 

THlA=C0NJGtTHl2) 

THA2=0.5*<  ZBS-ZB-ZZS+ZZ  1 * l Z-ZB m < ZB-ZZ3 *ZBS I 

AZZ=A2-ZZ2 

AZZB=A2-ZZB2 

BZZ*12Z-ZZB)/(ZB-ZZ) 

T H3  2 = 0 . 5* ( L.*ZS*ZZBS/AZZB)*(ZZB-ZZ)/(Z-ZZB» 

TH3A=C0NJGI TH32I 

TH5L=BZZ*(THA2-0.5*lZ-ZB)*(ZZ/ZZS-l.»/ZBS» 

TH62=CLUGI Z-ZZBl+CLOGIZB-ZZl 
TMl  l=CLOGI-ZB*Z2-A2-ZZS*ZBS> 

TM13=CLQG(-Z*ZZ-A2-2ZS*ZS) 

TM  14=C0N JG( TH1 L ) 

TML2=C0NJG(TM13) 

ZZMl=-2ZS*ZZ 
ZZMIB=C(JNJG(ZZHl  ) 

ZP=ZB*ZZB 
ZPM=  ZB  + ZZ 
ZMPM=Z*ZZ 
ZPMM=Z*2  ZB 
BZZB=ZZ-ZZB 

TM41  = 0.WDH*I  ZBl-ZZMIB)  /(  ZP*ZBS) 
TMA2-O.B*ZOH*(ZBI-ZZMI)/(ZPM*ZBS) 
TM3L=0.5*BZZB*ll.-ZS*ZZflS/AZZ8)/ZPMM 
TM32=-0.5*BZZ8*( 1 .-Z S*Z 1 S/AZZ ) /ZMPM 
TM3A=CONJG<TM32) 

TM33=CONJG( T M 3 1 > 

ZZJM-ZZ/ZZS-1. 

U JBM=C()NJG(ZZJM) 

TM51=BZ2B*tTM<,t-.5*ZDM*ZZJBM/ZBS)/ZP 
TM52=~BZZB*ITM42-.5*ZDM*ZZJM/2BS)/ZPM 
TM61-CL0GI ZPM) +CL0C< ZPMM) 

TM62-CL0GI ZMPM 1+CL0G1  ZP 1 
IF  (Yl-YJ)  160,161,163 

160  TH6  l=CtOG< l-Ll )+Cl UG( ZB-ZZBI 
B Z2  U= ( IZb—lZ ) / ( ZB-2ZB  1 
$P1=REAL ( TH61-TH62) 

$Pl=SPt*REAL(TM61-TM62) 

TH3  1 = 0.5*(U*ZS*ZZS/AZZ)*<ZZ-ZZB)/(Z'-ZZ1 
TH33=CONJG(TH31 ) 

TH5  1 = 0.5*IZBS-ZB-ZZ8S*ZZB>*U-ZB)/UZB-ZZBI*ZBS» 
TH52=BZ2P*(TH<»1“0.5*(  Z-ZB)*(ZZfl/2ZBS-l.  I/ZBS1 
GO  TO  162 

161  E P 1 =CD ( J I 

EPP=CD( 1*11 
£P2*C0( J+2 1 
FR  1 = EP2-EPP 
ER  2=  EP  l-E  PP 
E R 3 = E P2-F  P 1 
ER4=ER  3/2 . 

YL  PM  = 2.* (ER1*AL0G< ABSIERl) >-ER2*AL0Gt ABS1ER21I-ER31 
Yl PM=YLP  MZERA 
SP1=YLPM-REAL(TH621 
SP1  = SP  UREAL  (TM61-TM62  1 
TH3l=-0.5*Z*(ZZ-ZZB»/(ZS*ZS> 

TH33=-0.5*ZB*( ZZB-ZZ(/lZHS*ZBSi 
TH 41  = 0. 5*1 Z-28) *1 ZB/ZB  S-l.J/ZBS 
TH52=-A2*( ZZB-ZZ )*(Z-ZH) /<4.*ZS*ZS*ZS*ZS» 

162  SP2  = REAL(PK12*(  TH 1 1- T H 1 2* PK* I TH 1 3-TH l A 1 ) +PK  1*  ITH41-TM42 ) *PK*( THBl- 
ST  H3  2 > ♦ TH3  3-  TH34  ♦ T H5  l-  TH52  ) 
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SP2  = $P2+Rf  AL  l PK12  * I TM1 1-  TM  12*-PK*  I TM14-TM13U+PK  l*C  TH42-TM41  I*PK*l 
3TM31-TM32  )MM33-TM34*TH51-TM52  I 
SP=  PK*SP I -5P2 
IP  IYI-YJ)  30,30,31 

30  S0=  Y I /ASS 
GO  TO  32 

31  SD=YJ/ASS 

32  IL*(J-l)*NM 

IF  IJ-N1)  201,202,203 

201  B I S=  ASO* l SO-SP*ZP 1 ) *OEL 1 
GO  TO  204 

202  B 1 S=  A SU* ISD-SP*ZPl ) *DEL4 
GO  TO  204 

203  IF  IJ-NS)  205,206,207 

205  81  $=ASD*tSD-SP*ZPl)*DEL2 
GO  TO  204 

206  6lS=ASDMSD-SP*ZPl)*DEL5 
GO  TO  204 

207  81  S=AS0*(  SO-SP *Z P 1 ) *0E L 3 

204  IFU-Jl  38,39,  38 

38  SBIIU=blS 
GO  TO  69 

39  SBULI  = l.+BIS 

69  CONTINUE 

K S*  0 


THE  SIMO  SOLVES  FOR  UNKNOWN  SHEAR  STRESS  IN  THE  ADHESIVE 
CALL  SIHOISB.AB.N.KSI 


C 


c this  oo  loop  prints  out  shear  stresses 
DO  795  1*1, N 

795  WRITE  (6,76)  CDCIMl.AtMI) 

SP  = 0. 


c This  00  LOOP  computes  NORMALIZED  stress  intensity  factors 


DO  505  1=1, N 
IF  ( I-Nl 1 AO L , 402  » 403 
401  SR  = SR  F AB( I ) * AL ( I ) *OEL 1 /A 


GO  TO  505 

402  SP  3 SK*  ABI I ) * AL  ( I I *DEL4/A 

GO  TO  505 

403  IF  ( 1-N5  ) 406,407,408 
406  SR*SR+  AB( I ) * AL  l I I *DE  L2 / A 


GO  TO  505 

407  SR  3 SR ♦ ABI I )+Al tl ) *DE l 5/A 
GU  TO  505 

A 08  SR=SK+  ABtl I*ALI I ) *DEL3/A 
505  CONTINUE 

SZ=SR*ZP3*1. 


WRITE  (6,85)  SZ 
800  A = A + DA 
STOP 

SUBROUTINE  S I MQ l A , B ,N  ,K S ) 
DIMENSION  A(  l ) , B(  1 ) 

T OL  = 0.0 


K S=0 
JJ  = -N 

DO  65  J=t.N 
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JY= J* i 
JJ* JJ+N+1 
8 [GA»0 
I T* JJ-J 
00  30  I*J,N 
I J=  1 T 1 1 

t F C AB  S ( B 1 GA)-ABS( A ( IJ)I  | 2 0,30,30 
20  BlCAzAMJ) 
l MAX* I 

33  CONTINUE 

IF( A8$IBIGA)-T0L>  35,35.40 
35  KS=1 

RETURN 

40  IL-J*N*(J~2) 

I T= IMAX-J 

00  53  K3  J , N 

1 1=  I 1>N 

I 2* I 1 + 1 T 
5 A VE  = A ( r 1 ) 

A<  l 1 UAU  2) 

A(  I 2USAVE 

53  AM1)*AM  ll/BICA 
5 AVF  * B ( IMAX) 

B t f MAX)=B( J) 

BIJUSAVE/BIGA 
IFIJ-N)  55,70,55 
55  I0S=N*(J-1) 

00  65  I X= J Y ,N 
1 XJ  = I QS ♦ I X 

1 UJ-IX 

DO  60  J X=  J Y , N 
I XJX*NM  JX-1U  IX 
J JX* [ XJ  X+ I T 

60  AUXJXUAI  IXJXl-l  AC  1XJUAI  JJXt  I 
65  B(  I X I »B  t lX>-( Bt  J) *A(IXJ) > 

70  NY*  N—  1 
I T*N*N 

00  33  J*  1 , N Y 
IA- IT-J 

1 B*  N-  J 
I C*  N 

DO  80  K* 1 » J 

B(  I B ) =8  1 1 B I —A ( I At  *0 ( 1C) 

I A* IA-N 
80  ICMC-l 
RETURN 
END 


124 


SAMPLE  INPUT 


SAMPLE  OUTPUT 


distance  of  centerline  of  each  stringer  from  center  of  plate 

MUA  = O.6C00E*O5  EP  = 0.1030E*08  ES  * 0.1030E+08 

PLATE-THI  CKNFSS  * 0.063'  '“''STIFFENER-AREA  = 0.125 

ADHESIVE- THICKNESS  = 0.008  DE BOND-LENGTH  = 0.0 

0. 6000 


0.600 


= 

TUY  * HALF-CRACK-LENGTH 

0.2000E-01 

0.2S02E-01" 

0.4000E-01 

0.3398E-01 

0.6QQGE-OI 

0 . 3602E-01 

0.80006-01 

Oft  368  IE” 01 

0,  1000E  *00 

G.36966-QI 

0.12336400 

0 . 36  73 E - 0 1 

0-  14006  400 

0.  362 7 E- 0 I 

0.  16006  400 

0. 3566E-01 

3.  1 60  )E  400 

J.  3<t96E-  71 

0.2000E  400 

0.34206-01 

0.22006  400 

0.  3 3406” 0 L 

3 • 2 4 3 3 F 400 

). 32686“ 3 1 

0.26006 400 

0,31 76E-01 

0*28006 4 00 

0,3  0O4E-0 1 

3.  3 30  J )E  400 

3.3312E-01 

0.32006  *00 

0. 29 32 E- 01 

0.34006  *00 

0.28646-01 

3*  363)6 403 

3,  27786-  31 

0.  3800E  400 

0 , 2 700E- 0 1 

0.40006  *00 

0.25CKE-01 

3.44Q3E  400 

3. 2271 E- 01 

0.4  8006  #00 

0.21466-01 

G.5200E  *00 

0.2022E-01 

3. 560 36  *00 

3.  1 899 6- 01 

0.6000E  *00 

0.  L760E-01 

G.64G0E  *00 

0. 1 66  5F- 0 1 

3.6  8006 *00 

0, 1556E-01 

0.  72006  *00 

0.14526-01 

0.  7600E  *00 

0,13536-01 

0*  8 0006 *00 

0.  1 260  E- 0 1 

0.8 4 OOF  *00 

0,  1L7LE-01 

0.8 8006  *00 

0. 10076-01 

0.920QF +00 

0. 10026-01 

0.960GF  *00 

0. 90186-02 

0.  1 000  E +0i 

0. 76856-02 

0.  1 2 OOF  401 

£3.50366-02 
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0.  1 400E  401 

0. 3396E-02 

0. 1 600E  401 

0.220  IE-02 

0. 1 800E  401 

0.  1 162F-02 

0.20006  401 

0.7H35E-03 

i..?23.>r  di 

j.  iimr-m 

0.2400F  401 

0,  1 2 0 7 F - Cl  1 

0.26006 *01 

-0.6630E-04 

0.280JE  *01 

-3.2013E-03 

0.30006*01 

-0. 3076E-03 

0.32006  *01 

-0 . 40  78E- 03 

0.3430E  *01 

-3.62146-03 

0.3600F  *01 

-0.6  779E-03 

0.3800E  *01 

-0 • 9 368  E-03 

3.400  36  *01 

-3.1489E-02 

S.I.F  * C.5448927F*00 
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TWO  ADHESIVELY 

APPENDIX  C 

BONDED  LAYERS  WITH  A CRACK  IN  ONE  LAYER 

(SMALL  CRACK  LENGTHS) 
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A generalized  computer  program  for  this  problem  could  not  be  written 
due  to  the  fact  that  load-shedding  to  the  sound  layer  is  a function  of 
crack  length,  as  well  as  adhesive  and  adherend  properties.  The  conver- 
gence of  the  solution  is  strongly  dependent  on  crack  length.  Hence,  two 
programs  are  written,  one  for  a small  crack  length  (half-crack  length 
less  than,  or  equal  to  0,4  inch)  and  the  other  for  a long  crack  Length 
(half-crack  length  larger  than  0,5  inch,  and  smaller  than,  or  equal  to 
one  inch).  Similar  programs  can  be  easily  written  for  crack  lengths  lar- 
ger than  one  inch  by  taking  more  collocation  points  in  the  numerical 
analysis , 

The  computer  programs  are  written  for  the  case  of  no  debond  in  the 
adhesive,  or  the  debond  approximated  as  elliptical  in  shape.  The  programs 
have  the  option  to  input  an  elliptical  debond  with  a minor- to-major  axis 
ratio  of  0,0  (no  debond),  0,1,  0.2,  and  0,3,  Other  debond  shapes  and 
sizes  can  easily  be  easily  incorporated  by  changing  the  region  and  shape 
of  the  surface  over  which  integration  of  the  kernels  is  carried  out, 

INPUT  DEFINITION 

The  computer  programs  have  the  following  inputs: 

1,  modulus,  Poissons  ratio,  and  thickness  of  the  cracked  layer 
(El,  MU 1 , HI) 

2,  modulus,  Poisson's  ratio,  and  thickness  of  the  sound  layer 
(E2 , NU2 , H2) 

3,  shear  modulus  and  thickness  of  the  adhesive  layer  (Ml'A,  HA) 

4,  ratio  of  the  minor- to-major  axis  of  the  elliptical  debond  (DEB) 

5,  half-crack  length  a (a,  less  than  or  equal  to  0,4  inch) 

6,  number  of  increments  in  the  half-crack  length  (NAA) 

7,  increments  in  the  half-crack  Length  (DA) 
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DESCRIPTION  OF  THE  OUTPUT 


The  outputs 
Locations  in  the 


of  the 
bonded 


programs  are  shear  stress  (t  , 
region,  and  the  stress  intensi 


T ) at  various 

y 

ty  factor  (K/a/na), 


LISTING 


C 

£t**4**it*#e**#¥:t*9#***#*4t*************************  *********##*****#*♦ 

C PKOgR AM  DLVLLCjPE  D by  M.M. RATWAN I OF  NORTHROP  CURPOR  A T ION  , A I RCRA  FT  DIVISION 
C HAWTHORNE, CALIFORNIA.  T EL  tPHUNE  !2I3)  970-S285 

L PROGRAM  COMPUTES  STRESS  INTENSITY  FALTOkS  AND  SHtAR  STRESSES  IN  TWO 
C PLY  ADHESIVELY  BUNDED  STROClJkt  WITH  A CRACK  IN  ONE  LAYER.  THE 
C PROGRAM  HaS  OPTION  TO  POT  NO  Dl. BONO  IN  THE  ADHESIVE  UK  ELLIPTICAL 
L OEdLND  WITH  MINOR  TO  MAJOR  AXIS  RATIO  OF  0.1, 0.2  UR  0.3 
C THIS  PRUGRAM  GlVtS  CGNVtRGINb  RESULTS  UPTD  A = G.A  INCH  AND  RATIO  OF 
C AOHtSIVE  IliICKNLSS  TU  SHEAR  MUDULJS  - 0.1E-6 

C TELESCOPIC  GRID  IS  USED  IN  TriE  INTEGRATION  OF  KERNELS  WHICH  HAVE 


L LGoARITHIMIC  SINGULARITIES 

C THE  S 1 2 E OF  TELESCOPIC  GRID  DEPENDS  ON  CRACK  LENGTH 

c the  PKtmLfM  is  kiDoctu  tu  solving  a set  uf  integral  equations  in  which 
c shear  stklssf_s  tux  and  tuy  are  taken  as  unknowns 

C SOLUTION  UP  UNKNOWN  SHtiAR  S1RESSES  TUX  AND  1UY  1$  ACCOMPLISHED  BY 

c solving  a set  of  simultaneous  euuatiuns 

C X,Y  ARt  The  Cu— GROIN AlES, FkoM  THE  CE NTRL-L I NE , AT  WHICH  SHEAR  STRESSES 
C Akc  COMPUTtD 

C S.I.F  IS  THE  STRESS  INTENSITY  FACTOR  IN  THE  OUTPUT 

L NX  1 , NX2 ,NY  1 ARE  THE  NU.ibER  OF  COLLOCATION  POINTS  IN 

L ThL  TELESCOPIC  GRID  FUR  NUMERICAL  ANALYSIS 
C INPUTS  ARt  - El ,NU1 ,Hl , E 2 ,NU2 , HZ , MU A , HA , DtB , A ,UA,NAA 

COMPLfcX  2,Z2,Z.i,Zb2,ZdSt2S,ZZZ  ■ Z I , Z I d , Z JB  , PHZ 1 , PHZ2 , PH2 3 , PH l , PH2 , 
3F  1 tF2,lH21 1 Th2^,ZZ,  THS,ZZti,2Z2,  ZZ2A,ZZS,ZZb2  , ZZ2b»ZZBS»ZPH,  ZPHD, 
3ZZ2ZB,ZZI,ZZIu,ZZJrt,ZZJ,Thl I , T HI  2 , T Hi  3 , TH 14 , TH31 , THA2,Z21 , TH1.TH2, 
3TH32,TmA1,Z1Z,Z2Z,2L1,Zc2,Z3ZB,29ZR,ZPL,ZPLB,PhZA 
S.ZXY.ZAZ.ZIZd.ZiZb.ZLS.ZLA 

3,  Z J,ZlM,ZIbM,  ZJM,ZJuM,PH3,PHA  ,ZPHP,ZPHHP,  ZM,ZDM,Z3[iM,ZdMM,  (Mil, 

3TMlc,  IMli,  I M 1 A , Z P HP  ii  , Z P hP  P , Z b i tP  , Z LtP  P , lM31,TM32.lM33,TMjA,TH33, 
3Th3A,ZZIM,2Z 1RM,TmA3, 1 HA9 , T Hb 1 , Th62 ,TH53,THSA,ZZ JliM,ZZJM 
REAL  NU  1 1 NU2  » MUA 

DIMENSION  ZXY ( lUAl , FI 20B I ,Ctol A326A)  ,F 1(  IGA ) ,F2( 1CA) 

AS  FORMAT  l 2X , 1 6 , At  I 3 . A ) 

50  FORMA  I 1 1-10.3, 3F1U.A  1 

51  FORMAT t 2F1G. a, SX, IZ ) 


76 

87 

83 

69 


F URMA T i 
FORMA H 
3U  = ,F7 
FORMA  I ( 
3 NU  = * 
FORMAT ( 


2x,bLiJ.A) 

X , 2 r»H  SUUND 


LAYER-THICKNESS 


»F7. 


• A 1 

X , 27H  CRACKED  LAYER-THICKNESS  * ,F 
F 7 . A ) 

X , 2tH  ADHtSIVE  - SHEAR  MODULUS  - , 


A.2X.SH  E = ,E9.3, 3X.6H  N 
7. A, 2X»5H  E * ,E9.3,3X,bH 
E10-A. 5X, 13H  THICKNESS  * 


3 , F 7 . A ) 

100  FORMAT ( 2X , 1 1E9 .3) 
es  FuRMAT  (X.VH  S.I.F.  =■  f E15.7) 

£.6  FORMAT  ( 2X , SH  X - *dX,SH  Y - ,8X,7H  TUX 


,6X,7H  TUY 


, 10X.5H  A * 
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91 


3,F7.A)  .. 

FORMAT! X,33h  NO  DtBOND  IN  THE  ADHESIVE-DEB  - ,F5.3I 

FORMAT! X,5dH  ELLIPTICAL  DL BUND  IN  AUH t S 1 V L-MI NOR  TU 
310=  iFA.l) 


MAJOR  AXIS  RAT 
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c**************************** ****** ************************************* 

C CRACKPD  LAYER  PARAMETERS 

C El  * YUUNbS  MODULUS  OF  CRACKED  LAYER 
C NU 1 = POISSONS  KA 1 10  OF  CKACKCO  LAYER 

c hi  * thickness  jf  cracked  layer 

READ  IS, SO)  E 1 iNO It  Hi 
C 

^** **********************  *********************************************** 

C SOUND  LAYER  P Aft  AM  L T E R b 

C 62  = YUUp.GS  MODULUS  uf  sound  LAYER 
C NU2  * P01SSUNS  RATIO  uF  SUUNU  LAYtR 
C H2  s THICKNESS  OF  SOUND  LAYER 
READ  (5, SC)  fc  2 i Nu  2 t H2 
C 

c ************  **************  ********************************************* 

c aohesiye  parameters 

c ha  * thickness  uf  adhesive 
C MUA  * SHEAR  MODULUS  OF  AUH6S1VE 

C Dt. iJ  = RATIO  CF  MINOR  TD  MAJOR  AXIS  OF  ELLIPTICAL  SHAPE  OF  DEHONO 
C D to  * 0. | 0. It  C.2  UR  0.3 
READ  (5,50)  MUA t ti A • DE b 
C 

(,***********♦****.***************** ****************  ********************* 

C CRACK  LENGTH  PARAML1LRS 

C A * HALF  CRACK  LENGTH 
C OA  * INCREMENT  IN  HALF  CRACK  LENGTH 
C NAA  = NUMbER  UF  iNLKLMENTS  IN  HALF  CRACK  LENGTH 
READ  (5*51)  A i UA ( NA  A 
C 

P 1*3. 1A1592 
PI2*2.*P1 
P I22=P 1/2. 

PK1*< 3.-NU1I/I 1 .*NU 1 ) 

PK2*  ( 3 . — NU2 ) / ( 1 • +NU2 ) 

PPK1*P12* ( l.*PK!> 

PPK2=P 1 2* 1 1 ,*PK2) 

PKS*PK 1 *PK1 
PKU*(PM-l.)/2. 

PKlDsPKl/2. 

Gl=Ul/(2.*il.+NUl)  ) 

G2*t </ ( 2 .* ( 1 .+NU2) ) 

GA  = A . * 1 

012*2. *Gl 
G22=2 . +G2 

pkz*pki-i . 

PKP=PK1* i. 

H2R*H2*PPK2*G22 
HRR*H  1 *P  PK 1 
h1RshRK*G12 
hg=ha/mua 
PIN*!. -PI/2 . 

WRITE  lotSU)  H!,E1,NU1 
WRITE  ( 6, o7 ) H2.E2.NU2 
WRITe  (6.891  MUA, HA 

C THE  SETTING  OP  <X,Y)  CU-ORDINATES  OF  TELESCOPIC  GRIO  STARTS. 

NA3*NAA*L 

C THE  FOLLOWING  DO  LOUP  VARIES  THE  CRACK  LENGTH 
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QU  370  NK  J = l t NAD 
NX  1=3 
NX2  = Z 
NX3=9 
NX4  = 1 
NXS=2 
NY1  = 6 
NY2  = 9 
NY3  = 6 
NYA  = 2 
NY5  = 3 
A2  =A* A 

UY 1= A/ i 2 *NX 1 ) 

DY2  = 2 • *UY 1 
DY3=2.  *OY2 
DYA=2.*0Y3 

IF  (Util-0. 1)  570,571,169 

169  I F { DEB— G .2 } 170,170,171 

170  NRZ«NX 1-9 
WRITE  (b,91  1 DEB 

NR3=  1 

NPS=NX 1* 1— NR2 

N Y L = NY 1 +NK  3 

NXY=.4X  1*NY 1*NPS 

EP=DYl*DY2 

DU  179  1 KL  =NR2  »NX1 

I KP  = 1 KL-3 

EXM l*  IKL~1>  *DY1 

179  7XY(lKPl=CHPtX(EX,EPI 
GU  To  1B3 

171  NR2=NX1-1 
WRITE  (6,91)  DE3 

NR3=NX  1-9 

NPS=2*(NXl+l )— NR2-NR3 

NX  Y =NX  1*  NY  1+NPS 

EP=DY1+DY2 

DO  ISO  IKL=NR2.NX1 

lKP=IXL-6 

EX=12MKL-U*0YI 

180  ZXY(  Jk.P)=CMPLX(EX,EP> 
EP=tP+0Y2 

DO  181  IK.S-NK3.NX1 
IPP=1KP+ 1KS-3 
EX=1  c*  1K.S-1M0Y1 

181  ZX  Y ( I PP ) =CHPLX I EX , EP ( 

GO  TO  163 

670  tP=DY 1 
NR  l = 1 
NR  2=  1 
NR  3-  1 

WRITE  ( fc , 9G 1 DEB 
GO  TJ  67 2 
971  EP=DY1 

NR  1 =NX 1-1 
NR2=  1 
NR  3 = 1 

WRITE  16,91)  DEB 
672  NP  S = NX 1 ♦ 1— NR  1 
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non 


NYI=NY1+NK2*NR3 
NX Y=NX  1*NY1+NP$ 

DO  873  I =NR 1 1 NX  1 
EX=( 2*1-1 I *D  Y 1 

e?3  ZXY ( 1 + 1-Nkl I =CMPCX( EX.EP1 
133  DU  573  J=l,NYi 

NP= ( J— 1 ) *NX  1 + NP5 
EY  = E P* J*DY2 
DO  573  I = 1 * NX  1 
NQ=NP* I 

l:X  = l 2*  1 — 1 ) *DYl 
ZXY (NO  )=CWPLX( EX,E Y ) 

573  CONTINUE 

cP  = EX  + LiYl+DY2 
DO  574  J=1,NY2 
NS*NXY*(J-1)*NX2 
£Y=<  2*J-1 I *DY2 
00  574  1=1, NX2 
nm=n$+i 

HX*EP* l 1-1 )*DY3 

574  ZXYtNM  ) = CMPLX|EX»EY) 

JXY=NXY+NY2*NX2 

1XY=JXY+NY3*NX3 

EP  = EYOY3 
DO  575  J=I,NY3 
EY=EP* 1J-1 )*UY3 
NP=JXY+( J-1)*NX3 
DU  575  1=1, NX3 
fcX=(2*I-l)*DY2 
N0=NP-*-I 

57  5 2XY  (ftfc  ) = CMPLX(EX,EY1 
EP=EX+  DY  2 + OY  3 
DO  57b  J = 1 , N Y 4 
E Y x l 2 * J - 1 )*DY3 
NP= 1 XY ♦ ( J— 1 ) *NX« 

DU  57b  1=1, NX4 
Nw=NP+ 1 

EX=EP* ( I-i »*UY4 
57b  ZXY1NQ  ) = C4PLX  lEX.EYl 
KXY=1X  r«NX4*NY4 
EP=A]MAG(ZXY 1 1XY  1 )4DY2*DY3 
KN=KXY+NX5*NY5 
DO  5 77  J=1  ,NY5 
EY  = EPMJ-1J*DY4 
NP=KXY ♦ ( J— 1) *NX5 
DO  577  I = L ,NX5 
NQ=NP+  I 

EX=I2*I-1  )*DY3 
577  ZXY(NQ  )=CMPLX( tX,EY) 

C THE  SETTING  CjP  TELLSCQP1C  GRID  LND5 
NXN2=2*KN 
KNN=NXN2*KN 
KNM=KNN*iiN 

EVACUATION  UF  KERNELS  K II , K1 2 , K2 1 ,K22  OF  INTEGRAL  EQUATIONS  45  OF 
SECTION  2.2.3  STARTS 

FULLUWING  LUOP  IS  UUTER  LOOP  FOR  EVALUATION  OF  KERNELS  VARIES  Z 
DU  b*  1=1, KN 
IP  = I 
I Q=KN+ I 
Z=ZXYl I ( 

Z 2 = Z * 2 
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Zd*CO  JJG(  Z ) 

ZS>2  = Zb*ZB 
ZS-CSjKT l Z2-A2 ) 

ZfiS*CUNjG<ZS) 

III  = Z-Zli 
Z21=Z/ZS 
zaz=Z3/zbs 
c 

C F ( I P ) IS  RIGHT  HAND  SIDE  OF  FIRST  OF  INTEGRAL  EQUATIONS  *5  OF 

C SECTION  2. 2. 3 

L Fdul  IS  RIGHT  HAND  side  of  second  of  integral  equations  as  of 

C SECTION  2. a. 3 

c 

FI  IP) -PKZ*REAL( ZXYt  I ) J+PKI*REALt  ZSJ-REALI  ZBS ) +2 .* A IMA&I  ZXYt  1 ) ) * 
3AIMAGI ZAZ I J/GA 

F ( IQ  ) =(  -PKZ+AlMAGIZXY  ( I ) ) +PKl*A  I MAG  ( Z S ) -A  IMAG  I Z US  )-2.  *A  IMAG  I Z X Y I 1 ) 
3)*KEALI  ZAZt  l/G*. 

ZI=ZS-Z 
Z IH=CUNJGiZ l 
ZIZ=A-Z 
Z2Z=A-Zb 

ZLl=A*Z 
ZL2=A*Ztt 
Z J=Z/ZS-I. 

ZJB=oONJG(ZJ> 

Z IH=-Z  S + 7 
Z 1 :jM- 4 jJblt  U.  > 

ZJ‘isi  J 
Z J 6 'I  - L Jl 
pm/ 1=  u ♦; j )/i\i 
PHi  ?- u«z  1 >/Z_Z 

PrlZ3=IAiZi  D/ZLI 
fnZ**a  l A + l 1 b-'i)  /Z  Li 
PHl=( Pri2* «ZJ J/Z 1Z 
PH2  = f PHZ2 ♦/ Jd  1/Z2 Z 
PH3=  ( Pti/i  + Z J‘A  ) /ZL1 

Z°L=-PH2  1 ■*rJM-M’HZi  + ZZZi*{JHl*PHZ3-PKl  wPHZAtZZZ*Pri3 
ZPLj=-22Z*  PHL’-PHZ  2*  PKi*PnZ  1-7  1 1 * PFw  + P H7  A-W*l* PH 7 3 
Ull  ) = ZPL •* Z P Lu 
Fill) -ZPL-ZPLb 

INhcR  L J J P KJP  LVAlUATIUi  CjF  KERNELS  VARIES  ZO 
00  6*-  U = 1,K.N 
NP  = ( I L-l ) > *XNi * 1 
NQ=K,,*NP 
NS=6  M+\P 
NH  = NN  1 + NP 
7 Z =7  X V ( 1 1 ) 

IF  Ill-iXVI  669,667,670 
t-69  LIX1=0Y2 
GU  1 0 671 

670  IF  Ili-1XY|  672 , 672  « 67  3 
67?  0X1=  iiY  3 

GO  T 0 u7  1 
673  DXl  = .'YA 

GO  TO  671 
67  1 DYYsJXl 

DXY  = DX  1*DYY 

10  = 1 I 
IS=kN+  1 i 
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ZZj  * t L)'<  JG  ( 7 Z ) 

Hi  = ll*Ll 
ZZ2A  = Z2:-A2 

zzs  *ci*.an iitk) 

22B2  = 2 Zt>*Zi.  u 
Z 2 2 b = 2 2 L i_ -A  2 
ZZBS-C  J'4JGI  ZZS  ) 

ZPH  S l-Ll 

ZPHB=CCJNJGtZPH) 

22  Z Z 8 = ZZ-ZZfl 
ZZI  * ZZS-Z2 
ZZIB*CONJG(ZZI» 

ZZJ=ZZ/ZZ$-1. 

2Z  JB-CUF4  JG  ( ZZ  J ) 

Z PHP=Z ♦22 
2PHtsP=CQNJb(ZPHP  ) 

ZPhPb=Z-ZZa 
ZPhPP  ®Z  *228 
ZdHP-CLiNJGIZPHPB) 

ZbPP=GU<N  JG  I Z PHPP  ) 

ZM=Z*ZZ 

ZBMsCONJGUM* 

Z BbH  = Z*Z Z G 
Z8MM=COFiJG(Z3BM) 

7 Hi 3aLLCGIZM-A2+ZZS*ZSt 
TH14=CLGG( ZBBM-A2*ZZBS*ZS) 

TH1 1=CDNJG(Th13) 

7'Hl2-COlvJbtTMlM 
T-M12*CLC)G(  -Zb8M-A2*2B S* ZZ S > 

T M 1 3=Cl_ uu C -Z8-A2-ZZS*ZS) 

TM11  = Cl»UJGITM13) 

TMU=CGNJG(TM12) 

PH3=CLOG( 2 PHP ) 

pha=coi;jg(  PH3  t 

PHZ  1 *C  LCb  ( Z PHPtl ) 

PHZ2BCUriJGiPHZl  1 
PhZ  3 = CL  uG ( ZHHPP  ) 

PHZ A=C  O.'l  JG  ( PhZ3  ) 

Z L 1=C • 5 * ( l.-ZS*ZZBS/ZZ2B)*ZZZZ8 
TH31*-ZL1/ZPhP6 

TH31»0.5*{ l.+ZS*ZZBS/ZZ2Bl*ZZZZB/ZPHPP 

TH33®COFvJg(  TH31 ) 
n33»CGNJut  fM31> 

7832-  0.fj*{  1.  ♦Z8S*ZZbS/£ZZB»*ZZZZB/ZPHBP 
7M3**aCCiN  JG  17  M32  | 

2ll  = 0, S*ZZ2 
ZZI»Zi.S-ZZ 
ZZIBSCC.JJGIZZI) 

Z21M=-tZS«ZZ 
ZZIbMsCarjJUl  ZZ  IM  ) 

7H*2«ZLl* I Zlb-ZZ lbH)/(ZPHUP*ZBS ) 
7HA3aZCl*(Z IB— ZZI )/CZbS*ZBHP> 

THAA=ZU*(  2Ib-Z21M)/(2bS*ZBPP1 

ZLl*ZLl/ZbS 

2ZJ  = < ZZ/ZZS-l.  ) 

Z Z Jb  =CCjN  JG  ( Z Z J ) 

ZZ  JB8  = ZZ J J 
ZZ JM»ZZJ 

TH51*HH4  3-ZL1*ZZ  J)*ZZZZ(J/ZBHP 

TH5  2 = -( THA4-Z  L 1*ZZJM I *ZZZZB/ZBPP 
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TH5A=(THa2-ZL1+ZZJBM)*ZZZZB/ZPHBP 

TH22=CLDG( Z + ZS) 

TM1 =CLOG ( Z -Z  S ) 

Th21=PK11MPK1+Th2 2— C0NJG1 TH22 )— PK1*TH1+C  QNJG( TH  L ) ) 

IF  (I. £0.11)  GO  TO  77 
PH  1 =C LOG  t Z PH ) 

Pm 2-CON JG (PHI) 

TH2-ZPh/ZPHb 

TH32=-0.5M l.-Z&$*ZZBS/ZZ2B)*ZZZZB/ZPHB 
TH3*=C0NJG(TH32) 

TH53=-ZZZZBMTh9l-ZLl*ZZJB)/ZPHH 
THAl=ZLl+(  Zlb-ZZlB)/ZPHil 
GO  TO  78 

77  THA1=0.  b*UZ  + (ZB/2BS-1.  )/ZBS 

TH3^-0.5*ZB+ZZZZB/tZBG+ZZBS) 

TH3A=C0NJG(TH32  ) 

TH53  = 0.  2b*A2*ZZZZb*7ZZ/(  < Z B2-A2 ) * ( ZB2 -A2 ) ) 

0X2= ( OX  1 »0X1  )/9. 

DX 22  = 2 .*0X2 

DXPI=DX22+PI 

TH2=CMPLAIDXP1,0.) 

TH2=Tn2/0\Y 

PHl=DX22*l A LOG ( 0X22 ) +P 122—3. ) 

PHi=t‘rtl/OXY 
PH2  =PH  1 

79  ZL l=-t Pul  * PH2  ) +PH3+PHA 

ZL2=ZPHPb/28HP-ZPHPP/ZBPP 
ZL3=-tPMZl*PHZ2)+PMZ3*PHZA 
ZL*=TH2-ZPHP/ZPHbP 
Z1ZB=PK2*2L 1+ZL2 
Z2ZB=PK2*ZL3+ZLA 

Z3Z&=PKl*iLl*ZC^  + PMD*(  TH13+THIL-TH13-TH1 1 ) *0 .5 * t -TH  I2-PKS* TH  14* 
3TM12  + PXS*TM1A)+TH51— Ttt52*PKl*(TH3A— TM  39+THA 1— TH92  (-TH33+T  M33-TH<»3 
3*  T MAA  ♦ T tu  1 

ZAZB=PK1*ZL3*ZLA+  TH53— T H5A*PK10*t  TH1A+Th12—TM1A—TM12I+0.5*(-TH11- 
3PKS*Thl3*TMH  + PKS*T«l3)  *PK1*<  TH31-TH3 1 + THA3-THAA )-TH32*TM32-THAI  + 
3 THA2+TH21 

R1  1 = RF:  A L < ZIZ3)*R£AL  ( Z2Zi>  ) 

R12=— AIVAGI ZIZb) ♦ AlMAGlZi’ZB ) 

R2L=AIMAG(ZlZb}*AlMA&tZ2ZB> 

R22=RLALl ZIZb )-R£AL (Z2ZB) 

H 1 1 =RE A L I Z3Zo)+RtAL(ZAZb) 

Hl2=-AIMAG( Z3Zb)+AlMAG( 2AZB  ) 

H21  = AlMAGt23ZB  ] ♦AIMAGUAZD) 

H22  = kE  A L { Z3ZB)  —RE AC ( Z *tZ  3) 

51  l = (hl  1/H1R+R11/M2R)  +IjXY 
S12= ( H12/H1R+R12/H2R ) *DX  Y 
S21=t n21/H1k+R21/H£R)  *OXY 

52  2=IH2  2/HIR*R22/H2R)*L)XY 
IF  (l.LQ.ll)  GU  TO  789 

GO  TO  790 

789  SI  1 = HG  + Sil 
S22  = HG  + S22 

790  Cu(NP)=SIl 
C!i(N0)=S2I 
C B ( NS ) =S 1 2 
CB  t Nil  1 =S22 

69  CONTINUE 

E VACUA (ION  UF  KERNELS  ENDS 

S 1 Mo  SOLVES  FUR  UNKNOWN  SHEAR  STRESSES  TUX  AND  TUY 
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CALL  SIMg(CB,F,NXN2,KS) 

WRITE  (6,fl6)  A 
DO  889  1=1. XN 

t6S  WRITE  (6,76)  ZXY( II •FIII«FIKN*XI 
SUM=Q. 

C FOLLOWING  LUUP  COMPUTES  STRESS  INTENSITY  FACTORS  FROM  EQUATION 

C Ail  OF  SELT  1 UN  2.2.3 
00  269  1=1, KN 
IF  U-NXYI  268,266,271 

268  0X1=UV2 
GO  TU  270 

271  IF  (1-IXY)  272,272,273 

272  OX  1 = UY  3 
GO  TO  270 

*73  DX1=0YA 

270  0YY=0X1 

OX  Y = UX 1*0Y Y 

269  SUH=SUH+  (REALIFKI  ) ) * F lit  -A  IMAG(  F2  (1)  )*  FUN+I)I*DXY 
SUM  1 = 1 . — SUM/ t A*MRR ) 

WRITE  16, 35)  SUM1 
i70  A=A+OA 
STOP 
END 

SUBROUTINE  SIMg(A,B,N,KS) 

DIMENSION  A ( 11,311) 

TUL*0.0 
KS  = 0 
JJ  = -N 

OU  65  J=1,N 
J Y = J+  1 
JJ=JJ+N+l 
t>lGA=<) 

1 T= J J-J 

00  30  1 =J  , N 
I J * I T + 1 

IFtASSI  iilGA)— A8S(  At  IJ>)>  20,30,30 
20  3 1 GA=A 1 I J ) 

1 M A X = I 

30  CONTINUE 

I F ( A JS ( 8 IGA ) — TOL ) 35, 35, AO 

35  KS=1 

RETURN 

AO  n=J  + NMJ-2> 

IT=1MA<-J 
DC  50  K=J,N 
11=11 ♦N 
12=11+17 
S AVt.  = A ( 1 1 ) 

A(  1U-A  1 121 
A < 12)  = SAVC 

50  A 1 1 1 ) =A 1 1 1 I /B  IGA 
S AVt  = l>  l 1MAX  ) 

A t IMAX ) =8 ( J ) 

H | j 1 = SAVt/UlGA 
1FIJ-N)  55,70,55 
55  H)S=N*(J-1) 

DO  65  1X=JY,N 
1XJ  = 1 jSMX 
I T =J- 1 X 
DU  6C  JX=JY,N 
IXJX=N* tJX-1) +IX 
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J JX=I X JX* IT 

60  A( IXJX  >=M lXJX»-( At 1XJ)*AI JJX)  ) 
65  BUX»=3tlX)-lBt  J)*A(IXJ)  » 

70  NY=N-1 
I T = N*N 

DU  ao  J=1,NY 
I A = I T — J 
Ib=n-j 
I C =N 

DO  00  K = 1,J 

5 l i B ) =B  1 1 B )— A ( IA)*B(ICI 
ia=ia-n 
eo  icsjc-i 

RETURN 

END 
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SAMPLE  INPUT 




SHALL  CRACK  LENGTHS 

FMOip+ee* 

'’“'GY 

or 

■’t17ilC(ULA3 

*u. 

AD»  iOl  k-o. 

Tel  0«&M  mo. 

A h Hl V ! f 

0*  Tf 

3 36000000  1 1 

i-Li  li  i r i i a i 

i i i \ i i 1 1 j 7 
11(117  110  1 

Lt.Lf  l.T.I  * 4 11  7.1.1  0 1.1.1,*.!  I T t 1 0 7 1 J 1 1 * 7 t_Lf.  K 1.1.  * 1 1,1.1  1,11.1  1 * 1 1 7 B 1 0 

1 3 - 3 E + f> 

, 3 3 

. G & 3 

io  * + e 

,33 

, -063  , 

..  0...6E  + S. 

0 0.8 

0 . 0 

2-u.L- 



a— * ! 

OUTPUT  FOR  SAMPLE  INPUT 


CLACKED  LAYEK-THICKNLSS  = C . 06 30  E = 
SOUND  LAYER-THICKNESS  = 0.CI630  E = 0 

ADHESIVE  - ShEaK  MODULUS  = 0. 60006  + 05 
NO  DE9UND  IN  THt  AOHkSIVE-DcB  = 0.0 


Y = 

0, 1250E-01 


X a 

C.L250E-01 

0.37505-01 

0.6250C-G1 

0.S7SCE-01 

0,11256+00 

0.13756+00 

0.16256+00 

0 . 1 8 75E  +00 

0 . 1 2 5 05 -0 1 

0.37506-01 

0.62506-01 

0.87506-01 

G.1125E+C0 

0 . 1 3756  + G0 

0.162SE+00 

0.18756+00 

0.12  50E- 0 l 

0. 37506-01 

0.625CE-01 

0.3750E-01 

o.uzse+oo 

0 . 1 375  fc  +O0 

0.162  St +00 

0.13756+00 

0,12506-01 

0 .3  760E-0 1 

0.62506-01 

C.376CE-01 

0.1 1256+00 

0.1375E+C0 

G.1625E+OU 

C. 1675 6 +00 

0.1 2506-0 1 

0.37506-01 

0. 62506-01 

0.C75CL-O1 

0.1125E+00 


0. 1*50E-01 
G.125GE-01 
0.1250E-01 
0.  12  50c -0 i 
0.12506-01 
0.12  5CE-0 1 
0. 1250E-01 
C.375GC-01 
0.2750E-01 
0.37506-01 
0.3750E-01 
0.3750E-01 
C. 3760C-01 
0. 37iuB-0l 
0.2760E-01 
C • 62606-0 1 
0.6250E-01 
C. *2506-01 
0. 62506-01 
G.6250C— 0 1 
0.6250E-01 
0.62  506—0 1 
C. 62506-01 
0.6750E-01 
0.67506-01 
0 • P750E-01 
0.27606-01 
0.7750E-01 
0.8760E-01 
0.37606-01 
G,  £7  5CiE  — 01 
C.1125F.+CC 
0. 112  6E+GG 
0.  LI  25  E + 00 
0. 1 12  5 E + 00 
C-.  1 125E+00 


TUX  * 

C.2273E-02 
0 .66 13E-02 
0. 1C6AE-01 
0 . 1 c 1 6 c —0  ) 

0. 1L53E-G1 
C.9Lo9E-J2 
0. 1&32E-01 
0.1267k— 01 
0.160  L 6-02 
0.  J303E-02 
0.LAA56-0L 
0. 19136-01 
0.2 1 9 1>  t —0 1 
G.2AAL6-OL 
0 » 2 7 ti  A E — 0 1 
U . 3096L-G 1 
0.122^6-02 
0.96956-02 
0.  L763L-01 
0. 23531-01 
C.2ti5')E— 01 
0.32  5A6-0 1 
0 .36 17t— 0 1 
0.3795E-0 1 
0.1 6776-03 
0.98096-02 
0.19216-01 
0 . c 56  6E-0 1 
0.01blifc-01 
0 . 362 1 E-0 1 
0 . 3903L-0 1 
0.3  fcA2  L— 0 1 
-0  .A  1 56E-0 3 
G.9999L-G2 
0.201  lb-0 1 
0.2672  fc— 0 1 
0.3231 t-0 1 


C.  1036  + 06  NU  a Li. 3300 


1026+0(1  NU 
T H I C KiNE  S S 

7UY  a 

0. 173A6+00 
0.17266+00 
0. 16566  + 00 
0.  16O0R  + C0 
0. 1AAA6+  00 
0. 1207E+00 
0. 100aE+00 
0 • 5 6 A7F. -0 1 
0. 15a36+00 
0.15256+00 
0. LA67E+00 
0.  1371U+00 
0.  12566+00 
0. 11026+00 
0.66336-01 
0.  5w23K-&l 
0.13656+00 
0.  1357E+00 
0.  12  69  6 + 00 
0.  1201E+00 
0. 1065E+00 
0 . 9394E-01 
G.735VC-01 
0. A63ab-01 
0. 1201E+00 
0 • 1 1 95E  + 00 
0.  1 1306+00 
0.10A7E+00 
0.9390E— 01 
0. ->0726-01 
0.6A5A6-01 
0.  **6  536-0 1 
0.  10526  + 00 
0.  luA7E*Q0 
0 . 93 AA6-01 
0. 91236-01 
0. 6 1C  76— Cl 


0.3300 

0.0060 

A = 0.2000 
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C.l  375E+U0 
0 .16256  + 00 
0 • 16  752  *00 
0.1250F-61 
0.3750E -01 
0 .62E0E-G1 
*5  71*  t,L  “1/ 1 
C.l  list  + 00 
0.1375E+00 
o.ifcrst+oo 
G. 18756+00 
O. 12506-01 
C. 37506-01 
0.625  06-01 
0.6  750t -0 1 
o.H256+no 
C. 13756+00 
0. 16256  + Go 
0.16756+00 
0.12502-01 
0.37506-01 
0.52  506—01 
0.37506-01 
0.11*56+00 
0.13756 +00 
C. 15256+00 
0.16756+00 
0.22506+00 
0.275Gt+00 
0.22506+00 
0.27506  +00 
0.22506  +00 
0.27506  +00 
0.2  2 506  +GG 
0.2750E+G0 
0.2 5006-01 
0.7  5006-0 1 
0.12506+00 
0.17506  +00 
0.25006-01 
0.76006-01 
C . 1 2 506  +00 
0.17506+00 
0.2  5C0E— ill 
0. 75006-01 
0.12506+00 
0.1 7506  +00 
0.2  5006  -01 
0.75006-01 
0.12506+00 
0.17506+00 
0.25006-01 
G.75U0E-01 
0.12506  +00 
0.17506  *00 
V ..  5l‘i.  L --O  1 


0.  11256+00 
0.  11256  +00 
0.  i 12  56  + OG 
0. 13756+OC 
0.  13756  + 00 
0. 13756+00 
0.  13  752  + 00 
0. 13756+00 
0.13756+00 
0.  1375  E + 00 
0. 13756+00 
0. 1&25E+00 
0.  10256+00 
0.  lt.256  + 00 
C. 16256+00 
0. 16256*00 
0.10256+00 
0.1625L+OG 
0. 1 6 25  E *00 
0. 12 75E+OC 
0.  U75E  + 00 
0.1b756+00 
0.  lij 756  +00 
0.  IE75E  + 00 
0. la  75 6 + 00 
0.16756+00 
0.13756+00 
0.25006-01 
0.25  006— 0 1 
0.7500E-01 
O.75C0E-01 
0.  1250E  + 00 
0. 12506+00 
0. 1750E+0G 
0. 17506+00 
C. 22  506  + 00 
0. 22  5 dE  + 00 
0.2.506+00 
0.22506+00 
0.27506+00 
0. 27506+00 
0.27506+00 
0. 2750E+00 
C. 32606+00 
0.32506+00 
0.32  306+00 
0.32506+00 
0.  37 ‘>06  + 00 
C. 37506+00 
0.3750L+00 
0.37506+00 
0 . *2  50E  + 00 
0. *2506+00 
b*^  but *00 
wsor  *og 

i / Ml  L *uu 


Cu37S4i>0l 
0. 3 bout —0  1 
0.392*1-01 
-0.1  5**E  -0  2 
0.92556-02 
0.l969t-01 
0.26786-01 
0.32966-01 
0.36326-01 
0.38*36-01 
0.381*6-0 l 
-G.237dE— 02 
0.86786—02 
0.1769L-01 
0 .2  58*»t— 0 1 
0.30866-01 
0.3**66-01 
0.3/7*6-01 
0.37506-01 
-0.250*6-02 
0 .760 1 L— 02 
0.16696-01 
0.23676-0 l 
0.28376-01 
0.33106-01 
0.359*6-01 
O.  27*0E-G 1 
0.15606-01 
0 . 1603C-U 1 
0.33206-01 
0.25766-01 
0.35526-01 
0.31566-01 
0.36b  36—0 1 
0.23806-01 
0.532*6— o2 
0.17836-01 
0. 28906-01 
0.3**56-01 
0.3357t-02 
0.1*686-01 
0.2*606-01 
0.3 19  IE— 0 1 
0.9759L-03 

0.1208t-01 
0.22006-01 
Q.2813E— 01 
—0 .116 * 6—0  2 
0.76516-02 
0.176*6-01 
0. £ 3956— 0 1 
— 0 . 5**0L— 02 
0*51926-02 
0.15336-01 
O. *1966—01 
-o.  / ,>«j  n -ot 


0. 696*6-0 l 
0. 569o£  — 01 
0. *2  606-0 1 
0 . 9 l 7ot  —0 1 
0.915*6-01 
0.a587£—ui 
0.79-.76-01 
0. 708*  E— 01 
0,81056-01 
0. 5077E— 01 
0. *.0176-01 
0.79936-01 
0.79896-01 
0. 7592E— 01 
0.69*26-01 
0. 02166-01 
0. 5*23E-0 1 
0. *5506-01 
0.27176-01 
0.69*26-01 
0.09776—0 1 
0.66336-01 
0.61 1 5E— 01 

0.5*876-01 
0. *336 6-01 
0.*080E-01 
0. J*31E-C1 
0 .*6  3ofc— 02 
-0.63296-03 
0. 16226-01 
0.32996-02 
0.23  >7E-u l 
0. 75*76-02 
0. 2*28  E-01 
0.10926-01 
0.573*6-01 
0.52606-01 
0.*2  5*6-01 
0. J255E-01 
0. *3696-01 
0. *0766-01 
0. 3‘»05c— 01 
0.28*06-01 
0.3339L-C1 
0.32066-01 
0.2687,6-01 
0.21*16-01 
0. 25*7t-ol 
0.25626-01 
0. 21286-01 
U.  1 oo*  E-0 1 
0.19656-01 
0.20 J9E-01 
0. 16  *2  6 — 01 
0.1*5t.E-01 
U.IM  96  — O 1 
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U*  /OOUt  - jl 

0 » \ t l*  U i + w u 

o.: 4 ^c.^j 

^*f.  LwW 

o.r>u.o;  -o; 
o.i v 'J'*>  U 
C.SLGO*:  -C  1 
G.l  'wuf  O^ 

0.5CUGL-0 1 
0-1  6909 + 'j0 

I • F 


t 

t ^ T -"jL  + 00 

0*~/60t.  +00 

: - roc  r.-  -oi 

0 * 1 i*  v 0 +00 

L * 6 6 ^0  ^ + Gu 
v , V>  i'tC  z.  + 00 
0*6  0005:  + 00 
0 • 06  uCE  + GC 
tu  / 6u0'  + 0^ 
f),  /*  ._  09  + 00 


0*2  lb  2t  — O 4! 
wl  UIl-01 
0*1  /^2t-G  1 
0 • ^ 4 0 1 t “"C  1 
V.3610L-G  1 
0 * 1 t.~U2 

'J.liOSC-Cl 
0,4.02Vl-u2 
o * f co  2t “02 
0*3972^-02 
0.S7j0E^02 


u*  Oi 

0.1^19^-01 
C- 12AAt-ol 
0.  Jw*E-G2 

0,  i^dL-01 

0.  1102  lL— 

0. 1C20E-01 

u.tHHc-O*. 

0 * Oft to  l — 02 
U*6^^t;-u2 

0 *70Q7t-C? 


C.7£jiGM^t  + O0 


APPENDIX  D 

TWO  ADHESIVELY  BONDED  LAYERS  WITH  A CRACK  IN  ONE  LAYER 

(LONG  CRACK  LENGTHS  0.5  ^ A ^ 1.0  INCH) 
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INPUT  PARAMETERS 


The  computer  programs  have  the  following  inputs: 

1.  modulus,  Poisson's  ratio,  and  thickness  of  cracked  layer 
(El,  NU1 , HI) 

2.  modulus,  Poisson's  ratio  and  thickness  of  the  sound  layer 
(E2,  NU2,  H2) 

3.  shear  modulus  and  thickness  of  adhesive  layer  (MUA,  HA) 

4.  ratio  of  minor- to-major  axis  of  the  elliptical  debond  (DEB) 

5.  ha  If -crack  length  A (0.6  s AS  1.0) 

6.  number  of  increments  in  half-crack  length  (NAA) 

7.  increment  in  half-crack  length  (DA) 


DESCRIPTION  OF  OUTPUT 

The  outputs  of  the  programs  are  shear  stresses  (t  , T ) at  various 

x y 

locations  in  the  bonded  region,  and  the  stress  intensity  factors  (K/c/na). 


LISTING 


c 

C PRO  GR  A v CEVELCPEC  8V  M.M.PATWANI  OF  NCRTHPCf  CCRPCRAT  I CN , 4 1 RCR  A F T DIVISION 
C HAWTHORNE. CALIFORNIA, TELEPHONE  1213)  970-5285 

C PROGRAM  COMPUTES  STRESS  INTENSITY  FACTORS  AND  SHEAR  STRESSES  IN  TWC 
C PLY  AC-ESlVELY  5CN0E0  STRUCTURE  WITH  A CRACK  IN  CNF  LAYER.  THE 
C PROGRAM  HAS  CP T I C N TC  PLT  MC  CE80N0  IN  THE  ADHESIVE  CR  ELLIPTICAL 
C DEMONS  WITH  HINCR  TC  MAJCR  AXIS  RATIC  OF  C.1,0.2  CR  0.3 

c this  program  gives  converging  results  for  half-crack  lengths  between 

C C. 5 and  l.C  INCHES 

C TELESCOPIC  GRID  ts  USFD  IN  THE  INTEGRATION  CF  KERNELS  »H1CH  i-AVE 
C LhGAR  1 1HIV  IC  SINGULARITIES 

C The  PRCPLFM  IS  RECuCEC  TC  SOLVING  A SET  CF  INTEGRAL  ECUATICNS  IN  WHICH 
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c SHEAR  STRESSES'  I LX  AND  TUY  ARE  TAKEN  AS  UNKNOWNS 

'C  SOLUTION  CF  UNKNON  SHEAR  STRESSES  TLX  ANC  TUY  IS  ACCOM  PL I SHED  BY 
C SCI.  VINT,  A SET  OF  SIMULTANEOUS  EQUATIONS 

C X t Y ARE  THE  CO-CPC  I N A 1ES, FROM  THF  CEN TRE-LI NE , AT  WHICH  SHEAR  STRESSES 
C A3  F CCMPUTEO 

C S.I.F  IS  THE  STRESS  INTENSITY  FACTOR  IN  THE  OUTPUT 

C NX l » NX  2 ,f(  Y I ,N Y2 , ARE  THE  NUMBER  CF  COLLCCATI CN  POI NTS  IN 

'C  THE  TELESCOPIC  CRID  FCft  NUMERICAL  ANALYSIS 
C THF  SHE  OF  TELESCOPIC  GRID  DEPENDS  CN  CRACK  LENGTH 
C INPUTS  ARE  - El,MJl,Ht,E2,Nli2,H2,MUA,HA,D£E»A,DA,NAA 
"C 

COMPLEX  Z ,7  2 ,ZB  ,Zfi2,ZBS, ZS.ZZZ, Z 1 .2  18  ,ZJB  ,PH21 ,PHZ2, PH 23, PH1.PH2, 

3F1,F2,TH21.  TH22,2Z,ZL4,ZZB,ZZ2, ZZ2A, ZZS.ZZ82,ZZ2B»ZZBSiZPH, ZPHB, 
3ZZZZE,ZZI,Z21B.ZZJB,ZZJ,TMll,TH12,THl3,TM14,TH31*TH42,221.Tnl,TH2, 

3TH32* THAItZ IZ.Z2Z.ZL1 tZL2tZ3ZBt Z A ZB  , ZFL  ,ZPLB  , PHZ4 
3.ZXY.ZAZ, ZlZE,Z2Ze,ZL3 

3,ZJ,ZIM,ZIBM,ZJM,ZJBM,PH3,PH4,ZPHP,ZPHEP,ZM,zeM,Z36M,ZBMM,TMll, 

3 TM 12, TM 1 3 ,T  M 14 , ZPHPB , 2PHPP, ZBHP  ,ZBPP  » TM31  ,TM32,TM33»  TM34#Th33  » 

3 T i-  34,  Z 7 I M , Z Z (BM  , THAI,  Th4  A ,TH5  1 * TH52  » TH53,  TH54,  ZZJBM,  ZZJM 

REAL  NL1.NU2.MUA 

0 1 MENS  ION  Z X Y I IAS) ,F( 298) ,C B I 88 8 04 ) , F 1 1 149 ) ,F2 ( 149) 

A5  FORMAT  ( 2X,  It,«E  13.  A) 

50  ” FORMAT  IE  I C. 2 » 3 F 1C.AJ 

51  FORMAT ( 2F  10. 3,  SXT  12) 

76  FCRMATiJX.8EI3.AJ 

87  FCRMAtTx,’ 25H  SOUND  L A YER- TH ICKNE S S = ,F7.A,2X.5H  E = , E9.3.3X.6H  N 
3U  = '.F7.AI” 

88  FCRMATIX. 27 H CRACKED  L A YER- TH IC KNE SS  * ,F7.A,2X,5H  E * ,E9.3,3X,6H 
3 NU  * . F 7 .4  ) 

89  "FORMA  Tl  X,  28H  ACHESIVE  - SHEAR  MCCULUS  = , E 1 0 * A , 5X . 1 3H  THICKNESS  = 
3, F7.A ) 

ICO  _ FORMAT ( 2 X , 1 LES.3) 

85  FORMAT  IX, SH"  S.I.F.  = . E15.‘7) 

86  FORMAT! 2X.5H  X * .8X,5H  Y = .8X.7H  TUX  * ,6X,7H  TUY  = ,10X,5H  A * 
3, F7  .A  I 

90  FCRMATIX. 33H  NC'CEBONO  IN  THE  ACH E S I V E-CE B = . F 5 . 3 I 

91  FORMAT! X. 58H  ELLIPTtCAL  OEBOND  IN  ADHESIVE-MINOR  TO  MAJOR  AXIS  RAT 
310=  , FA.  1 I 

CRACKED  LAYER  PARAMETERS 
El  = YOUNGS  MCOLLLS  CF  CR  AC  K EC  LAYER 
NU I = “CtSSCNS  RATIO  CF  CRACKED  LAYER 
H_l  _**  THICKNESS  CF  CRACKED  LAYER 
■"READ  (5,50  El.NUl.Hl 


**ts(l***fc*»****<I*»**» 


SC UNO  LAYER  PARAMETERS 
E2  = YOUNGS  MODULUS  CF  SO  UNO  LAYER 
NU2_=  POISSONS  PATIO  CF  SOUND  LAYFR 

>2"  = Thickness  cf  sound  layer 

R E AC  15,50  E 1 .NU2.H2 


C_ 

C**«**a***e**4**e 


, * * * , • A * 


C ACHESIVE  PARAMETERS 

C_  MUA__=  SHEAR  MCOLLLS  OF  ADHESIVE 

c ha  * Thickness  cf  adhesive 

C CES  = RATIO  OF  MINCE  TO  MAJOR  AXIS  CF  ELLIPTICAL  SHAPE  OF  CEBOND 
C DEB  = 0.,  0.1,  C.2  CR  0.3 
RE  AC  l 5, 5CT  MUA  , HA, DEB 


C 
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C "'CRACK  LENGTH  PARAMETERS 

C A = HALT  CRACK  LENGTH 

C NAA  = NUMBER  QF  INCREMENTS  IN  HALF  CRACK  LENGTH 

C CA  s INCREMENT  IN  H£Lc  CRACK  LENGTH 
READ  15,5  1*  A , C A , N 4 A 

C 

Pt=3* 1V15S2 
Pf2=2.*PI 
P 122=91/2. 

P K l-n.-NUl  )/(  I.*NUl  * 

PKZ=13.-NL2 1/t  l 
PPK  I^PI  2*1  U+PK13 
P p K 2 ~ P £ 2+  ( L • ♦ P X 2 * 

PKS=PKi*PKl 

PKM*|PKl-U  1/2. 

FKIC=FK 1/2. 

' G1  = E 1/12.*  ( I.  + NUUI 
G 2=E  2/ ( 2 .*  I 1.+M2J  ) 

C4^4.*GI 
C l 2=  2 • *G  1 
G22=2.*G2 

PKZ  = PKl~  W 

PKP-PK'1M« 

H2R=K2*PPK2*G22 
HRfi-F  1*PPX  i 
1R  = HRR  *G1  2 
EG = HA/MU A 

P I N = I ,-P I /2  * 

WP  ITE  (6,68)  HI # E 1 , N U l 
WRITE  (6,97)  H2,E2,NU2 
WRITE  16,691  M U A , H A 

c the  setting  of  ix,y)  co-ordinates  of  telescopic  g^ic  starts, 
nab-n  a a ♦ i 

C_  THE  FOLLOWING  DC  LCCP  VARIES  THE  CRACK  LENGTH 
CO  270  IMKJ - 1 , N AB 
NX  1=  16 
N V Is  1 
N X 2-  1 
NY  2=  4 
NX?=6 
N Y 5 
K X 4 = 1 
NY  A-  2 
N X 5 = 3 
NY5  = 2 
A 2 = A * A 

OY1=a/(2*NX11 
C Y 2“ 2 * *DY  L 
DY  ? = 2 • * D Y 2 
CY4  = 2 ,*DY3 

JR  (CEfi-O.  1 ) 570,571, L69 

169  IF  (DEB-0.2  ) lit, 170, 171 

170  NR2  = NXl-l 

WR  I TF  ( 6, SI  ) CEB 
NR3*NX 1-6 

NPS^2*(NXUU-NP2-NR3 
NR  5 ” 1 
NR  4 = 1 

NY  UNYUNR4  + NR5 


144 


NXY*NX  1*NY  1 + NPS 
EP*DY  t*DY2 

00  ISC  IKL=NR2,NX1 
IKP=1KL-19 
EX=I2*IKL-1 ) *D  Y 1 

1 80_  2 X Y t IKP)=CKPLXIEX,EP) 

EP  = 'FP+DY2 

DC  131  IKS=NR2,M1 

1 P P * I KP* l KS- 9 

E X = I 2* 1 K S-  1 l+CYl 
181  ZXYI  IPP )=C*PLXf  EX.EP) 

_ GC  TC  183 
171  NR3=NXl-2 

WRITE  16,91)  0E8 
E P = CY  1+DY2+CY2 
NR9*NX 1-9 
NR  5 = N X 1-9 

DC  139  (KL*NR3*NX1 
I K P = I KL—  1 3 
EX= I 2*  IKL-1 ) *C Y 1 
189  ZXYI  IKP  ]*CMPL  XI FX,FP  ) 

NPS=3*(NX  1 + 1 J-NR3-NR9-NR5 
NXY=NX  l*NYl«NPS 

SP  = EP  +DY2  

DC  I8C  IKS=NR9,NXl 
(PP- I KP+I KS- l 1 
EX=| 2* IKS-L  1 *CY  1 
188  ZXYI  IPP  ) =CYPl X l EX  ,EP) 
EP=EP»CY? 

00  186  IKL=NR5,NX1 
tKP=  | PP+ I XL -6 
E X = I 2*  IXL-1  ) *CY  1 
186  Z XYI  IKP  )=C*PIX(  EX  ,EP  ) 

GC  TO  183 

570  EP-CY1 
NR  1*  1 
NR  2*  1 

NR  3=1" 

NR9  = 1 
NR  5 = 1 

WRITE  16,90  OEB 
GO  TO  572 

571  FP  = D Y 1 

' NR  1*NX 1 

NR2** 

NR  3 = l_ 

' NR  9= 1 
NR  5=  l 

WRITE  (6,81)  CEB 
5 72  ~ NP  S*  2* INXlf  1 ) -NR l - NR  2 ' 

NY  1 = N Y l *N  fl  3 +NR9+NR5 

NX Y=  N X 1* NY  I *NP S 

DC  87  3 I =NP  1 ,NXl 
EX* 1 2*1-1) *CY  l 

873__  ZXYI  I ♦ 1-NR  1 ) =CYPLX(  E X,  EP  > 
‘ EP  = EP*DY2 
NR  3 = NXl+  2-NR  1-NR2 
DO  6 79  I *NR  2 ,N  XI 
EX* ( 2*  I-  1 ) *CY  l 
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6 74  2XY  I I *NR  2) =CVPLXf EX, EPI 
_ l 8 3 _J7 C 57  2 J*1,NY1 

NP=  1 J-l  ) *N  X l*NFS 

ey=ep*j*cy; 

OC  57  ? I = 1 ,NX  1 
' NC  = N P ♦ I 

£X  = t 2*  I-  1) *CY  1 
_ _ 2XY  ( NO UCMPl X I EX.EY) 

5T5  continue 

EP  = EX  *cv  HCY2 

00  5 74  J ■=  l • N Y 2 

N S=N  XY ♦ l J- 1 ( ♦ N X2 
F Y = f 2*J- 1 I * C Y 2 
DO  574  J = 1 , N X 2 
NM=NS»I 

EX  =E  P + ( I - 1 ) *0  Y 3 

574  2XY(NM)=CVPIXIEX,EY) 

JXY=NXY*NY2*NX2 

1 X Y = J X Y + N Y 2 * N X 3, 

E P-E  Y *0 Y 3 

CO  575  J= 1 , N Y 2 

EY-EP  + IJ-l  ) *0  Y 2 
N P = J X Y + ( J — 1 ) + N X 3 
DO  575  [ =1  ,7.  X3 
"EX  = I 2* I- t) *CY2 
NG^NP* I 

575  2XYINQ  )=CMFLXIEX,EYI 
EP=EX+CY2+CY3 

CO  576  J = l , N Y 4 

E Y=  I 2*  J-  1 ) *03 

NP=  IXYM  J- l >*NX4 
CO  576  I = 1 , N X 4 

NCf  NF*  I_ 

EX  = FP+[  I-1MCY4 

576  ZXY{ NO )*CMPLX( EXtEYJ 
K X Y“  IXY4-NX4YNY4 

E P = A I^AGIZXYt  IXY) I*DY2*DY3 
KN  = K X Y *N  X £ *N Y 5 

00  _5  7 7 JM  ,NY5 

E Y * E P +"(  J—  1 J » C Y 4 
N P = K XY ♦ I J-  I ) *N  X 5 

00  5 7 7 _I  = l , N X 5 

NO -N  P ♦ I 

EX= I 2*  I-  1 ) *CY3 
_ 577  __ 2 XY( NC I -CM  PL  XI EX,EY) 

C Th-E  setting  OF  TELESCOPIC  GRID  ENOS 
NXN2  = 2*KN 
KNN=NXN2*KN 
KNFsKNNfKN 

EVALUATION  CF  KERNELS  K 11 , K 1 2 , K 2 L , K22  CF  INTEGRAL  ECUATIGNS  45  OF 

SECTION  2.2.3  STARTS 

C FOLLOWING  LOOP  IS  CUTER  LCCP  FOR  EVALLATI CN  OF  KERNELS  VARIES  2 

OO  6 5 I»  1,  KN 

IP=  I 

I C =K  N* I 

2-ZXYUI 

22  = 2*2 

2e=CCNJGI 2 ) 

2E2=2fi*2B 

2S=CSCPT(22-A2) 
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Z9S=C0NJG( Z S ) 

III  = Z-ZB 
12 l=Z/Z$ 

Z4Z  = ze/ZBS  _ 

F(IP>  IS  RIGHT  HAND  SIDE  CF  FIRST  OF  INTEGRAL  ECUATICNS  A5  OF 
SECT  tCN  2,2 .2 

F f 1C),  IS  RIGHT  HAND  SIDE  CF  SECCND  OF  INTEGRAL  ECUATICNS  CF 
SECTION  2.2.2 

FUP)*(-PKZ*PEAL  (ZXYII  > )*PKl*REAL  I ZS  ) -PEAL  I Z6S ) *2. * A IP  AC  I ZXY ( [))* 

3 A IMAGIZAZ I ) /G A 

FltC1-l-PKZ*Mf’«GtZXVm)«PKl*AirAGIZ!>-AIPAGIZBS>-2.*MN'AG(2XYtl> 
3)  *».EALI  ZAZ)  I/GA 
Z !=ZS-Z 
Z I E = CON JGI Z I ) 

Z IZ=A-Z 
Z ZZ  = A-ZB 
Z L 1 = A ♦ Z 
Z L 2 = A + ZB 
Z J=Z/ZS-1. 

ZJ0=CCNJGIZJI 
ZIP=-ZS*Z 
Z I BP  =CCN  JG  I Z IM  ) 

Z JV=ZJ 
Z J 2M  = Z J8 
PHZ  1 = I A+Z 1 1 /Z IZ 
PHZZ=( A*Z IB )/l2l 
Pi-22-  I A«Z  IM  l/ZL  I 
PFZA=IA+ZIBM)/ZL2 
_P  Hl= (PHZ  I+ZJI/Z  1Z 
PF2-I P^ZZ  + ZJB  I/Z22 
PH3={  PHZ3*Z  JN1  l/ZL  1 
PHA  = I PHZ  WJBMJ /ZL2 

ZPl3~PMZ  l + PK  1*PHZ2  + ZZZ*FHI+FHZ3-FM*PHZA+ZZZ*PH3 
ZPLe=-ZZZ*PH2-PHZ2FPKl*PHZl-ZZZ*EH4+PHZA'FKl+PHZ3 
Fill  )=ZPL+ZFLR 
F 2 ( I ) = ZPL-ZFLB 
C INNER  LOOP  FOR  EVALUATION  OF  KERNELS  VARIES  ZO 

CO  6S  I 1 = 1 , KN 

N 0 = I 1 1-U*NXN2«I 
NC -K  N+  NP 
N S = KNA  *NP 
NP-KNM+NP 

H=  zxyii  iT  *■ 

IF  ( I I-NXYI  66S.669.67C 
669  DX  I =OY 2 
GO  TO  671 

67C  IF  IIl-IXY)  472,672,673 

672  DX  1 = CY  3 
GO  TC  671 

673  OX  1 = D Y A 
GO  TO  671 

671  C Y Y=  0 X 1 

C X Y=  OX  l*OY Y 

IO=I  1 

I S = KM  II 
ZZB  = CONJGIZZ) 

_ 112  = ll*ll 
ZZ2A  = ZZ2-A2 
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ZZS  =C$C«nzZ3M 
_ ZZ92  = ZZP*Z  Z8 
Z Z 26  -ZZB2-A2 
ZZBS=CCNJG  < ZZS  ) 

_ZPH  =z~zz 

Z °H  B = CCN  JG ( ZPH  ) 

in  is  = u~ue 

zzi  = zzs-zz 

ll  IB^CGNJG ( ZZ  f ) 
zz J=ZZ/ZZS-1. 
ll Jft-CCNJG ( ZZ J I 
“ ZP(-P  = Z+ZZ 

ZPhBP=CONJG(ZPhPl 
ZPHP8-Z-ZZB 
ZPhpp=z*Z  ZB 
ZB(-p*CCNJG(ZPhfb  I 
ZBpp=CC\JG(ZPHFP) 

' ZM  = Z * Z Z 

Z BW  = CCNJG(  ZM 
IBB'*-!*?.  IB 

ZBWsCGNJG ( Zfipv ) 

Thl  ? = CLnc I Z^-flZtZZS^ZSI 
rHU  = CLCGIZF!i?M*2»ZZeS*ZS) 

Th  l l^GCNJG  t Th  111 
Th  12=C0\JG I TH \ t) 
rvi2-CLOG|-zev“-A2-Z8S*ZZS) 

T v 1 3 = CLDG  l-ZV-/S2-ZZ$*ZS) 

tv  i t=cc^jr.  I tv  i ri 

TV  l^CGNJG  {TV  12  ) 

P V 3 = CLCG ( ZPhP) 

PH  MCCN  J G ( PH  3 J 
P H Z 1 - CL  CG l ZPHPE  ) 

PH  Z 2 = CONJG  I Fh  Z 1 ) 

Ph Z 3 = CLCG  1 Z Pt-PP  I 
PHZA=CCNJG (PHZ  ») 

ZL  1*0  .5*1  l.-Z5*ZZ95/ZZ2B)*ZZZZB 
Th  3 l = -ZL  1 /ZPHPS 

TM3  l = C.5* { l.  + Z S*ZZSS/ZZ20 I ‘ZZZZe/ZPHFF 
TH33*crNJG  I T H 3 1 ) 

T v 3 3 = CON JG ( TV  3 1 } 

TV22=  0.5*1  L. *Z ES*ZZB S/ZZ2S ! *ZZ lid/ ZPHBP 
TM3^=CCNJG ( TV32) 

ZL1=0.5*ZZZ 
ZZ  I = ZZS-ZZ 
ZZ  IB=CGNJG  (ZZ  t ) 

ZZ  IMx-ZZS+ZZ 
ZZ  1BV  = CCNJG  f ZZ  IH) 
TH42=ZLl*(Zie-ZZ[BM)/lZPHeP*zqsi 

Th^3=ZL l*(Z IB-ZZI )/( ZBS*ZBHP) 

TH44*  ZL  l*( Z le-ZZIV) /( ZBE*Z8PPI 
ZL  1 = ZL1/ZBS 
ZZJ= l Ll/ll S-  1 . ) 

ZZJ8=CCNJG ( ZZJ  ) 

zzjbv=zzjs 

ll JV=ZZ J 

TH51=(TH43-ZLI*ZZJ)*ZZZZB/ZBHP 
T h 5 2 3 “ I TH44-2L  1 *ZZJM  *zz  ZZB/ZBPP 
TH  5 4 = ( T H 4 2-  2 1 MZZJBM)  *ZZZZB/ZPHBP 
TH22=CL0G( Z*ZS  ) 

TH1=CLCGIZ-ZSJ 
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Th2  l = PK  1 l* ( PKL  *TH2  2-CCKJG(  TH22  )-PK  l *TH*CGN JG(  TH1 ) 1 
1E  ! I.eQ.I  II  GC  TO  77 
Pi-  l=CLCC(  IP*-  1 
PH2  = CCNJGt  Pf- 1 ) 

TH2=Z°H/ZPHB 

Th32  = -0.5*{  l.-ZPS*ZZB5/ZZ2ft)*ZZZZB/ZPHf! 

T H 34 =CCNJG l TH3  2 I 

TH53  = -ZZZZBM  TH41-ZL 1*ZZJB)/ZPH8 
TH4  l=ZLl*  1 1 Ifi-ZZ  IB)  /ZPH8 
GO  TO  78 

IT ' T t-  4 1 = C.5*ZZZ*(ZE/ZF3$-1«  ) / Z B S 
TH32=C.S*ZB*ZZZZP/IZBS*ZZ0S) 

I>?<i=Cr,NJG  [ Th’2  ) 

THS3=C.25*A2*ZZZZ8*ZZZ/(  I ZB2-A2  ) MZB2-A2)  I 
0X2= (DX 1*DX  1 1/4. 

0X22=2. *0X2 
' DXP  I =DX22*P I 
TH2=CVPL  X( CXPl ,C« I 
TH?  = TH2/DXY 

' PH 1 = DX22*<  ALCG (0X22) *PI 22-3. ) 

PH  1 = PH  l/DXY 
PH2=PHl 

78  ZL  l = - I PH  l*PH  2 ) +PH3*PH4 

Z12=ZPHPB/Zet-P-ZPHPP/ZBPP 
ZL  X=-l PHZ  t*PHZ2)+PHZ3*PHZ4 
ZL4=TH2-ZPHP/ZPH0P 
Z l Z 0 = PK  2*  Zl.  I ♦ Z L 2 
Z ?.l  P = PK2*ZL  3 + Z L 4 

Z3ZB  = PKl*ZLI*ZL2*PKlD*ITHl3+THll-TM13-Twil)+0.5*(-TH12-PKS*rHl<t* 
3TM12*PKS*TM  14) *TH51.-TH52*PKl*lTI-34-TH24  + TH4i-TH42)“TH33*TM3  3-TH43 
3*TH44  + TH2  l 

Z4ZB=PKI*ZL3+ZL4*TH53-Th54+PX1D*ITH14*TH12-TH14-TM12)*0.B*I-THII- 
3PKS*Th13+TM 1 1+p«S*TM13) *PKl*l TH3l-TM3l*TH43-TH44)-TH32+TM32-TH41* 
3 TH42*Th21 

R 1 1=RE AL  ( Z 1 Z6  I *PEAL  I Z2ZB  » 

R 1 2 = - A I M AG ( l 1ZP  l+AIMAG(Z2ZB  ! 

R21=A  (MAGIZ  lZe)^IM4GIZ2ZB) 

' R22=RFALt ZIZe )-REALI Z2ZB) 

HI  l =R  E AL 1 Z 3 7 6 ) + RE6L I Z4Z0 I 
H12=-AIMAG(  Z 3 Z B ) *AI MAGt Z4ZB) 

~ h21=A1HAG(Z3ZP ) »A IMAGt  Z4Z0 ) 

H22=RE AL ( Z3ZEJ-REAL ( Z4ZB ) 

SI  1 = (H11/H1R*R l 1/H2R I *OXY 
S 12  = 1 H 12/H1R+R 12/H2R J *DXY 
S2t=(M21/Hl^+«21/H2R}*DXY 
S22=(H22/hir*r;;/h2R)*Dxy 


IF  < 

1 *EC  . 1 

r l j 

GO  TO 

79C 

789 

SU  = 

hG  * 

S 1 I 

%?2  = 

^0  + 

S 12 

790 

ce (np  t = sn 

C0(N<;)=S21 
CB (NS ) = S 12 
' C 0 ( N N )=S22 
b9  CONTINUE 

C EVALUATION  CE  KERNELS  ENOS 

C ' “SIHO  SOLVES  EDR  UNKNOWN  SHEAR  STRESSES  TUX  AND  TUY 
CALL  S IMQ ( C e ,F ,NXN2, KS) 

WR  ITE  ( fc»  86  > A 
00  88S  1 = 1 ,KN 
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889  WRITE  (6,761  Z XTM  [ ) , F (I)  ,f  ( KNt  II 
SUM^O. 

C FCLLChING  LCCP  CCWPUTES  STRESS  INTENSITY  FACTORS  FRGf  ECUATION 

C AB  CF  SECTICN  2.2.3 

CO  269  1 = 1 , KN 

IF  ( I-NXY)  268  ,268,271 

268  OXl=DY2 

GO  TO  27C  _ 

27 1 "IF  II-IXY)  272,272,273 

272  C X 1=  D Y 3 
GO  TO  2 7 C 

273  DX  l = OY  A 
270  CYY=OXI 

C X Y = C X 1 * DY  Y 

269  S(JM=  SUM*  (REAL  (Fill))*  F ( II  - A IM  AG  ( F 2 l I ) I * F I KN ♦ II) * DXY 
SUHl=l.-SU«/tA*HRR) 

WRITE  (6,_85)_SUY1 
370"  A=A+CA 
STOP 
END 

SUBROUTINE  S I M C ( A , B, N ,K S I 
D l^FNSICN  A ( l)  ,E( LI 
TC  L = 0 . C 
K S = 0 

J Jf_-N 

DC  65  J=  l'.N 
JY=J+ 1 
J J = JJ  * N*  L 
B IGA-0 

DQ_  3C  J»J_»N 

I j-  I T + I 

IF  I ABSIBIGA  l-ABSIAI  I JM  1 20,30, 3 C 

2 0 B I G A = A t IJJ 

"Iv'axM 
30  CONTINUE 

IF ( A8S( BIGA l-TCL ) 35, 35, AC 
" 3 5 ”KS*1 

RETURN 

AC I I = J + N*  I J-2J^ 

I T« I M AX- j 
DC  5C  K*J,N 
I 1=  I I + N 
I 2=  1 1 + 1 T 
S A V E = A ( I L I 
_ A ( tn=A(  121 
A ( 12  1=  S A VE 

5 C Al  I L l=AI  ill /g  IGA 
S A V 6 = F ( I MAX  ) 

01 IVAX)=B(J) 

El  J ) = SAVE/8  IGA 
TFIJ-M  5 5,  7 C , 5 5 
55  I 0 S = N * I J- 1 1 

DC  65  I X" J Y , N 

I X J = ICS* I X 

I T = J-  IX 

DC  6 C J X = J Y , N 

I X J X - N * 1 JX-  1 I *!  X 
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JJX  = IXJXMT 

60  " "At  tXJX]=A l 1XJX  }-( A I [ XJ) *A( JJX) ) 
65  BtIX)  = BUX>-  let  J>  *Al  I XJH 

7_C N Y = N-  1 

I T = N * N 

00  80  J=L,NY 
U=  IT-J 

IB  = N- J 
ION 

DC  80  K=1,J 

B ( 1 6 ) = B ( I3J-AI  I A ] * 0 ( I C ) 

1 A - IA-N 
80  IC=IC-1 

' ""RETURN 
END 
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SAMPLE  INPUT 


ids  TiTTl 

I.ONC  CRACK  LENGTHS 

tHOiNft* 

f*  Of 

HI 

N7*  fOD  hi 

*0*  O«0*t.  NO. 

anal 

- - . 

1 i“k  l 1 ■ 1 

QDOjpgeoei  t I t 1 } 
t UllHULt  Tll  9 0 1 i.  3 4 S 

l 

A 

f l > I J I J » ? J J J i 

■ *: » * > it'iit 

J 1 

O'  T 

) j I 1 1 1 1 i i « t t 1 i 1 i 

2 1LH.LI  1 O l J 1 * t 

1 I 3 1 M t 1 J 1 1 4 | M ( • I M 

t T T + 7 T 'T't, 

■7  1 J I f . t T » 9 0 

. , .1.0,,3.E  + 6 

. 3 3 

.06  3 

. A 0 - 3 E + 6 

..  3 3 

. 0 6 3 

0..6E.+  5 

. OQfl 

0..  i 

0 

Q.,6 

a.  o 

,0 

OUTPUT  FOR  SAMPLE  INPUT 

CMCKEC  LAYER-THICKNESS  = C.063C  E * 0.103E+08  Mj  = 0.3300 

SflUNO  LAYER-TMCKNESS  = 0.C63C  E =•  C.1C3E+CE  NL  = 0.3300 

Tohes i've -"shear ‘mcoulus  = c.6cooe+c5  ' thickness  *■  o.cceo 

FLLIFTICAL  CEECND  IN  ACHES  IVE-HI NCR  TO  MAJCR  AXIS  RATIC=  C.l 


X a Y 

- 

TUX  = 

TUY  = A * C. 6 000 

o.sene+oc 

"6TI 8 75E—  Cl~ 

c.AAeeE-ci 

C.7451E-01 

0.  1687E  + CC 

0.5625E-CI 

C * 29  7 C 1 

C.24876+C0 

0 • 2062  E + 0 C 

G.5625E-C1 

C. 3 7 2 7E -C  1 

0.23566+00 

0. 2437E+QC 

0*5625E-Ct 

0.4725E-C1 

C*2273E+Q0 

0 . 28 1 3 E +00 

0.5625E-01 

0*5C1€E-C  l 

C*  2 2 22  E + 00 

0.318 7E+CC 

C.  5625E-C1 

C*5<tdCE-Cl 

C.2133E+0O 

0.35626*00 

0.56  2 5E-C 1 

0.5S2CE-CI 

C* 1 993  E + 00 

0.393BE  +00 

0.562  55-  01 

C.5972E-C 1 

C*  1 8 7 7E+  CO 

0 . 4 2 1 26  +00 

0. 5625E- Cl 

0.6  1616-01 

C*  I 79  7E  + 00 

0*  468BE + CC 

C. 562 5 E- Cl 

0.626  6E-C  1 

“ C.  16C2E  + 00 

0.50625+0C 

Q.5625E-EU 

C. 64 1C6-C1 

C.  1426E+0Q 

0.54376+00 

0 * 562  5E - C l 

0.6  7566-C1 

C.  5 7 77 E- 01 

6.58136+00 

0*5625E-Cl 

0.66  19E-01 

C.  5842  6-01 

0. 1875E-C1 

0*9375E-C1 

-C.  37566-02 

C.  2C70E  + C0 

0.56256 -Cl 

0*9375 E- Cl 

0.7 076 E- 02 

C. 21  886  + 00 

0.93756-01 

" 0.9375E-C1 

C".  UA3E-C1 

C.2219E+00 

0 .13126+00 

0.9375E-C1 

0. 17546-01 

C.2215E+C0 

0. 1687E  + CC 

0 .9375E-  Cl 

0.26506-01 

C. 21596+00 

0.2062E  +00 

0 # 9 3 7 5E- C t 

0. 26276-01 

C*2073E+00 

0.2437E+00 

0*9375 E- Cl 

C .4  6 2CE-C1 

C*  1959E  + C3 

O.20I3E  +00 

C.9375E-C1 

0.4523E-C 1 

C* 1 9C5E  *00 

0 * 3 16  7 E * CC 

0*9375 E- Cl 

0.54266-01 

C.  1 82  OE  + 00 

0*3562E  +CC 

C.9375E-C1 

C.56666-C1 

C.  16 82  E + CO 

G*  3938E+O0 

0*937  5E- C 1 

0.61516— Cl 

C.  16006  + 00 

0 • 4 3 1 2E  + 0 0 

0.937  5F- C 1 

“ C .6  5 C 7E— C 1 

C.  14526+00 

0 . 46  88  E + CC 

0.9375E- Cl 

C.6656E-C  1 

C.  1 212E  + 00 

0 .5062E  +00 

0.93755- Cl 

0.7 1CAE-C1 

0. 1155E  + C0 

0.5437E+OC 

0.9375E-C1 

0 * 7456E-C1 

C.  6859E-C1 

0.5613E+CC 

0.9375E-C1 

0.7244E-C1 

C. 4654E-01 

0. 18756-01 

0 .131 2E  +CC 

-0.7  11  16-C2 

C. 17656+00 

0.5625E-01 

0.13L2E  + CC 

C. 5C77E-C2 

0. 18826  + 00 

0.9375E-CL 

0.131 2F  +CC 

0.62456-02 

C.  1915E  + 00 

0 . 1312E  + CC 

0.1 31 2E  +00 

0.  1 5 5 6 6-  0 1 

C.  19196  + 00 

0.  1687F  +CC 

0.1312E+CC 

C. 2 7 276-01 

C . 1 8 8 t £ + 00 

0.20626+00 

0. 13126  + 00 

C.36456-CI 

C.  18C4E+00 

0.24376+00 

0 .131 2E  +CC 

C.4  552E-C1 

0.16886+00 

0. 2813E  +C0 

0. 1 31 2E  +00 

C.454C6-CI 

0. 1 64  CC  + CO 

0.3187E+0C 

0. 131 2E  +CC 

C. 55046-01 

C.  1555E  + 00 

0 .3562E+0C 

0 . 131 26  + CC 

0.60  7CE-CI 

C. 14456+00 
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0 - 39  33  E + CC 
0 .431  2E+CC 
0 .4688  £ + C 0 
0 .5062E  + 0C 
0 .5437E+0C 
O.5013F+OC 
0*  1 0 7 5E- C 1 
Q.5625E-QI 
0*9  376F-01 
0 . 13  I2E  + 00 
0* 1687E+CC 
0 -2062E+O0 
0-2437E+CC 
0 * 281 3E+QC 
Q * ■!  18  ?E  + CC 
0 . 35*28 +0C 
0 . 3933E+CC 
0.431 2 E + CC 
0.-468  BE  + CC 
0 .50628  *00 
0-  543  7^  + 0C 
0.581 3E+0C 
0.  18758-01 
0.5625E-01 
0 *93758-01 
0 . t 2126+00 
0 . 1687E+C0 

0 -2Q52E  + QC 

0. 2437E+0C 
Q. 2013E+QC 
0.3187E+0C 
0.3562E+CQ 
0 - 39  30  E + CO 
0 -43 12E+0C 
0 -4688E+CC 
0 . 5062E+QC 
0*543 7E+CC 
0 - 53I38+CC 
0 .63766  + QC 
0 .6  3 756  *CC 
0 *6375E+CC 
0.6375E+CC 
0 .3750E-C1 
0 . 1 I25E+0C 
0, 1C75F+CG 
0-2625E+00 
Q-2375F+0C 
0 * 4 1 2 5E  +C  C 
0 .4675E+CC 
0*56256+00 
0 .375QE-CI 
0 . 1125E  + 0C 
0. 1075E+OC 
0. 26256  + 0 C 
0.2375E+CC 
0.4  1256+00 
0.4075E  +00 
0 ,56  256  + CC 
0* 17508- CL 


0*131 2E  + CC 
0*131 2E  + CC 
0-131 2E  +CC 
C - 13 1 2E  + C C 
0.1312E+CC 
C*  I 31 2E  + CC 
0.168  7E  + CC 
C * 16  B 7E + CC 
C.  16876  + CC 
C-  1687G+CC 
C-  16S7E  + CC 
0. 1687E+CC 
0*  160  ?£  + CC 
0 . 1 6 8 7 E + C C 
0*168  7E + CC 
0. 160  7E  + CC 
0. 16878+CC 
0.1687E+CC 
C*  168  7E  +CC 
0.168  78  + CC 
0*  168  7E  +CC 
€.  160  7E + CC 
O.2062E  +C0 
0-2062E  +CC 
C . 2062E  + CC 
C - 2062E+CC 
C.2062E+CC 
0.2062E  + CC 
0.2062E  + CC 
C. 206 2 E+CC 
0.2062E  + 0C 
C ,20626  + CC 
0 - 206  2 E + CC 
0.2062E+CC 
0.2062E  +0C 
0,20628 + CG 
0 .2062E+CC 
0.206 2 E+CC 
O. 3750 E- Cl 
0. 1 125E  + GC 
0.1875E+CC 
0.2625E+CC 
0.33  75E  + CC 
C . 3 3 75E  + C C 
0.3375E+CC 
0.3375E+CC 
0-33  75E  +CC 
C.3375E+C0 
0.3375E+CC 
0.3375E+CC 
0.41256 +CC 
C-4L26E  + CC 
0.4125E+CC 
0.4125 E+CC 
0.4125E  + CC 
C.4L25E  + CC 
0 .4  I25E ♦ CC 
C.4  12  5E  +CC 
0.4875E t CO 


C.62E5E-C1 
C . 6 6 4 4 E—  C l 
0.7  179E-C1 
C.7474E-C1 
C. 76846-01 
C.7177E-C  l 
C.1C14E-C1 
C.2U4E-C2 
0-4  126E-C2 
0*12 25 8 -Cl 
C. 2623 E- Cl 
C.2E44F-CI 
C*4  5CIE-C  1 
C * 5 C 4 1 E * C 1 
C.555CE-C1 
C-6CC  7E-C1 
0*66  74E-C1 
C- 7C78F-C 1 
0.7  364E-CI 
C.7402E-C1 
C.  75 5 EE- Cl 
C.  7 4 1 1 E — C 1 
0. 1286E-C1 
C*  ICC  IE-C3 
C-21636-02 
C.  1CECE-CI 
0 . 2 i 6 1 E - 0 1 
C.2659F-C1 

0*46  1 2F-CI 

C.44C5E-C1 
C.5  7C5E-CI 
C.6  2 15E-C1 
C.6E 65E-C 1 
C.7453E-C1 
C.  75  36E-C1 
C. 7 6 596- Cl 
C.76C9F-C  1 
C-7  2ECC-C  1 
0*4CC3E-Cl 

C .6  1 7 3E-C1 
C.66C3E-C1 
C.72C2E-C1 
C.74C3E-C2 
C.7224E-C2 
C.2C2CS-C 1 
C.4227E-C1 
0.558 3E-C l 
0.6266E-CL 
C * 71 5 7 E - C 1 
C.772CE-C 1 
0.1251E-C1 
0.9  1 1CE-C4 
C.14C6E-C1 
0.3  3 1 56-C 1 
C. 42  5 58* (1 
C.56C6E-CI 
C.6  76  7E-C1 
C.72  2 >r-ci 
0. 1 665F-C1 


C* 1338E+CG 
C. I L54F+00 
C-  I 071 E* 00 
C.  9277E-C1 
C*  71396-01 
C. 51 59E-0L 
C. 1512E+00 
C. 163LE+OQ 
C.1666E+00 
C*  1678E  + 0O 
C.  16  5CE  + 00 
C-  1641E  + 00 
C-  L401E+0O 
C- 14  I9E  + 00 
C-  L357E  + CO 
C.  12  796  + 00 
C. 11  19E  + Q0 
C.  1C12E  + 00 
C.9C40E-01 
C.  7656E-C1 
0*  62 1 9E-01 
C.4364E-Q1 
0.  13C4E  + 00 
C.  1426E  + C0 
C.  1465 F + QG 
C. 1432E+00 
C. 1468E+00 
0- 1454E+C0 
C. 1297E+C0 
C.  1 2 72  E + 00 
C*  1 1 54E  + 00 
C* 1C46E+QQ 
C-  595LE-01 
C. 83E6E-0L 
C.75C9E-0L 
G.6  356E-C1 
C. 531 7E-01 
C.41 26E-C1 
C - 2 722  E- 02 
C- 12638-01 
C.6376F-02 
C. 19688-01 
C*  EG58E-C1 
C.  51 1 5E - 01 
C.51S9E-0L 
C. 7724E-01 
C.6780E-01 
C. 5 74  4 E- 01 
C. 382BF-C1 
C.2624E-01 
C. 5667E-01 
C. 67698-01 
C.69  L2E-01 
C.626CE-01 
C-51 308-01 
C*  4238F-01 
C.28C7F-0L 
C- 1 732E -01 
C. 3557C-01 
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0 - 1 1 2 3E  OC 

0 * 4 6 7 4£  + CC 

-C.9  542E-C2 

C*  46E-C1 

o . ie  73F  + cc 

g.4B7se+cc 

C.  1 la&F-Cl 

C*  52S9E-0L 

0 .262FE  +00 

0.4675E+CC 

C-2  1 LSF-CL 

C.  5 35  IE-01 

C - +0C 

0.4S75F+CC 

C.5C£  IH-C1 

C. 3889E-CI 

0.4  L2^E  OC 

0*46  75E  + CC 

" c.^ecte-ci 

C*  3356E-01 

0.4P7SC+CO 

0 *48  75E  + CC 

0. 5 5 3 1E-C l 

0*2C@iF-01 

0 -662SF +C0 

0*407 5 C+CC 

C . £ E 7 2F- C l 

C.  1255E-GI 

0 . 3 7 S0E - C 1 

0.5625E+CC 

-C.2C  CEE-C1 

C*2eC5E-0l 

0 . 1 I 25  e +0C 

0 * 562  5E  + C C 

-0. 1 C 71E-C 1 

C*  3965E-01 

O*1875E*0C 

0.5625E  + 0C 

C.5S£5E-C2 

C-4328E-0L 

0. 26  25 F +0C 

0.5625C  ♦ C C 

c.  lectF-c i 

C.43C5E-C1 

Q , 3 3 75E +0C 

O . 5 6 2 5 E + 0 C 

0.355 2E-C1 

C, 3C1  7E-01 

0. 4 125E  *0C 

C .5625E+CQ 

C-^E  25E-C  1 

t*  27CCE-CI 

D.4P75F  + GC 

C.56Z5F+CC 

0.555  l^-Cl 

C- 1960F-0I 

O,56?5F+O0 

0.5625E+CC 

C.6t  7 3 6- C 1 

C.96  15F-02 

0. 2750f-Cl 

0 * 6 3 7 6E  * CC 

-C.22C6E-CI 

C* 1955E-01 

0 . M ?5E  + CC 

G.&375E+CC 

— C . 1 2 5 4E- C l 

0*31696-01 

0 * 18  7 5^ + 00 

0 * 6 3 7 6 E + CC 

0 . 1 U4E-C2 

C.  3 9 7 1 E - 0 1 

0 .26256  +0C 

0-6375E  + CC 

c.  ifssf-ci 

C*  3 5 35  5 - C 1 

0 *3375€+CC 

0*63  75E  + CC 

C.2S3SF-C  t 

C. 31 03E-CI 

0 .41  2SE+C0 

0 *63  7 5 E *CC 

t.<*  \ 5 7 E — Cl 

C*2463E-C1 

0 .4P75F  + 00 

0.637 5 F+CC 

0 . 5 i Uc- Cl 

C- 1820E-01 

0 .56256+00 

0*63  7 6 E + C C 

' C.fc55CF-Cl 

C.9243E-02 

0 .6750F  +CC 

0*  75O0F -Cl 

C.5t57E-Cl 

-C.6CC0E-C2 

0 * 6 7 50r + 00 

0 * Z 2 5 C E + CC 

C.6  525F-CI 

C* 7 l C l E-  03 

0. 75 COE-01 

0.750QE+CC 

" -0.1  5 1 3F-C1 

C*23  30E-01 

0, 225QE + 0C 

0-750CE  + CC 

C.9  552F-C2 

C.2977F-CI 

0 . 375CE+CC 

0*  750CF  + 0C 

0.3‘26F-C1 

C . 2 1 7 7E-01 

0 * 7 SCOT- 01 

C.90OCE + CQ 

-C. 1 lEtE-Cl 

C- 1 74GF -01 

0 * 2Z60E  + CC 

C .9QQGE + CC 

-0.5  1 -3F-03 

C.2915E-C1 

0 , 3750F+CC 

C.90OCE  + CC 

C.3CE5E-C  1 

C.25C9E-0L 

r.P.  - C.'.)73flC'<EOO 
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